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1. — Introduction. 


The motion of a fermion in various types of external fields is described by 
‘a four-component wave function y, satisfying the Dirac equation 


H, is the free particle Hamiltonian, 
H, = pm + ap, 


1 - Supplemento al Nuovo Cimento. 
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where f and a are the well known Dirac matrices, and units are chosen so 
that h—c—1. The representation of 8 and a is the usual one, with f diag- 
onal: f,,= Be, P= pa—— 1 he interaction with the external field 
is represented by the term gH,, where g is the coupling constant, for example 
the electric charge of the fermion in the case of external electromagnetic fields. 

Following FoLpy and WoUTHUYSEN (1), we distinguish between even and 
odd operators. An odd operator in the Dirac theory is a Dirac matrix which 


has only matrix elements connecting upper and lower components of the wave © | 


function, while an even operator is one having no such matrix elements. It 
can easily be shown that a necessary and sufficient condition for a matrix to 
be even (odd) is that it commutes (anticommutes) with B. This allows one 
to write any matrix, say M, as a sum of an even and an odd matrix in a 
simple way: 

MMM, 
where 

M°-}:(M+MB)= even part, 


M°=}3(M—-Mf)= odd part. 


If the Hamiltonian were an even operator, the Dirac equation would sepa- 
rate into 2 two-component wave equations, one equation for the upper com- 
ponents, and one for the lower components. In their paper FoLpy and 
WoUTHUYSEN demonstrate how the Hamiltonian, by means of successive 
canonical transformations, can be converted to an operator, which is even to 
any desired order in the expansion parameter 1/m. The Hamiltonian is put 
in the form 


H=fm+6+0, 
where € is even and @ is odd. Then a canonical transformation is performed, 


yp = exp [iSy], 


7 : . 0 ha ne 
HH’ = exp [as| (4 a ai exp [— iS] + eres 


3 SCIE DE : 
=H+i|3,a—i5|-5]s, past) 
With S=—:80/2m, the lowest order odd term in 1/m will be of first order, 


i.e. of one order higher than in the original Hamiltonian. By a sequence of 
further transformations, the generator of the transformation at each step being 


(1) L. L. FoLpry and S. WourHuyrsEN: Phys. Rev., 78, 29 (1950). 
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CANONICAL TRANSFORMATIONS OF DIRAC’S EQUATION TO EVEN FORMS ETC, 3 
chosen as S =— (i8/2m)x (odd terms in the Hamiltonian of lowest order in 
1/m), the Hamiltonian can be made even to any desired order in 1/m. The 


resulting transformation is the product 


ee. =... exp [45,] exp (78,1, 
S,=— (i/2m)pO, S, = — (i/4m?) E é—i 5; , 
and the transformed Hamiltonian to third order in 1/m is 
H* = fm + 6 + (1/2m)B 02 — (1/8m?) o, È é — 5] —_ 
= (1/89) 86, IA =| - (1/8m*)BOr+ .... 


When the particle is free, the transformation can be written in closed form, 


U, = exp[i8], 


where S is odd and Hermitian, 


w 


S = — (i/2p)Bap tg (p/m) , BS 


3 | 


The transformed Hamiltonian becomes 


& = BE, = pV m2 p?. 
A eanonical transformation in closed form, making the Hamiltonian even, 


by time-independent external fields, is given by ERIKSEN (?). He considers 
the equation 


| (6, H")=[8, UHU*|=0, 


which is a necessary and sufficient condition for H" to be even. The con- 
dition is equivalent to 


[ORB A | == 0. 
A possible choice for U*8U is 


(2) U*pU = H|/VH*=A, 


(2) E. Eriksen: Phys. Rev., 111, 1011 (1958). 
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1.€. 


B CET H*/V(H*)? : 
or 


He = VHS). 


This means that the transformation makes the Hamiltonian even and energy- 

separating. That is: the wave equation for the upper components of yp" has 

positive energy eigenvalues only, and the wave equation for the lower com- 

ponents has negative energy eigenvalues only. Obviously, the condition (2) 

is also a necessary one, for the Hamiltonian to be even and energy-separating. 
As a solution of eq. (2) ERIKSEN takes 


U=VBA. 
This is a square root of a unitary operator: 
UZ=%yzUuz uu 


For the definition of this symbol we refer to the Appendix. 
U has the property (see Appendix), 


pu = U's, 
which implies that U can be written as \ 
U= exp[i8], 
where S is odd and Hermitian. This shows that U, in the general case, is 
different from U,w given by (1). According to a theorem by HAUSSDORFF (*) 

for the product of «symbolic » exponential functions, 
exp [A] exp[B] = exp [A+ B+ 3[A, B]+ higher order commutators] , 
we have to third order in 1/m 
Uyw = exp [iS + i8, + iS; + 4[S8,, So] +...] . 


[S,, Ss] is an even operator, and generally non-vanishing. Hence, the trans- 
formations are in general different. They also produce different Hamiltonians. 


(?) F. Haussporrr: Ber. Stichs. Ges., 58, 19 (1906). 
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To lowest order in 1/m the difference is (1/16m?)B[€,{[, ©?]], which should 
be added to Hfy in order to give Eriksen’s result. When the particle is free, 
the transformations are equal. 
When the interaction term is odd, d.e. 
H=fPm+ 0, 
the transformation with 4 gives a simple result: 


H"=BVm+ 0. 


Applying this to an electron in a magnetostatic field, H = curl A, we get 


6) H*(® = 0) = B[m°+ (p— e4)°:— eo], 


which result has also been obtained by CASE (*), by a different approach. 
When, however, an external electric field is present, there is no simple 


- form for the transformed Hamiltonian. Neither is there any simple expansion 


in powers of the electric charge. In general, if one wants to have H” expanded 
in powers of some coupling parameter, one will find extremely complicated 
results already in the first order term. This is due to the fact that the zero 
order term of 6A (the square root of which is U) is no c-number quantity, but 
equal to the operator (m+fap)EH, °. 

The main purpose of the present paper is to introduce various types of 
canonical transformations, which, apart from making the Hamiltonian even 
and energy separating, have the property that they can be expanded in powers 
of the coupling constant as simply as possible. 

Since the condition U*B8U= A is a necessary and sufficient one in order 
to have the Hamiltonian even and energy-separating, we shall first. consider 
the operator A, and expand it in powers of the coupling constant g. An expan- 
sion of 2 in powers of 1/m will also be made. 

Then we shall present three unitary operators in closed form, satisfying 
the equation U*BU=, and which are appropriate for an expansion in powers 


_ of g. Two of them are in the free particle case different from the Foldy-Wout- 


huysen transformation function, and in some respect of a simpler nature. 

Transformations in powers of g will also be found by an iteration proce- 
dure analogous to the Foldy-Wouthuysen method. In this connection we 
derive advantage from a pleasant feature of the iteration method, namely 
that the transformed Hamiltonian to order 2n+1 is being expressed in terms 
of the » first generators of the transformation. 


(4) K. M. Case: Phys. Rev., 95, 1323 (1954). 
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Then we shall proceed to the general case with explicitly time dependent 
external fields, and it will be shown that the generalization mainly consists 
in replacing 2 by an operator A, which is a constant of motion in the Dirac 
picture. 

The transformed Hamiltonians to first order in e will be written out for 
time dependent external electromagnetic fields. Certain terms of first order, 
which can be expressed as commutators with BH,, can be removed to first 
order by an additional even unitary transformation. The transition to the 
classical case is performed. As an example, the Born approximation in one 
of the new representations (pictures) is calculated. 

In the new representations the Dirac particle shows, in a way, a behaviour 
characteristic of a particle with extension. The charge distribution of this 
extended particle, in the various representations, is expressed in terms of well 
known functions. 


A) Time Independent External Fields. 


2. — Expansions of the operator 1. 
21. Hxrpansions in powers of g. — We shall make the series expansion 


A= H/VH®= A+ g4%+9g°4,+..., 


where 
H = H,+ gH,. 
We obviously have 


(4) ea 
(5) [A, H,+ gH,|=0, 
(6) Ao = H o 


As will be seen below these properties are sufficient to determine the expansion 
coefficients ,. 

In the following we will make use of the fact that any operator F can be 
written as the sum of two parts, one part (F°) commuting with 2,, and one 
part (°°) anticommuting with 4): 

eee 
(7) F°=HF+ FA), 


eee 


f 
4 È rl 
Se ee 
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According to (4) 


Il 


(Ap + 9A + ga, + ...)>= 1, 


(8) Aghy oi Ado =0, 
| MA, de AAA 0% “eto. 


By means of these relations 4}, can be expressed in terms of the /,’s of lower 
~ order 


ZA =0, 
Aa =—4A,4,, etc. 
or generally 
n—-—1 
(9) Xi, SS to > Ax Link 
A k=1 


The anticommuting parts are determined by (5), which is equivalent to 
[Ags tol > [Anca Ha) = 0 ; m= 1, 2 
The anticommuting part of this equation is 
pA Se Anna Hi =. 
Putting H,=/,H, and multiplying by /,, we get 


(10) (4, B,], = bs 


nm? n? 


da ca Aol An—1s H,| ° 


This equation has one solution only, since the homogeneous equation is only 
satisfied by Z$= 0 (EH, being positive definite). 
It is possible to write 2; in the following way: 


0 
(11) As ex Lior lb, exp LE ri dr. 


—o 


This way of writing operators has also been used by MoRPURGO (°) (eq. (51) 


(®) G. Morpurco: Nuovo Cimento, 15, 624 (1960). 
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in his paper). The proof is simple: 


0 


[Ares ales Lee = exp [E,7|(H,b% + bh E,) exp [H,t] dt = 
cs: 0 
= | a(exp [E,t]br exp [E,7])= bi. 


Hence, 2, can be expressed as follows 


0 
1 ni al 
(12) An ET 9 o > Ax nt SE DI exp [E, T]I[Ao, [AS H,]] exp [E,t] dt 7 
sd k=1 


where the sum is zero when n= 1: 


(13) ev | exp [2,7] [o Hy] exp (H,c] dr. 


—co 


For the evaluation of these integrals it may be convenient to make the fol- 
lowing transformation: 


r'=exp{E,t]rexp[-Erl=r+qt[B,rlt = 


=r-+it[H,, 0/0p)=r—icE;p,; 


exp [E,t] F(r) = exp[E,t] F(r) exp [— E,t] exp [#,t] = F(r’) exp [E,7] . 


Hence 
0 
A= ia[. H,(r’)] exp [2H,t] dr . 
Example a). Electron in a constant electric field E, i.e. H,=—Er: 
1 NE Meni © Cho 
Mes A, (r )] = [Aes — Er a it E; Ep| = — Eldo, rie iE i 3 
0 
x Ody. f È 1: 0À, 
(14) A, = AE a [exp [2E,r]drt = ; id, Ha? de = 


= 5 Bs (imbaE — o( px E)). 


i 
e 
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Example b). Fourier component: H,(k) = exp[ikr]. 


exp [E,t] exp [ikr] = exp [ikr] exp [—tkr] exp [E,t] exp [ikr] = 

| = exp [ikr] exp [E,t], E, =[m?+(p +k)*}. 
Here we have used 

exp[— ikr]p exp [ikr]= p+ k. 
Hence 


0 


3 I(k) = fox [E,t] exp [ikr] exp [B,T]dT = 


= 0 


— exp [ékr] | exp [(E! + E,)t]dr = exp [ikr](25+ 2), 


A, (k) == (Kc) = Ay I(k) Ay a 
= exp [ikr]{1 — exp[— ikr] /, exp [ikr]A,}(Z,+ E) = 
= exp [ikr]1—- XA)(E+E;}!, 
where J, = (H,+ak)/E,. 


Example c). Oscillator potential: H, = r°: 
exp[E,t][%, rr] exp[H,t] = [Ao, (r —itH,‘p) exp [2,1] (r + itH, 'p)| 


= [ Ao, r exp [2H,t|r —it[H, ‘p exp [24,7], r]] 
= /4,,7r exp(2H.7Ir| , 


A= tale rexp[2E,r]r]dt = DIA rag tte a Be 


2°2. Haxpansion in powers of 1/m. — In order to find the coefficients 4,,, in 
the series 


A= ho + (A/m)Aay + (A/m)?Ag) + è 


we use the same technique as above: 


(16) nel 
(17) [A, HP=0, H=mB+4h=mB+€+0, 
in (18) ho = 6. 


= 
a 
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By means of (16) the even part of 7, can be expressed in terms of the coef- 
ficients of lower order 


n-1 


agree 1 e 
An Ta +p > Ace) Ain ’ ia = 0. 
k=1 


The odd part is determined by (17), which gives the following equations 


[Ano PI+ [Amay, AJ = 0, n=1,2,.. 


1.0. 
Kin = $B Ams, h] . 
Hence 
n-1 
(19) Any ae Mey a5 Kin) a 3B([Za-v; h] a > Aix) Mn) È 
k=1 


The lowest order terms are 
hay = 0,. hey i 3B[0, é| a BO? ’ 
As) = t[¢, [é, ©] a i[0, LO, cal FEO, 
We are now in a position to make a formal proof that the Foldy-Wouthuysen 
transformation gives an energy separating representation in the case of non 


explicitly time dependent external fields. Consider the operator UfyfUrw(=Z). 
We obviously have 


VAR = il and Z => b + (1/m) Za + DISCO 
Further, we know that Hfw is even, i.e. 


[B, Ur HU = 0 ’ 
or 


[Z,H |= 0. 
Thus Z satisfies the equations which define 2, i.e. 
L= User U eee =A ’ 


but this is just the condition which U must satisfy in order to give an energy 
separating representation. 
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3. — Transformations in closed form. 


As mentioned in the introductien, a necessary and sufficient condition for 
the transformed Hamiltonian, H"= UHU*, to be even and energy separat- 
ing is 


(20) UAUt = B, where UU* = U*U=1. 


A solution of this equation is the operator U= Vf, which, however, gives 
extremely complicated terms when expanded in powers of g. 

Our problem is to find a transformation which brings 4 over into 6. Our 
device is that this can be effected in two (or more) steps, by introducing an 
intermediate operator f, f =f* =, so that 4 transforms into f according 
to the scheme 


A>f>f. 


The first transformation, -f, can be carried out by means of the trans- 
formation function Vf, 


VAAAVAE =f, 
and the second transformation, f >, correspondingly by 
VBI IVIB=B- 


The resulting transformation, being the product of two unitary transfor- 
mations, is of course also unitary: 


(21) OSV piv iA; 
UAU* = VBf VfAaAv af Vip = VBitV {B= B. 


Obviously the method can be generalized to an arbitrary number of trans- 
formations, 


U = VBfaViati Vi Vira « 
Now consider the definition (Appendix) of Vf}, 
(22) VA = (1+fA)@2+f4 +20). 


A simple expansion in powers of g (when expressed in terms of 4;) can be 
obtained by requiring the zero order term in the square root to be a c-number 


tt 
a 
oa 


a 
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quantity. Letting f be of zero order, we eae 

(23) fay + def = c-number , Dpef E dat 
31. Transformation I. — As is well known there exists a Dirac matrix, 


(24) n= By; = — LAZIONE 


with the properties 
y= =7, [ns pl=0, (7, ele = 0- 


Hence this matrix is a possible choice for f in eq. (23), where the c-number 
quantity then becomes zero, 


[n, dol =[n, Bm +ap],H,* = 0. 
The corresponding transformation in two steps is 
(25) Ur = Vn Vn, 


which transforms 2 into 6 via the Dirac matrix 7. Due to the relation 
Bn+nP=0, the square root ./f7 is simply 


— 1 4 
VBI = Jyh + Ay) 


== 


Vi (1 + ips) Vd. 


To first order in the coupling constant, 

(26) U, = (1 + tys)(1 + A.) [1 + g(4[m, 41] + 44041) +]. 

In the limit of vanishing external fields, 

(27) U, > Vo=4(1 + tys)(1 + 74) = $(1 + tys) Hy (EB, + Ho) - 


It should be noticed that this transformation does not coincide with the Foldy- 
Wouthuysen result, which can be written as 


(28) Uy = E,(E, a m)|* (E, ci BHo) 0 


The transformed Hamiltonian is in either case equal to BE,. 
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In order to make calculations and expressions simpler, we rather express 
U, in the following way: 


U1 = VB VndNnB V Bn = VBA VEN « 


Here, and in the following, we let a dash on a symbol denote the symbol trans- 
formed by £7, 


(29) F'= Vpn FV np =4(1 + ys) FA — dv), 
i.e. 


ny=B, B=—n, H=-nmtop. 
Then the transformed result can be interpreted as 
ta VBI A VIB, © H= Vin vn. 


The transformation from.H to H' means going over to a Dirac matrix repre- 


sentation whose characteristic property is to make the free particle Hamil- 


tonian an odd operator. The next transformation, H '—> H”, is in fact Frik- 


sen’s transformation, U, =, in this representation. 
In this notation, and to first order in the coupling constant, we get | 


1 
V2 
(31) AY = BE,+4g(A[0, Hy), + BLA, Br). + blr, Hol) + + 


(30) VB = Ty (1 + BM) + 23 0 (BAY — 4M + BA); 


As usual the indices e and o denote even and odd part, respectively. 
There are two special cases in which the Hamiltonian Hy is considerably 


simplified : 


a) H, is even. This is equivalent to 
[H,,6]=0, ie. [Hin] = (Ah, tBys] = 0. 


According to (13), 4, is given by 
0 
M= is exp [H,t] [Z6, Hi] exp [#7] dr, 


= 60 


. : . VSRTO 
which is an even operator, since 4, is odd. Hence 4, =0, and we have 


(32) H = BE, + 49(H, + AHA) +. 


999 
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b) H, is odd. Then the Hamiltonian H' is odd, . 
[D',8]:=0, te ([H, yh =([H, By.) =0. 
Thus, the characteristic feature in this case is the fact that there exists a Dirac 
matrix which anticommutes with the Hamiltonian. Since A’ is odd, we get 
— ; 1 ; 
VBA'= (1 + Ba')(2 + BA + AB) = St BI), 
(33) Hi 4(1 + BAa') H’ (1— #2')=}1+ BA’) N°, 
BA'H' = BVH", 


where H'= }(1+-1y,) H(1 — ty). 
Consider the case of a fermion with an anomalous magnetic moment x in 
a magnetostatic field H = curl A, 


= pot dr 


My is the Dirac part = e/2m, and 4, is the intrinsic part of the magnetic mo- 
ment. Then one can use a Hamiltonian amplified with a Pauli interaction 
term, 


H=fm+a(p—eA)—y;feH, 
which is an example of a Hamiltonian which anticommutes with èYs, 
H' =— nm + a(p — eA) + uw, noH, 
H\' = B[m? + (p — eA)? — 2muoH + u,po(H x p— px H)+ WH}. 


When yu, = 0, this is the same result as given by eq. (3), which was obtained 
by a different transformation. 

Because of the term 3g8[4/°, 7], the transformed Hamiltonian to first 
order in the general case (31) is by no means any simple operator. One should, 
however, notice that the term can be written as 


bgBl Ay’, Hy, ag ESZAZIO PE) ’ 
which is possible due to the relations [4,, 6], = 0 and [4% a\],=0. Thus, 


to first order the term can be removed from the Hamiltonian by means of an 
additional even unitary transformation 


@ = exp [— 494A] , 


1.00 


ee cee ee ENN gi wee : det Pn 
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giving the Hamiltonian 
tro tr 104! ti 
GH7G = Hi — 3glA hg By]+.... 


To first order we therefore have 


(34) Hy > BE, + 494, [40 Hy), + 39804, HY, +. 


3°2. Transformation II. — We can make a different simple choice for f in 
eq. (23), in order to obtain a c-number for the expression fo + Zof. If we put 


{= Ae; 1.6. fi + Aff = 2, 
we get a two-step transformation 


(35) Un = VB VAÀ; 


which transforms 4 into f via the operater 4. The transition Z,->f is in 
fact a free particle F.W. transformation: 


VBA, >, Un = [2E,(E, + m)} (E, “in BH) . 


In the case of vanishing external fields we obviously have 


In the general (stationary) case we obtain the following expressions for the 
transformation function and the transformed Hamiltonian, to third order: 
(36) Oy, = U1+ 49404 + = 92(SAg As + Ay) +7 9° (640%, — 2A Ay + Ag Al) + +) è 
(37) Hy = PE, + U(9Hî+ tol, Ail, + 39%%[4, Hil, +...) Uo 


The index a (or c) means anticommuting (or commuting) part with respect 
to A,, as explained in eq. (7). 


3°3. Transformation III. — This transformation is carried out in three steps, 
A+>h->n >, 
that is 


Um = V By VÀ Via À ; 
(38) or 


Um — Vi VA cd , 
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where V is the free particle transformation function which was found as the 
limiting case of U,, eq. (27), de. 


Wo= 314 tys5)(1 tr No) . 


Due to the formal similarity of U,,, and U,, the following expansions can be 
written down by replacing U, with V,in eqs. (36) and (37) respectively: 


(39) Uy, = Vo(1+ 49444 £9? (4Aq Az + A) + 9? (640 Ag — 2A do + Agar) + --) 5 


(40) HY, = BE, + Vo(gH; + 19°%[4, Hi], + 49° Ag, Hi], + ---) Vo - 


In order to simplify calculations, we put U,,, in a more practical form 
(as was done for U,): 


Un = V Bn Vaà V dod = V Bn Vado VA VB V Bn = VBA VAA VBA - | 


The dash denotes transformation with Vf7, as explained in eq. (29). Then 
the transformed Hamiltonian can be written as 


(41) = VBMVANH'VXIRVAB, B= Ven HV of. 


The transition from H to H’ corresponds to a change in the Dirac-matrix rep- 


‘resentation, as mentioned earlier, and going from H' to H" is in fact trans- 


formation II applied to the Hamiltonian H’. 
From (41)—and remembering (8)—we then get to first order in g 


(42) His = BE, + 39 VBAodo[do, Hil V U8 , 
where 
Ay = (— nm + ap)E,* , 
Hy = ¥(1 + iy;)H:(1— ts) 
By means of the relation 
V Bho = FA GAM, 


we can now show that to first order H*, is equal to Hy’, given by (34): 


Hy, = BE, + 49(1 + BAYAA05 Ail, (1 — BA,) 
ae = PE, + 4g(1+ PADA(1+ BAy)[Ao, Hy], + 
| | + 49(1+ BA)A,(1 — BA)(A, HT} 


| = BE, + 39443 Hy), + 39004, Hy], - 
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4. — Transformations by iteration. 


In this section we shall proceed along the following lines: 


a) By an iteration procedure, analogous to the F.W. method, the Hamil- 
tonian can be transformed to an operator, THT*, which commutes with 4 
to any desired order in the coupling constant, 


(44) [THT*, A] = 0. 


b) Then a free particle transformation is made, which transforms o 
into S. This can be effected either by the operator U,, eq. (28), or by the 
operator V,, eq. (27): 


UU hi, Aa == 
Then we obtain two different Hamiltonians, both of them being even 
(45) HY, = UTHT*Uî, 
(46) Fhe ce GA A a eh i 


The notation used is due to the fact that these transformations are related 
to the transformations (35) and (38), respectively: 


Un U.Va, Une Vo apd: 
The even character of the transformed Hamiltonians is easily demon- 
strated by 
[Hy , 6] =[UTHT*U,, Ud Us| = U[THT*, 4A]U0 = 0, 


and quite analogously for H¥,.. 


Iterative solution of the equation |THT™*, 2,|= 0. In the present theory the 
role played by 4 is analogous to the role played by f in the F.W. theory. 
Instead of writing the Hamiltonian as the sum of an even and odd part, we 
now write it as the sum of a c-type part (commuting with /,) and an a-type 
part (anticommuting with /,): 


dl. Hy gHy= H+ 9H, + gHi, 


H; = 3H Ay Hy Ao) ’ 
Hi = 3(H,— AH; Ay) ° 


2 - Supplemento al Nuovo Cimento. 
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In order to remove the a-type terms of first order, we make a canonical trans- 


formation, generated by the operator gs,, $1 being an Hermitian a-type oper- | 
ator 


(47) [s1, dole = 0. 
Then, 
Hs H® = exp [igs,|\H exp [— igs;| 
= Hy + gH; + 9H; + tg[s, Ho+ gHi + 94s) +... 


We require the first order a-type part of this expression to vanish: 


(48) i[s,, Hol +Hi=0. 
With the substitution 


8 = 40x40; 
we get the following equations for x,: 
(o Hel [es HJ =O, Tos Aeon 
These are just the same equations as those which defined /, in a unique way: 
[Ai Hel iA l= oO, [Age Agi Oe 
Hence the only solution of (47) and (48) is 
(49) ee Ty he 
To third order the transformed Hamiltonian generated by this operator is 
(50) H® = H,+ gH; + gold, Hil, + 39404, Hi] — 
— 19° [A [A HSI] — &9° [4 [A Wi] +. 
The next transformation is one which removes the term 39*4[A,, Hi]. Re- 
quiring s, to be an a-type operator (corresponding to being odd in the F.W. 
theory), one can show that the only possible choice is to put 


53 = 4454, 
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In general, the n’th generator s, shall satisfy the equations 


Sn a Sn 3 
(51) [Sn, hol+ = 0 ’ 
[8n, Ho] = thy, ’ 


where A, is some Hermitian a-type operator which depends on the preceding 
transformations. Multiplying the last equation from the right by A, we get 


[8n, Hy], = ihn Ag ’ 


whose only solution, according to (11), is 
0 
Sn =] exp [E,r]thnA, exp [H,t] dt . 
Thus the iteration procedure outlined above can be carried out to any order, 
‘and the restriction, [s,, A)|, = 0, makes the expansion unique. 

We should like to point out a pleasant feature of the iteration method. 
Let us assume that n transformations have been carried out. Then the result, 
H™, is a c-type operator to order n. By the next transformation, generated 
by the operator g"+!8,,,, the following expression is added to the Hamiltonian 


H+0 — AM — tg" S47, HH] E - 3 gan+a [en wage Her] dl ae 


The lowest order c-type term in this expression is of order 2n+2. That is, 
up to order 2n+1 the c-type part of the Hamiltonian is not affected by this 
transformation. Obviously, this part is neither affected by the succeeding 
transformations. Hence, to order 2n+1 the transformed Hamiltonian is equal 
to the c-type part of H”. The transformed result to third order is therefore 
the c-type part of H™, given by eq. (50), 


THT* = H,+ gH, + t97A[A1, Hil,— j9° [Ax [4,; Hî]] rosta 
According to eq. (45) and (46) we then have the following even Hamiltonians 
(52) Hy, = BE, + Ui(gHy +49°4[4,, Hi), — 39° (A, (4, Hil] +...) To , 

(53)  Ayy= BE, + Vo(gH{ + 49°A[A,, Hil, —F9° [A [4 Wil] + -.)Vo » 
where 


U, = [2E,(E,+ m)|-*(Z,+ BH) and V,=43(1+ iy.) ED (E, +nH)- 
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To second order in the coupling constant these results coincide with the results 
from transformations II and III respectively. 

In contradiction to the transformations found in a closed form, (I, II 
and III), the iterative solutions are found without making use of the equation 
UAU* =f, which is the condition for the transformation to be energy sepa- 
rating. We shall show that this relation is valid for U,, and U,,,, as far as 
the expansions in powers of g are convergent. We make use of the property 


Ao, PELL [= 0” where Lr] [ig?s,] exp [igs] , 


which is equivalent to 
[T*A,T,H]|= 0. 
Obviously we also have 


(T*AT) =1 and (I AoE Varo order — ho si 


These are however the same equations as (4), (5) and (6), which defined the 
expansion of 4 in powers of g. Hence 


ToT =A, 
or 


DATE == fle 


Making a unitary transformation of this relation, by means of U, and Vo 
respectively, we get: 


UT Ut = UA, 05 5 VITAE 


Uy AU Cei B ’ U sir AU rr mi B . 


B) Generalization to Time Dependent External Fields. 


5. — Notation and introduction. 


The Dirac equation will now be written as 


(54) Ky = (Kt ghi)y = 0, 
where 
(55) K,= Hin, iGo eee - 
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Making the transformation y*= Uy, the Hamiltonian is converted to 


Hu gru + it. 


The necessary and sufficient condition for this expression to be an even oper- 
ator is 


[B; UKU*] =0 ’ 
Te. 
(56) [U* BU, x) =| U%pU, B43] =o, 


Hence, the operator U*$U has to be a constant of motion in the Dirac picture. 

In the case of stationary external fields, the operator U*BU was put equal 
to A= H/VH?. This is, however, not possible in the present case, since in 
general 


What we want is a generalization of 4 to an operator A, with the following 
properties 


(57) iA=[A, Ky =|4, H-iz|=0, 

(58) JA =A? Se 

Further, in the case of vanishing external fields we require 
(59) A>h= H,E;’, i.e. A=À- 
Then U can be found as a solution of the equation 


(60) U' BUS e UAUS= $,, 


which is quite analogous to eq. (20). 


6. — Expansions of the operator /. 


61. Expansion in powers of g. - We want to determine the operator coef- 
ficients A; in the expansion 


ei gr; 
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where 

Ay == ho 
Requiring the equations 

DS 5 


(4. Hi 4+gH)| = 0, 


to be satisfied to each order in g, we obtain the following equations, 


(61) ¥ Ap Ans =, n= ioe 
k=0 

: a 

(62) [4 hes iz) + [A,a, i) =0, n= 


By means of (61) the c-type part of A, (the part commuting with 4) can be 
expressed in terms of the A,’s of lower order: 


Ay =—44,43, 


; n-1 
Ae =—3A, Si Asie. . 


k=1 


Next, let us consider the a-type part of eq. (62): 


aw i! 
[at Ho— 55] =— 34140 Cea Fall 


Multiplication with A, gives the equation 


RAI 
(63) [An, E,= ilo tap = Bi ’ 


where B* is the explicitly time dependent operator 
B;(t) = $[Ao, [An-(t), H,(t)]] : 


The following solution has been found 


0 


(64) A, = = (1+ A,) | exp [2,7] B3(t —iz) exp [H,1] dr + Herm. conj. 


_o 
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The proof goes as follows: 


OAn 
Ei; El}: iA, — ot = 
0 


3(1 + Ay) } exp [E,t]{[Bs(t— ir), Hy], — i OBy(t — ix)/6t}- 


— 0 


‘exp [E,r]dt + Herm. conj. = 


0 
41+ A) fexp [E,t]}{E,Bs(t— it) + OBR (t— ir)/dv + Bi(t— ir) E,}- 


—o 


‘exp[H,t]dr + Herm. conj. = 


0 
4(1+ 4) afexp [E,t]Ba(t— it) exp [E,7]} + Herm. conj. = 


—-o 


= 3(1 + Ay) Ba(t) + 3 Br(t)(1+A4) = Bild). 


I 


. When the operator B* is explicitly independent of #, then this will also 
be true for the solution which we have found for A; (*). In fact, 


0 
An >3(1+ Ay) I exp [E,t] Bi exp [E,t]dr + Herm. conj. = 
= o 
2 exp [H,t] Bi exp [H,r]dr , 


-—o 


where 


Bes | ey Agata) 


Comparing this expression with the previous result, 


0 
(11) Ai = |exp [E,1t]b;, exp [E,r] dr, 


—_ 00 


be = 4[A, hae 19 H,]], 


we see that for stationary external fields we obviously have A= 4. 


(*) This property of A; will be spoilt if one adds solutions of the homogenous 
equation. 
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In the general case the operator A, can be written as 


A, = $(Ap — 1) ¥7 + (Ao + DE 


where 
0 si 
IF=|[exp[E,t][Ao, Hilt + 17)| exp [E,7]dr. 


00: 


We shall give two examples: 


a) H,=—Erf(t), where E is a constant and f(t) a c-number function. 
i) 
#2 = [exp (HytI(%a, — Erlf(t + tr) exp [B,r]dr = 


0 
= /exp [H,7] 74 ES f(t + iz) exp[E,t]dT = 


== 00; 


= (4K - e (ta + ir) exp [2E,t]dr. 


Putting for instance f(t)= exp[iwt], we get 


0 
I= =i exp [iot] E ch (2H, —- ww)", 
op 
A, = exp [iwt]H,'(4E% — #2) 1(2imPaE — 20(p x E) + iokE; Ep — iwaE) . 
When w — 0, 
A, > /, = 3H,’ (impaE—o(p~x E)) , 
which is the previous result, eq. (14). 


b) We consider a Fourier component H,(k) = exp[ikr —iot]. By means 
of previous calculations, leading to eq. (15), we have 


 exp[E,t][4, exp [ikr]| = exp [tkr](A, — 4) exp [E;7]; 
where 


E,=[m?+(p+k)?]}, = (Hy. +ak)/E,. 
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Hence 
0 


IF =| exp [H,t][A), exp [ikr]] exp [— iwt + wt] exp [E,T] dr = 


55 0 


= exp [ikr — iwt](A6— in) exp (FE) = 


— co 


= exp [tkr — iwt|(A.— 4,)(#, +H, +o), 
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A; = $ exp [ikr — it] {1 —4,)(B, +B, + o)*+ (1+4,)(B, +B, — 0) } (1 — Ad) - 


When w +0, 


A, > A, = exp [ikr](1 —Aj/,)(H,, + E;)-, 


which is the previous result, eq. (15). 


62. Expansion in powers of 1/m. — We want to find the coefficients 4,, 


- in the expansion, 


AEM = le) ATE, 


when A is an operator satisfying the equations 


een 

(65) “e SE 
’ ot Pe ’ 

Ao=f; 


H=mB+h=mB+6+ 0. 
The commutator equation can be written as 


rar mp + h‘] MZ mB + C+ 0]= 0, 


where 
6'= E — = 
= ve 
Then (65) is formally equivalent to: 
(16) X=1, 
(17) [A, mb+6+0]=0, 


(19) 92 ho = B È 
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Hence we can make use of the previous result (19): 


n= 


Am = $B([Aw-», h'] i >. Aa Ann) * 
k=1 


Since the operator 0/0t will always appear in commutators, the operators Any 
will not contain differential operators in #. 
As a curiosity we would like to mention that the operator, 


KjVE? = (x Ma). 


where i(0/dt) is being considered as being of zero order, gives an expansion in 
powers of 1/m which is just the one we have found. 


7. — Transformations in closed form. 


As has been shown before, eq. (60), U can be found as a solution of the 
equation 


(60) UPU=A 

which is quite analogous to 

(20) | U*pu =A. 

Obviously, the previous transformations can be generalized to 
Us = 77 pAy 
Ur =vVBbynvnA, 


Un = V BAW AAs 
Un = Vpn VA, VA,A - 
Introducing the expansion 
A=A,+g9A,+..., 
the transformations I, II and III can be expanded in powers of g. The results 


are quite analogous to the previous ones, given by eq. (26), (36) and (39), 
respectively. The only modification being necessary is to replace A; by Aj. 
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For the transformed Hamiltonian we find: 


, ! U e I o 10 I 0Ay 
(66 ) Ay= BE, 39 (4014, Hy la + BA, di if + BLA, oh iA at) + ’ 


(67) Hy, = BE, + U.(gH; + 3g°A4[A1, Hi] + 4g*A4[ Az, Hil +...) Uo - 
(68) Hi,= BE, + Vo(gH° + 49°AslA,, Hil + tg°Aol 43, Hi}, + ---)Vo, 


where as usual the dash denotes the transformation F’= VnBF vVBn. 

By means of an additional even transformation of Hy, generated by the 
operator S,= (i/2)gAi°A,, the terms containing A," can be removed to first 
order: 


exp [iS,] (2; — i 4) exp [— iS] = Hi a dis 39 Aj’ Ao, BE,— 


2 
ot ot tai 
fd OAL 
di ara 
gue ai GAS 


= Hi - to; a 5 Go da 2 op[AP, Hole 


As— g9IApA BAH] +... = 


That is, 
(66) Hi ud BE» La 392% H 1 RL z5 4 9BLAy, Hi ie o 
The results, given by (66), (67) and (68), are to first order in g equal to 


corresponding transformation results for time-independent external fields, given 
by eq. (34), (37) and (40), respectively. 


8. — Transformations by iteration. 


We shall briefly discuss how the iteration solution in powers of g can be 
modified in the case of time-dependent external fields. 
The Dirac equation is written as 


Ky = (K.+ gk, + 9Ki)y = 90, 
where 


Ke H,— += R=, Kr=e=Hi. 
By successive transformations we want to determine a total transformation 7’: 


y>Ty, 
* 
H > THT*— iT 


AG 
=TET*+ ig 
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such that the Hamiltonian is converted to a c-type operator 
* . 0 

GRICE + are Ale AU , 
ET A: E=08 


In order to remove the a-type terms of first order in K, we make a trans- 
formation generated by an operator g8,, 


Si == Se , [So Al — 0 . 
Then, 
K -> KW = exp [tgS,|K exp [—igS,] = 
= Ko gKi ie gk ig[ Si, Ko de gKi Si gKi] DE DCO 
Requiring 


181, Rol Ki = 0, 


and making the substitution 


i A 
S, = 5 XA, 4 
di 


we get the following equations, 


‘ Co) 
x, Hy i 5, PAS] Ole [SX Apia OF 


which are satisfied by X,= A,. Hence 


The corresponding operator in the time independent case was 


(49) 


To third order the operator S, generates the result: 
K> KY = Ky + 9K + tg°AlA,, Ki, + $9?Ao[ Ai, Ky] — 


— 39° [A FA, Kal] = A; (41, He] J, S 
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The next transformation is generated by the operator (i/2)g?A}A,, which re- 
moves the a-type term }g?A,[A,, Ki]. 

According to a previous statement, the transformed A to third order in g 
is equal to the c-type part of K°: 


TKT* = K,+gE + i@A[A4, Ki]. —49° [A (Av, Kel] + 
pee n rr i 
: gti toro 
== (HO 1a gHi+ 1gA[A, Hi 39° [An [A,, Hsl] Ped 


This Hamiltonian, which commutes with A, can now be transformed by 
either of the free particle transformations, U, or VW, to an even operator: 


(69) HY, = BE, + U,(gHi + 49?A lA, Hi — 39° [Ar [Ar, Hi] + i 
(70) Hi = BE, o V(94H° + igrA[4, Hi a 59° [A,, [A, Hal + DOG 


To first order in g, these results are exactly equal to those which were found 
by iteration methods in the case of time-independent external fields, eq. (52) 
and (53). 


C) Electromagnetic External Fields. 


9. — Dirac particle in electromagnetic field. 


91. Transformed Hamiltonians to first order in e. — We consider the Hamil- 
tonian for a Dirac particle with electric charge e in a general time-dependent 
electromagnetic field, 


H = bBm+ap+e(6—aA) = H,+ eH, . 


According to eq. (67) and (69), the transformations II and Il' give the same 
Hamiltonian to first order in e: 


H® = pH, + eU,H:Uf, U,=vBi=[2E,(E,+ m)|*(E,+ BHo) 
By means of the relations, Uh =U, and A,U; = UP, we have, 


UH? Us = $U,(Hy + Ao Hy 4g) US = 4U, HU; + $8 UH. Us B = (Ur Hi U,)° = 
= (e/2) [E,(E,+m)]*{(E,+ m) P(E, + m) + (ap) D(ap)— 
— B(ap)(aA) (EB, + m)— (Hy + m)(aA)(op)B} [Hy (E, + my" « 
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This gives the following Hamiltonian: 


‘ +5 | P P 


il a di 
Ss pren: Mea di er 


P P m 
elia) +m 


where 


S=}0 and F=—VD+i[E,, A]. 


According to eq. (68) and (70), the transformations III and III’ also give 
equal Hamiltonians to first order in e: 


H*, = BE, + eW,H°Vi - 
As was demonstrated in eq. (43), this can be expressed as 
Hix = BE + 3edo[A Hy] + 3eB[A,, Hr hh, 


which is equal to the Hamiltonian HY given by eq. (34), when this has been 
modified by an additional transformation. Introducing 


H=3}1+iy,)(D—aA)1—iy)=D_ad, 
we easily find: 
(72) Htm= BE, + 4e(® + AOA) — teB(MaA + aA A) = 
< BOE na {E,DE, + (op + iBm) Oop — ifm) — 
= BE, + E DI + Poe a Sn Pee px y= 


i me 5, 1 1 e(p 
edi z — pS 5 lg VxA tx Ag) BS (F A+A FY, 


where 


S=}c, and F=_-VO+i[fE,, A]. 
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We now want to show that the quantity i[8£,, A] in F can be replaced 
by —0A/0t by means of simple, even, unitary transformations. Then the oper- 
ator F is replaced by the total electric field, E=—V@ — 04/0. In general 
we put 7 


A 


H™ = exp [ies] (we — i) exp [— ies] + ii = H* + te È HY — 15 +..., 


where s is some zero order even and Hermitian operator. The first order equa- 
tion is 

Tr, tr 9 os 
(H h= (H + t[s, BE,] — aE 

Hence, if some part of (H*), can be written as i[8Z,,s], then this part will 
be replaced by — 0s/0t. We easily see that we can choose 


i È, 
s=—S ws (A x p— p XA) i in eq. (71) 
and 
pS) A A + 2mpBA : in e (72) 
8s = on, | x P— Px + m E |. (4 


After these additional transformations we now get the following Hamiltonians 
for a fermion in a general time-dependent electromagnetic field, to first or- 
der in e: 


È ef. m P P 
(73) E Dt 


»p+m E, +m 
Baw Exp PXE) za 8 (H+ ali 
i IR 
° (74) T= BB, + ORSAY: 


ira (Ex p— PXE)z + PS ES — 
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Hi is even and time-independent: 
H©.= U, HU = BE, +49(H, + HA) + > 


e.g. for a fermion in an electrostatic field, Hij=:D'2@ 


= Ra (047 of + Lor) 95 2x P) 7 PSE OT. 


92. Low and high energy limits. — At low energies (p< m) and to second 
order in 1/m we get 


2 e ali 
(75) Tr Bm 4 - - eb 3 S3 (EX p— px E) 
© sH-8- (pA4+ Ap) 
ue, Ali + Ap), 
(76) Him = pm+ 2+ @— “85 (Ex p— px E) 


€ € e 
+8 SE> fi SH— fi (pan see 


To the order considered, expression (75) coincides with Foldy and Wouthuysen’s 
result. Expression (76), however, contains twice the spin orbit term, i.e. without 
the Thomas-Frenkel correction. On the other hand there appears a static 
coupling to the electric field, B(e/m)SE. If this term is considered individually 
it suggests an electric dipole moment (-operator) equal to — f(e/m)S, i.e. equal 
to the magnetic moment, but of opposite sign. The so-called Darwin terms, 
(e/3m2)V®D and (e/4m?)V2®, are to be discussed later, in connection with the 
charge distribution. 

At high energies (p>m) all Hamiltonians contain the same leading terms 
to first order in e: 


Po vr e 

(77) Him Bp o (042 P oF p(Ba+ a2) 7 
2° \p Pp 

The characteristic feature is the dependence of the scalar potential on the 

direction of the momentum. 


9°3. Transition to the classical case. — This is accomplished by making all 
dynamical variables commute. Even though that corresponds to putting 
hk = 0, we keep the spin variable S, just in order to see how it enters into the 
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classical expressions. We also put 6=1 for positive energies. From eq. (73) 
and (74) we then obtain the following classical Hamiltonians, to first order 
me: 


el Selle ei, @ i Eat e 
fr i tmp Hg 

m é e 
Ho = E —L eb — pi SEX p) + e qa SE— 7 SH— 7 Ap. 


We will also express the energy, being equal to the Hamiltonian, by the par- 
ticle velocity v. Introducing conventional units, and dropping the term eA 
in p, (second order!), we get to first order in e 


È ME e 1 e e e 
a y 400-5 8(L Exe) + % SE y o SH VA: 


where 


9°4. Born approximation in one of the new representations. — As has been 
pointed out, the main feature of the present transformations is that the wave 
equation splits in 2 two-component equations, corresponding to positive and 
negative energies. Let us now consider a Dirac particle with positive energy 
scattered by an electrostatic field ©. Choosing transformation II, we get 
from (73) the wave equation 


|P @ P LS = - pxE 1 Pol? 
dica E,+m ' E,tmP Em \it+gete 


where g is the two-component wave function. 
As can easily be verified, this can also be written as 


. op PE Bis um 
ue) Lot, DPR (og Var pu 


which form is convenient for the calculations below. We introduce free par- 


3 - Supplemento al Nuovo Cimento. 


1019 


per 


34 E. ERIKSEN and M. KOLSRUD 
ticle states which are eigenstates of H,, p and o,=op/p: 


COS — 


2 
|#, pit) = 6 exp [ip’r] , Op = 1, 
sin = 
2 
= cals 
|E, p'’,-—> = 6 exp [ip’r], op=—1. 
COs 5) 


E= (m?+ p"”)*; 0 is the angle between p’ and the z-axis (in the xz-plane). It 
should be noted that to zero order in e: p"= p and = Oops 

These states are simpler and more «natural » than in Dirac’s representation. 
On the other hand the transformed Hamiltonians are in general far more com- 
plicated, but seem—in return—to lend themselves more easily to physical 
interpretations than the mathematically simpler Dirac Hamiltonian does. 

Let us calculate the matrix element of (H77),, as given by eq. (78), between 
an initial state with momentum p,= pk along the ¢-axis, and a final state 
with momentum p;= p(k cos0+i sin 0) in the «2-plane. 

As E, and op are quantized in both cases, we obtain 


Etktm 


dal m _ id 
<E, Pn |HT|E,pt>= 3(1 i SCE O ®;; gea\l- sin 5) = 


i m i 6 0 
<E, pr,+|HT|E, pi, >= 5(1 DIC ie & Di AE = Dx C085, 


0 
=— 0,7 sins. 


By forming the sum of the squares, i.e. averaging over the spins along p;, we 
get the first order Born approximation for the relativistic scattering cross- 
section (apart from a factor), 


2 2 
o(0) ~|@ ,|? (cose ae a sin? 3) =D} f _ = sin? 5) : 


10. — Charge distributions of the fermion in the new representations. 


In their paper Fotpy and WouTHUYSEN obtain the following transformed 
Hamiltonian to second order in 1/m (ref. (1), eq. (35), corrected) : 


H"= Am + eb4+ È (p — eA)? — p= 0H- 


é 
Sm? 


E ot, Are Me 
o(E x (p— eA) — (p— eA)XxE) spa VET 
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The last term is a well known correction to the Pauli theory, due to ©. G. 
DARWIN (1928). This term—« which has previously been regarded as of rather 
mysterious origin, can now easily be understood, since it comes from the fact 
that the electric charge in the new representation is spread out over a finite 
region »—of the order of its Compton wave length 4 (= h/me) (*). 

In order to see that the potential energy of an extended particle gives rise 
to a term of this kind, we consider a spherically symmetric charge distribution 
0(s) about the position r. The potential energy of this distribution in the field 
P(r) is given by 


Vir) — feta) ® r+ s)d?s = x a fara o(s)(sV)" D(r), 


n=0 


which is easily found to be 


(82) VD(r) = eD + =<) VOR 
e<s®% is the 20th charge moment, given by 


ecs) ala 0(s) d's 
Hence, in the F.W. case 


<s? > ew = 34° 


In order to obtain more detailed knowledge of the charge distribution, 
one may continue the iteration process, selecting terms of the type eV'!2, 
eVSD, ete. Then, according to (79), higher order charge moments can be 
determined. To fourth order, for instance, we get 


Vow (r) = eD+ 1eX?V?D+ TE e7iViD. 


These calculations are rapidly becoming very tedious for increasing order of 
1/m, and therefore greatly limit the possibilities of knowing the F.W, charge 
distribution accurately. 

We notice that the energy terms of interest are proportional to the charge e, 
consequently they are all contained in the first order term of the transformed 
Hamiltonian, when this is expanded in powers of e. In the F.W. theory, however, 
the transformed Hamiltonian does not lend itself easily to an expansion of 
this kind. 


(*) F. and W. connect this feature with their new « mean position operator ». 
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Let us rather examine the charge distributions of the fermion in the repre- 
sentations which are found in this paper. Then we want to pick out from the 
first order transformed Hamiltonian the part which depends on ®, and which 
is free from the spin variables. These terms are, according to the expres- 
sions (74) and (73), respectively, 


ie Vee be 
(80) Wim(r, p) = 3 (2 a E, Arp a E, oF) ; 


(81) Wulr, p) = wise (+ aa? i Da DICE +5 


In order to select the co-ordinate function part of these operators, we apply 
them to a state with zero momentum, e.g. a constant. The result will then be 
the function which is interpreted as the potential energy V(r) of an extended 


particle: 


Wir, p\l.p=®=Vr)|p=®. 


In this manner W,m and W,,, from eq. (80) and (81), give respectively 


(82) Vil) => (1 ~ =) P(r) = ePsn (2) CALA: 


Pp 


(83) Var) =e V s(t ca 4 O(r) = ePu (2) P(r) . 


The operators P(p/m) are to operate upon @(r) only. 

Incidentally, from these expressions we may obtain the two sets of charge 
moments, simply by expanding the operators P in powers of p/m=—iZV. 
We get in the respective representations 


Vim(r) = e® + +e4°V?0 Pi 3 CAVID + .. 
V(r) = eD +4 eA2V2@ + th ehV'® aes 
As to the transformation U, = VÀ, it does not give H™ in powers of ¢ 
in a simply way, and hence we do not know W(r, p). It is, however, pos- 


sible to find the result of (H™), applied upon a constant, and from that expres- 
sion we may pick out the term we are seeking. Putting A= 0, we get in fact (*) 


(Hes PS 0) = V,(r)*|p= 0), 


(*) The details will be published elsewhere. 
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where 


ne 
2 (E,/m)— i P(r) = e Pa( p/m) D(r) 


(84) Va(r) = e NI 1+ (m/E,) (E)m+1 


| 
: 


Expansion in powers of p/m yields 


V(r) = eB + Leh2V2® + mR RVIDÒ.. 


To find the analytical form of the various charge distributions corresponding 


to these potential energies, we proceed in the following manner 


V(r) = eP(p/m) Dr) = PUY) | Or") d(r’— r) der’ = 


= i P(r’) P(AV)(222)- “foca exp [i(r’— r)k] = 


a = «| 3p D(r 2a) | ae PR exp [i(r’— r)k] = 


= ars Dr + s) o(s), 
where o(s) is interpreted as the charge density, given by 


sap | POO exp [iSK]|d?K , 


Si sid. 


Looking at the expressions for P(p/m) in eq. (82), (83) and (84), we see 
that P(co) #0. This means that o is partially a 6-function, i.e. a point charge: 


T Op Gan c0)] exp [tSK] dK , 


The Fourier integrals above can be expressed in terms of well known func- 


tions. We obtain the following charge distributions for the fermion 
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In picture I and III: 


IL 


Oxi 2 6(s) + de Es — K,(S8). 


In picture IT: 


ye A! 
gu = 7a + gr § dI 8 


[FQ)W_2a28) — TOM:C9]. 


In picture E: 


On = (V2— 1)e 6(s) + = [116 — fit) dl rs 


+ STOW 1208) — Temes 


27 = 
K, is the modified Hankel function, and W,,, is Whittaker’s function. 

In addition to the 6-function singularities, these charge densities have 
singularities 1/s? in s= 0. 

The asymptotic forms when s — oo, are: 


—6/4 


exp [—s/A], 


Crim A0S 
gi 


Cras exp [— 8/4] ’ 


—9/4 


Cn vs exp[_s/4]. 


APPENDIX 


Square root of unitary operator. 
Let uw be a unitary operator, uu* = u*u = 1. Then the operator 
Vea (lena 4 449 
is a solution of the equation X? = w. This expression is well defined, provided 


2+u +u* is a positive definite operator. Let us examine whether this is true. 
The eigenvalues of w are on the unit circle in the complex plane: 


| u|d) = exp[id]]ò), 
VA 


u*|ò) = exp[— id]|d). 
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Hence: 
(2+u+u*)|d) = 2(1 + cos 6)|d>. 


This eigenvalue is positive definite, except when cos é=—1. Thus, the op- 


erator (2+u-+u*)-? (positive root) is well defined when — 1 is not an eigen- 
value of u. In the following we assume this to be true. Then, 


Jue =(1 + uy (2+ u+uryt= (1+ Qu + w\(2+ ut Uy, 
= u(u* + 2+ u\(2+u+u*)*=u, q.e.d. 
The hermitian conjugate of V% is 
(Vu)* = A+) (2 + wt + ut Ve: 
| As is easily seen y/% is unitary: 
va(va)t= (Vu) Vu = 1. 
Of special interest is the case: 
ries Di, 


| where P= F* — F-1, G = G* = @" (for instance, F=f, G= A). 
Then, 


VEG 
(VFG)* = VGF = (14 GF)(2 + GF + FG)*. 


= (1+ FG)\(2+FO4 GF), 


By means of these expression one can easily verify the following relations: 
FVFG=VGFF, 
VGFFVFG=dG. 

| Putting F=f and G= 7, we see that the operator U=~VBA satisfies the 


| relations 
' BU=U*B, U*BU=A, CoS OU = 1; 


which have been used in the text. 
When F and G anticommute, the square root is a very simple expression: 


A = al 
VG == sO FG): 
Vv 2 
Putting F=f and G=ipy;, we get 
== yar a 
VBiBy; = Viys===(1+ 5) - 
Bipy sae y 
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DEL NUOVO CIMENTO 4° Trimestre 
Statistical Weights in Multiple Particle Production (*). ; 
È 


J. SHAPIRO 
Fordham University - New York, N.Y. 


(ricevuto il 28 Luglio 1960) 


4 ee SAR een orm 


Contents. — 1. Introduction. — 2. Statement of the problem. — 3. Deri- 
vation of the formula. — 4. Tables. 


— rt 


1. — Introduction. 


statistical weights enter because of the usual assumption of conservation of 
total isotopic spin. These weights are the product of two factors. One of these 
factors is the number of independent degenerate states (due to particle ex- 
changes) for a final state described by giving the number of particles of 
different types and charges. The other is the probability, for a specific state, 
that a measurement of the square of the total isotopic spin of the system will 
yield the value corresponding to the initial two particle state. i 
The purpose of this paper is to derive a simple analytic formula for these ; 
weights, valid for arbitrary numbers of particles in the final state, provided 1 
that they have isotopic spins 0, 3, or 1. Extensive tables are given for the — 
second factor for those reactions which are experimentally of interest. From I 


In treating the process of multiple production in high energy collisions, | 
i 
{ 


these tables, a simple computation gives the statistical weights (equation (5)). 
With regard to notation, it should be noted that upper case letters 7, 7's, 
will be used for isotopic spin operators and lower case letters, t, m, for the cor- 
responding quantum numbers. If a symbol is to refer to only a part of the 
system of particles, a subscript will be appended. | 


(*) This work was supported in part by the U.S. Atomic Energy Commission. The 
computer at Fordham is supported by the National Science Foundation. 
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2. — Statement of the problem. 
The reactions being considered are of the type 
A+B+ many particles . 


The cross-section for such a reaction is proportional to the square of the tran- 
sition matrix element between the initial and final states 


(1) aol F If. 


The initial state is either a pure ¢ state (e.g. pp collision) or a simple mixture 
of such states with known probabilities (e.g. a px® system is in either a t= 3 
state with a relative probability 3 or in a t= 3 state with a relative proba- 

bility of 3). Hence, the treatment is restricted to the case of an initial pure 
i state. If necessary, one can average over such states at the end. The initial 
state, of course, is always a pure m state. 

The final state consists of many non-interacting particles. Consequently, 
it has a product isotopic spin wave function with a definite value of m, but 
not of ¢. (Effects due to antisymmetrization of the fermions will not be con- 
sidered in the present paper.) This wave function can be expanded in terms 
of normalized simultaneous eigenfunctions of 7 and T;, i.e. 


(2) lp im = >t, my. 


t 


Tf one assumes that # is conserved (charge independence), then the reaction 
ean go only to that term in the expansion corresponding to the ¢ of the initial 
state. Hence, from (1) and (2) 


(3) GOG = dams 


A formula for the quantities P,, has been obtained by CERULUS (1) in terms 
of Clebsch-Gordon coefficients. For certain cases he has also obtained closed 
formulae for the P,,, (2). A simple closed formula will be derived in this 


(1) F. Cerutus: Suppl. Nuovo Cimento, 15, 402 (1960). Poe m, Of this reference 
is the same as P,,,(«, f) here. +P _.m, of CERULUS is the same as Cie 

(2) F. CerULUS (preprint): This paper gives a general procedure for obtaining 
closed formulae which yields results that, in appearance, differ appreciably from those 
given here. He obtains explicit formulae for the case of n pions, n pions and 1 nucleon, 
and for n pions and 2 nucleons. The author is indebted to Dr. CeRULUS for a pre- 
publication copy of this paper. 
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paper, which is valid for any number of particles in the final state, provided 
only that they have t,= 0, 4, or 1. 

The quantity P,,, refers to a single product wave function. Experimentally, 
however, one is usually interested in a final state described by giving the 
number of particles of various types and charges, (e.9. p+n+2K*+2K +n+4-3n° 
from a p+p reaction). 

There are many independent product wave functions for such a state, cor- 
responding to particle exchanges. Consider the particles as separated into 
classes corresponding to the different isotopic spin multiplets, i.e. nucleons, 
antinucleons, pions, ete. Let n(t;) be the number of particles in class è and 
n(t;, m;) be the number of particles in class è with a specific value m, for T3,, 
so that 


(4) mt) = > n(t;, mi) - 


Then, under the assumption of charge independence, one obtains degenerate 
independent wave functions under any exchange of two particles in the same 
class with different values for m,. Hence, the complete statistical weight for 
a state described as above is given by 


[n(t,)]! 
Gin Uh, m,)| = TI => = Pipl (ti, mi)}, 


i 


| i ([n(ti, mi)] o) 


mj= bt; 


— 
DI 
— 


where G is a function of all the quantities n. Tables have been given in the 
literature for P,, and G,, for limited numbers of particles (**). These tables 
have been calculated by using known expressions for the Clebsch-Gordon 
coefficients, and combining the angular momenta two at a time. Extensive 
tables of P,,, are given at the end of this paper. The reason for tabulating Ws 
rather than G, is that (as will be shown) it depends on only t, m and 2 other 
integers, whereas G depends on the complete set of numbers n/(t;, m,). 
Consider the various reactions leading to the same multiplicity in each 
class, but differing in the charges of some of the particles (e.g. in a pp colli- 


(3) R. H. MiLBURN: Rev. Mod. Phys., 27, 1 (1955); E .FERMI: Phys. Rev., 92, 452 
(1953); 98, 1434 (1954). There are a few errors in the tables in these papers (see 
Tables IV and V of MiLBurNn). The corrected values, with the authors’ values given 
in parentheses, are as follows: final state (n+ + ——) Gs , _; = 80/35 (87/35), Gy, 4 = 
= 2(123/60); final state (n+00—), Gg, y= 136/35 (129/35), Gy, _x= 14/5 (165/60); final 
state (pp+0—), G,;= 12/5 (154/60); final state (pn+ + —), G,,,= 3(175/60) and 
final state (pn+00), G,,,= 21/10 (121/60). 
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sion leading to 2 nucleons and 3 pions, the possibilities ppr'r°m, ppr?n’n?, 
porn, pnrtrtm, and nnntxtnx®. If one makes the statistical hypothesis 
that the final states (i.e. independent wave functions) for these reactions have 
equal weight, then the quantities G,,, are equal to the branching ratios between 
the above processes. The sum of the statistical weights G,, over the above 
set of reactions equals the total number of independent isotopic spin wave 
functions which can be constructed leading to the specified values of ¢ and m. 
Alternatively, this sum is the total number of multiplets with the given value 
of t for a system with the specified number of particles in each class. Tables 


of these quantities have been published in the literature (4°). These can be 


obtained from explicit formulae (4) or from simple recurrence relations obtained 
from the rule for combination of angular momenta (°). 
3. — Derivation of the formula. 


A closed expression for P,,, will now be derived, making use of angular 
momentum projection operators. Since the result is invariant under a re- 


flection in the 1, 2 plane in isotopic sipn space, one can assume m > 0. The 


values for negative m are given by 
(6) a nl [n(t;, = m,) | = Plots m;)| . 
An alternative expansion of the final state wave function in terms of ¢ states is 


7) Pel = 0a 
t 


where 0,, is the projection operator giving a state described by ¢ and m from 
state described by m. Comparing with (2), one obtains 


(8) Gili M> = Oi, |M) + 
Since the projection operators satisfy the relation 
(9) O''m = Om) 

and since the states are normalized, it follows that 


(10) Pom = |C¢\? = (m|Om|M> . 


(4) Y. Yervrn and A. pE-SHALIT: Nuovo Cimento, 1, 1146 (L955). 
(®) V. S. BaraSenxoy and B. M. BarBaSEv: Suppl. Nuovo Cimento, 7, 19 (1957). 
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The explicit formula for the projection operators is (*°) 


zx i. 7?%—y(y+1) 
(11) Om = lig +1)— rv +1)’ 


where 7 is the maximum possible value of #, which here is just the sum of 
the t, values for the individual particles in the final state. It has been shown 
by LOwprn (8) that formula (11) can be rewritten as a sum involving the usual 
raising and lowering operators 


(12) Ow a 


(2¢ + 1)(¢+ m)! (1), POT 
(t— m)! seo vi(2t+r+1)! 


This is valid for m>0. For m<0, one has the above expression with m 
changed into —m and with the lowering operators to the right. Substituting 
(12) into (10) yields 


(13) Pa = ee i (1) Cm | Tm 
tm (t—m 0 v!i(26+v LAI 


The problem is now reduced to evaluating the diagonal matrix elements 
in equation (13), between product wave functions. The (t,, m,) values for the 
individual particles, will be assumed to be restricted to the 6 sets (1,1) (1,0), 
(1, —1), (2) 3), (4, —3) and (0,0). 

In equation (13) one can think of the operators as being expanded in terms 
of those for single particles, i.e. 


(14) Dora 


i 


Since raising a state of maximum projection annihilates it, and since the 
raising operators appear to the right in (13), it follows that P does not depend 
on any particles with m=t. Thus P depends on only the number of par- 
ticles of the types (1,0), (1, —1), (3, —}). Further, the wave function in 
isotopic spin space for a (1, —1) particle behaves, under rotations, exactly 
like the product wave function for two particles of type (4, — 4). Since one 
is considering product wave functions here, each (1, —1) particle is equivalent 


(°) P. A. M. Dirac: The Principles of Quantum Mechanics, 3rd ed. (Oxford), p. 32. 

(7) P. 0. LòwpIn: Phys. Rev., 97, 1509 (1955); Adv. Physics, 5, 1 (1956). 

(*) P. O. Léwpin: Angular Momenta Wave Functions Constructed by Projection 
Operators; Technical note from the Quantum Chemistry Group of Uppsala University, 
May 10, 1958. 
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to two (4, —4) particles. Thus, if one defines 
29 2 ’ 


a= n(1,0), 


B mi n(3, 3) adibito), 


(15) 


(here n(1, 0) means the total number of particles with t,= 1, m;= 0) then P 
for the actual particles in the final state is the same as for a set of « particles 
of type (1,0) and f-particles of type (1, —1). Henceforth, this equivalent 
set of particles will be used and Pim will be written as P,,(«, 6) to show that 
it depends on only « and f. 

Next, one expands 


(16) ToTtatloo 


where the operators on the right are the sums of the single particle operators 
for the x and f particle classes, respectively. Using the binomial theorem 


i 


and recalling that one has product wave functions and that each particle must 
be raised and lowered the same number of times, one obtains for the matrix 
elements in (13) the result 


(18) TATA AP > = 


Ses Vem rm BITE" BEB 


v 


where |x) and |f> are the product wave functions for the two types of par- 
ticles. 

Again, think of the operators in (18) as being expanded in terms of those 
for individual particles. Since, for both classes of particles, m;,=t,—1, one 
is concerned only with single raisings and lowerings (or no raisings, which just 
gives a factor unity). Since 


t.-t,+|1, 0) =2 4, 0), 
(19) 


tti4|4, i) ="; 4), 


it follows that to evaluate the matrix elements in (18) one has to merely count 
the number of possible ways of raising and lowering once, and multiply by 
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the appropriate power of 2 or 1. Consider the « matrix elements, assuming 
the operators expanded. One is raising and lowering 7 particles of type (1, 0), 
giving a factor 2°. There are + ways of selecting è out of « particles. For 


each combination one has i! possible permutations of the factors in 7°, and i! 
possible permutations of the factors in 7". Henee 


ila! 


(20) (a|Ti-Ta+]® = Hie 
Similarly 
(21) BV, eh (t—m-+yv—i)!f! 


(B—t+tm—v+i)! 


Combining (13), (18), (20) and (21), one obtains finally 


02) Pate p= MEDE Catimini A ja 


((_m v!i(2t+v+1)! t-mtv— 


The limits on the summation are such that the arguments of all factorials are 
greater than or equal to zero. Equation (22) can be written in a more sym- 
metrical form by the substitution 7=t—m-+y»—zi leading to 


(23) Pim(%, B) "= na 
$$ (=1)* [6 + 41 a " (8) 
Al sé A—t tm) it + a ie 


Here the limits of summation have been written explicitly. The lower limits 
can be written as zero, if desired, since the additional terms vanish because 
of the poles in (i+4—t+m)!. 

For large numbers of particles, (i.e. large « and f) equations (22) and (23) 
are inconvenient to use because of the large number of terms involved. A 
more useful formula can be derived as follows 

Recall the well known property of the hypergeometric functions of argu- 
ment unity, 


(1° + 1)(t +m)! 
m 


I'(c)I'(e—a—b) 
la learn =bp% 


(24) Fa, b, 6, 1 


provided 
Re (c—-a—b)> 0 c4a0,,—1L-2..% 


Considering (24) as a complex function of ec, one can readily show that for 
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p, n, k integers satisfying kK>n>p>0 


(k—n)... (k —p—n-+1) 
hii) .k=p +s) 


F( Pi N k, 1)= 


If k<p-+n, the right hand side of this expression vanishes. Using the poly- 
nomial expansion for the hypergeometric function, one readily obtains 


A (bp) (—1)?(k&—p)'\(k—n)! 
(25) sa te tem (k—p—n)!p!n! 


Setting p=n=tT—m, k=i+A+t—m, equation (25) becomes 


(26) : (—1)° "LG + A) ae : < (1) “n A+ s)! 
È (E +À-t+ mm)! mM (80° m_ a)! 


This is valid for i>0, j>0. For 0<itA<t—m, both sides of equation 
(26) vanish. Hence, this can be substituted into (23) and the lower limits 
on the sums over i and % taken as zero. One obtains 


Pim(&, B) = (21 +1)(t + m)1(t—m)!- 


ia; © we i VI Al a\{P DI 
2, 2d —m— it msi +i en) 3) 


In this, the sum over 2 can be performed using (24) yielding 


P(e, B) = (2t + 1)( + m)!(¢—m)!- 


wey 2 (—1)**(¢+m-+f6—s)!(i+s)! CARI 
Zeri m— Hitt m— Net mt RA 1 sali 


Here, the sum over i reduces to a hypergeometric function of argument 2, 
yielding finally 


(2t + 1)(t + sfizio rim (1'(0+m+f—39)! 


C+m+B+1)! 22 ani m —s)!(t--m—s)! 
-P( OL 8, Ve MA b Za 


(27) Pin(%, B) = 


Equation (22), (23) or (27) gives an analytic expression for the probabi- 
lity P,,(a, B). The statistical weight factor is then given by equation (5) 


(5) G tml ME: 5 m,)| = Il È ì Pim(xB) b) 
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which has been rewritten to emphasize that P,,, depends on only « and f, 
which are defined by equation (15). 
The above results are valid for m>0. For m< 0, if one defines 


(28) B'= n(3, 3) + 2n(1, 1) 


and as before, the probabilities are given by 


(29) Ps mi (% B')= Pi, Imi(% B= B) . 


Since the hypergeometric functions can be obtained easily by recurrence 
relations, (25) is by far the simplest form for computing the probabilities. 
Furthermore, in those reactions observed up to the present, one particle is 
always a nucleon, so that the maximum value of t is $. Thus, one is interested 
only in the two cases m=t and m=t—1. For these cases, (27) reduces to | 


DI 
(30) Pita, p) = 14 TR a, 1, 2t+2-+ 8, 2), 
(ai). Pula I 
2241 AB en 2t—1 | 
di Feo, 1, Shee 302) scs a, 2; 2t PASO 2) | 


Perhaps a few words are in order regarding possible particles with higher 
isotopic spins. Although the results given here are restricted to t;<1, they 
include higher spin particles provided they have maximum or minimum pro- 
jection (i.e. m;=t; or m;=— t;). This follows because in these cases the wave 
functions behave as product wave functions of spin + particles. For example, 
consider a (3, 3) particle. For m positive, it’ has no effect on P,,; for m nega- 
tive one simply increases f’ by 3. . 


4. — Tables. 


Tables have been computed for P,,,(«, 8) for the six sets of (t, m) values — 
(1,1), (1,0), (0, 0), (3, 3), (3, +) and (4, 4). These were computed as decimal 4 
fractions on the Bendix G-15 computer at Fordham University, using equa- | 
tions (30) and (31). Table I gives P,,(a, 6) for these six states for the range — 
of variables 0<a<9, 0<<18. Some tables of P,, as rational fractions | 
are given in reference (1). 

A simple example will illustrate the use of these tables. Consider the 
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TaBLE I. — Values of Pim, B). — Pim, B) for state (t, m). 
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B (ile, Uy) (1, 0) (0, 0) (8, 2) (3, 3) (3, 3) 
0 1.0 0 1.0 1.0 0 1.0 
1 15 5 C 8 33333 | .66667 
2 6 5 33333 | .66667 | .4 5 
3 5 a6 25 67143 | «4 4 
4 42857 | 4 2 so 38095 33333 
È 378. 35714 | .16667 | .44444 | .35714 | .28571 
SIC 33333 | .32143 | \14286 | 4 33333 | .25 
ber. 7 3 29167 | .125 36364 | .31111 | .22222 
8 27273 | .26667 | .11111 | .33333 ] 29091. | .2 
9 25 24545 lox. 30769 | .27273 181 82 
10 23077 | 22727 | .09091 | .28571 | .25641 166 67 
11 21429 | .21154 | .08333 | .26667 | .24176 153 85 
12 2 19780 | .07692 | .25 22857 142.86 
13 1875 18571 | .07143 | .23529 | .21667 13333 
14 17647 | - .175 06667 | .22222 | .20588 | .125 
15 16667 | .16544 | .0625 21053 | .19608 11765 
| 16 15789 | .15686 | .05882 | .2 18713 11111 
17 ST, 14912 | .05556 | -.19048 | .17895 | .10526 
i 08 14286 | .14211 | .05263 | .18182 | .17143 | .1 
0 5 0 0 6 666 67 33333 
1 45 Îh 16667 | .53333 | .46667 | .33333 
2 a È 16667 | carsl9 |a 3 
3 36714 | 36 (a 42857 | .36190 266 67 
4 32143 | .31429 | .13333 | .38889 | .33333 238 10 
2 29167 | .28571 | .11905 | .35556 | .30952 21429 
6 26667 | .26190 | .10714 | .32727 | .28889 194 44 
7 24545 | .24167 | .09722 | .30303 | .27071 17778 
8 22727 | .22424 | .08889 | .28205 | .25455 163 64 
9 21154 | .20909 | .08182 | .26374 | .24009 15152 
10 19780 | .19580 | .07576 | 24762 | .22711 14103 
11 18571 | 18407 | .07051 | .23333 | .21538 13187 
12 175 17363 | .06593 | .22059 | .20476 12381 
13 16544 | .16429 | .06190 | .20915~| .19510 11667 
14 15686 | .15588 | .05833 | .19883 | .18627 110 29 
15 14912 | .14828 | .05515 |~-.18947 | .17819 10458 
16 14211 | .14138 | .05229 | .18095 | .17076 099 42 
17 13571 | .13509 | .04971 | .17316 | .16391 09474 
18 12987 | -.12932 | .04737 | .16601 | .15758 090 48 
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TABLE I (continued). 
x | B OES) (1, 0) (0, 0) (8, 3) (3, 3) (3, 3) 
2 | 0 A 0 333 33 466 67 266 67 .333 33 
| 1 35 3 166 67 .419 05 333 33 .266 67 
| 2 31429 3 133 33 38095 323 81 .233 33 
| 3 28571 27857 116 67 .349 21 304 76 .209 52 
4 26190 25714 10476 322 22 28571 .190 48 
| 5 241 67 238 10 095 24 .29899 268 25 .174 60 
| | 6 924 24 221 43 08730 .27879 25253 16111 
| 7 209 09 206 82 08056 | .26107 23838 | .14949 
8 195 80 193 94 07475 .245 42 225 64 .139 39 
ae 18407 18252 069 70 .231 50 21412 .130 54 
(SEO .173 63 17233 065 27 .219 05 203 66 lea 
| 11 .164 29 163 19 06136 .207 84 19414 .11575 
| 12 .155 88 15494 057 88 .19771 18543 .109 52 
13 Lie 214828 147 48 05476 .18851 17745 .103 92 
14 .141 38 140 69 05196 .180 12 17011 .098 86 
| 15 .135 09 13448 049 43 .172 43 163 33 .094 26 
16 .129 32 128 79 04713 .165 37 157 06 .090 06 
17 .124 03 123 56 045 03 .158 86 15124 .086 22 
18 11914 11873 043 11 15283 145 83 .082 68 
(Ea MOTI O E 6 0 37143 | 4 D 
1a 027807 3 ni) .342 86 314 29 2 
| | 2° | 20714 25714 PI .317 46 28571 18571 
SEL 238 0 23571 092 86 .295 24 266 67 17143 
4..:| 22143 21905 08571 275 76 25079 15873 
5 | .20682 204 76 079 36 .258 59 236 80 147 62 
6 | .19394 19221 073 81 .243 36 224 24 137 88 
7.0 18252 18106 068 94 .229 77 21290 129 29 
8 .172 33 17110 064 65 .217 58 202 60 121 68 
942-103 19 16214 060 84 .206 59 193 21 11489 
10 15494 15405 057 44 .196 64 184 62 10879 
11 Wika ae 146 70 054 40 .187 58 176 73 103 30 
12 | .14069 14001 051 65 .17931 169 47 098 32 
13, |. 18448 133 89 049 16 a7 16276 093 79 
14 .128 79 128 28 046 90 164.75 156 55 089 65 
15 .123 56 12310 044 83 .158 30 15079 085 86 
16 11873 11833 04293 .152 34 145 42 082 37 
17 114.26 11390 04119 .146 81 140 42 07915 
18-5 4) 211017 109 79 039 58 .141 66 13574 07617 
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TABLE I (continued). 


a B (1,1) | 29 | (0,0 | 49 | 40 | GH 
Se ame .257 14 0 2 31429 .228 57 2 
| 1 .235 71 21429 n: .292 06 .25714 .171 43 
| 2 .219 05 21429 08571 .273 02 .247 62 15714 
3 .204 76 203 38 07857 .256 28 .234 92 .146 03 
4 .19221 19048 073 02 24141 .222 80 13651 
5 .181 06 179 65 068 25 .228 13 .211 69 .128 14 
6 .17110 16991 064 07 .216 18 .201 55 .12071 
7 .162 14 16113 06035 .205 39 .19230 .11406 
8 .15405 153 18 057 03 .195 60 .183 82 .108 09 
9 .14670 14595 05405 .186 68 .176 02 .10270 
10 .140 01 139 36 05135 .178 52 .168 84 .097 80 
11 .133 89 13332 048 90 .171 04 16221 .093 34 
12 .128 28 12777 046 67 16414 .156 06 .089 26 
| 13 .123 10 122 66 044 63 ABP TT .150 34 .085 52 
| 14 .118 33 11793 04276 .151 86 145 02 .082 07 
15 .113 90 11355 04103 .146 38 14006 .078 88 
16 .109 79 109 47 039 44 14128 .135 42 .075 93 
17 .105 96 105 67 03797 13651 .13106 .073 19 
18 .102 38 10213 036 60 .132 05 .126 98 .070 64 
Bilge.) .214 29 428 57 0 .269 84 28571 142 86 
1 .202 38 21429 07143 .253 97 .238 10 .142 86 
2 | .19048 19048 07143 .239 54 .222 22 134.92 
3 .179 65 17857 067 46 226.55 .210 68 .126 98 
4 .16991 168 83 063 49 .214 84 20058 11977 
5 | .16113 16017 059 88 .204 24 .19142 .113 28 
| 6 | .15318 15235 056 64 .194 61 .183 04 .107 42 
7 | .14595 145 23 05371 | .18581 | .17534 | .10212 
8 .139 36 13873 051 06 .17776 .168 23 .097 30 
Bim) | -133:32 13277 048 65 .170 37 .161 66 09291 
100: (0 127.77 127 28 046 45 .163 55 .155 57 .088 88 
11 | .122 66 12223 044 44 .157 24 .14991 .085 18 
aoe) St > 11793 11754 042 59 .151 40 .144 63 .08177 
13 .113 55 113 20 040 89 145 96 .139 70 .078 62 
14 .109 47 109 16 039 31 14090 .135 09 .075 70 
do 710567 105 40 037 85 .136 17 .13077 .072 98 
16 | .102 13 101 88 036 49 .13174 .126 71 .070 45 
dibdo 09881 098 58 035 22 .127 59 .122 89 .068 09 
Tigri: 159.095 70 095 49 034 04 .123 69 .119 29 .065 87 
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TABLE I (continued). 
x B (1, 1) (1, 0) (0, 0) (3, 3) (È 3) (3, 3) 
6 0 .190 48 0 142 86 238 10 .190 48 142 86 
1 .178 57 166 67 07143 225 11 .206 35 12698 
2 .168 83 166 67 063 49 213 56 .199 13 119 05 
3 .160 17 159 09 059 52 203 13 .190 48 11255 
4 .15235 15152 056 28 193 64 .182 26 106 78 
5 .145 23 144 52 053 39 18497 .174 66 10157 
6 .138 73 138 11 05078 17702 .167 63 096 82 
7 .132 77 132 23 048 41 169 71 6s 092 49 
8 .127 28 12681 046 24 16297 .155 09 08851 
9 12223 12180 044 26 156 73 .149 48 08486 
10 .117 54 ILA 042 43 15094 .144 24 08148 
11 .113 20 112 86 040 74 145 55 .139 36 078 36 
12 109 16 108 86 039 18 14053 .134 78 07547 
13 .105 40 105 12 037 73 135 84 .130 48 07278 
14 .101 88 101 63 036 39 13144 .126 45 070 26 
15 .098 58 098 36 035 13 127 32 .122 65 067 92 
16 .09549 -| .09529 033 96 123 44 .11907 065 72 
ili .092 58 092 40 032 86 11979 .115 69 063 66 
18 .089 85 089 68 03183 11635 .112 49 061 72 
7 0 .166 67 333 33 0 21212 .222 22 Halo 
1 .159 09 166 67 055 56 202 02 .19192 11111 
2 .151 52 15152 055 56 19270 .181 82 106 06 
33 1144 52 143 94 053 03 18415 .17405 10101 
4 .138 11 137 53 05051 17630 .167 06 096 35 
5 . 132/23 13170 048 17 169 08 .160 61 092 07 
6 .126 81 126 34 046 04 162 40 .154 62 088 15 
7 .121 80 121 39 044 08 156 22 ..149 06 084 54 
| 8 ALTI 116 80 042 27 150 49 .143 86 08120 
| 9 .112 86 112 53 040 60 145 15 .139 01 07811 
10 .108 86 108 56 039 06 140 16 .134 47 07524 
11 .105 12 104 85 037 62 13551 .130 20 07257 
ID .101 63 10139 036 29 13114 .126 19 070 08 
13 .098 36 098 14 035 04 12705 .122 42 067 75 
14 .095 29 095 09 033 88 123 20 .118 86 065 57 
15 .092 40 092 21 03279 119 57 .115 49 063 52 
16 .089 68 089 51 03176 116 14 .11231 061 60 
Ny .087 11 086 95 030 80 11291 .109 30 059 78 
18 .084 68 084 54 029 89 109 85 .106 44 058 07 
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TABLE I (continued). 
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ot B A) (1, 0) 
8 0 .15152 0 

i .143 94 13636 
2 .137 53 136 36 
3 .13170 13112 
4 126 34 125 87 
5 .12139 12098 
6 11680 | .11643 
n SRI2G30 ENTI 
8 .10856 | .10827 
9 .10485 | .10459 
10 10139 | .10115 
11 .09814 | .09792 
12 .095 09 | 094 89 
13 .09221 | .09203 
14 .089 51 089 34 
15 08695 | .08680 
16 .08454 | .08440 
17 .08225 | 08213 
18 .08009 | .07997 
0 13636 | .27273 
1 .13112 136 36 
2 .125 87 125 87 
3 .12098 120 63 
4 116 43 11608 
5 11221 11189 
6 .108 27 107 98 
7 10459 | .10433 
8 .10115 | .10091 

\ 9 .09792 | .09770 
10 .094 89 094 69 
11 .09203 | .09185 
12 08934 | 08918 
13 .08680 . .08665 
14 08440 | .08426 
15 .082 13 082 00 
16 .079 97 079 85 
17 M2. 4. OTT SA 
18 075 97 075 87 


0,0 | @® | &5 | 44 
ILIIL o! .19192%./, 16162) CITDII 
05556 | .18337 | .17172 | .10101 
05051 | .17560 | .16628 | .09596 
04798 | .16845 | .16006 | .09169 
04584 | .16185 | .15416 | .08780 
04390 | .15573 | .14864 | .08423 
04211 | .15004 | .14349 | .08092 
04046 | .14475 | .13867 | .07786 
03893 | .13981 | .13416 | .07502 
03751 | .13518 | .12992 | .07237 
03619 | .13085 | .12594 | .06990 
03495 | .12678 | .12218 | .06759 
03380 | .12295 | .11864 | .06542 
03271 | .11934 | .11530 | .06339 
03170 | .11594 | .11213 | .06148 
03074 | .11272 | .10913- | .05967 
02984 | .10967 | .10628 | .05797 
02898 | .10678 | .10358 | .05636 
02818 | .10404 10101 | .05484 

0 17483 | .18182 | .09091 
04545 | .16783 | .16084 | .09091 
04545 | .16131 | .15385 | .08741 
04371 | .15524 | .14825 | .08392 
04196 | .14961 | .14312 | .08065 
04033 | .14436 | .13834 | .07762 
03881 | .13945 | .13385 | .07480 
03740 | .13486 | .12964 | .07218 
03609 | .13056 | .12568 | .06973 
03486 | .12652 | .12195 | .06743 
03372 | .12271 | .11843 | ,06528 
03264 | .11912 | .11510 | .06326 
03163 | .11574 | .11195 | .06136 
03068 | .11253 | .10896 | .05956 
02978 | .10950 | .10613 | .05787 
02893 | .10663 | .10344 | .05627 
02813 | 10389 | .10087 | .05475 
02738 | .10130 | .09843 | .05331 
02666 | .09883 | .09611 | .05195 


eR) Pa has 
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reaction 
n+p—>n-4 3p+ 2p+ 3xnt + 279437. 
Here a= 2, B=3+2(3)=9. From the tables P,,(2, 9) = .182 52 
Po (25.9) = .069 70. 
Hence P,,(2, 9) = $(.182 52) + 4(.069 70) = .12611 and 


At ot St 
@ = gii Si 319191 (12611) = 2240(.12611) = 282.49. 


* OR OK 


The author would like to express his appreciation to Rev. Professor J. 
McConnEL for many: stimulating discussions and suggestions and also for 
checking the derivations in this paper. x i 


x 
SUPPLEMENTO AL VOLUME XVIII, SERIE X Nel, 1960 
DEL NUOVO CIMENTO i 4° Trimestre 


Dislocation Relaxation in Silver, Gold, Palladium and Platinum. 


P. G. BORDONI 


Istituto di Fisica- Matematica dell’ Università - Pisa 


M. Nuovo and L. VERDINI 


Istituto Nazionale di Ultracustica « O. M. Corbino » - Roma 
(ricevuto il 15 Novembre 1960) 
Contents. — 1. Introduction. — 2. Results — 3. Discussion. 1. Activation 


energy. 2. Relaxation spectra. 3. Frequency relaxation. 4. Influence 
of anneal and cold work. 5. Subsidiary peaks. — 4. Conelusions. 


1. — Introduction. 


The existence of a relaxation effect due to the dislocations in silver, gold, 
palladium and platinum has already been reported in a preliminary commu- 
nication to the Third International Congress of Acoustic (Stuttgart, 1959) (*). 
In the present paper a larger amount of experimental evidence concerning the 
same metals is analysed in order to obtain a more complete knowledge of the 
effect, as it was previously done with copper (?). 

In addition to the fundamental parameters (activation energy and char- 
acteristic time) whose values are discussed in the light of the present theories 
of dislocation motion, many other aspects of the effect have been considered, 
such as the width of the relaxation spectrum, the relation between the total 


(1) P. G. BorponI, M. Nuovo and L. VERDINI: Relaxation of dislocations in face- 
centered cubic metals, in Proc. of the Third International Congress of Acoustics (Stuttgart, 
September 1959), in press. 

(2) P. G. BorponI, M. Nuovo and L. Verpini: Nuovo Cimento, 14, 273 (1959). 
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number of dislocations and the characteristic time, the influence of anneal 
and cold work upon the total number of active dislocations and upon the width I 
of the relaxation spectrum. Much attention has also been given to the sub-. 
sidiary or Niblett and Wilks peaks, which seem to be a general if not veri 
evident feature of the relaxation effect. 

It may be hoped that the considerable amount of experimental evidence. 
which is now available shall provide a good basis for a more detailed model 
of dislocation motion, and for a better theoretical treatment of the relaxation 
effect. 


2. — Results. 


The measurements have been made on circular polycrystalline plates of 
chemically pure silver, gold, palladium and platinum whose main character- 
istics are listed in Table I together with their thermal and mechanical his- 
tory (*). The experimental technique is the same used in the previous in- 
vestigation on copper (*). 

The energy dissipation and the resonant frequency have been measured 
as a function of the temperature 7 in the range (60 ~ 300) °K for the two 
lowest flexural symmetrical modes of the plates, having respectively one and 
two nodal circles and no nodal diameters. On the same specimens the dissi- 
pation and the resonant frequency have also been measured for one or two 
longitudinal modes. These measurements have proved more difficult than the 
corresponding ones on copper, owing to the coarser grain structure of the 
present specimens. It has been found that the average grain size may be of 
the order of some millimeters when the material is machined from the cast. 
When the wavelength is of the same order of magnitude, the amplitude vs. 
frequency curve is no longer bell-shaped but exhibits side peaks and the meas- 
urement of dissipation is no longer possible. The above difficulty has been over- 
come by obtaining a finer grain structure whenever possible by cold-rolling and 
intermediate annealing (silver no. 3, Table I) and in every case by measuring 
the dissipation for longitudinal waves with wavelengths larger than 4.5 mm. 


(*) The authors wish to thank the « Ditta E. Longobardi - Metalli Preziosi» for 
its co-operation in the preparation of the samples and in their analysis. They also — 
express their gratitude to Prof. D. GeNTtILONI and to Mr. MACCARRONE of the Istituto | 
di Scienza delle Costruzioni of the University of Rome for their kind help in the 
mechanical treatment of the specimens. 

(4) See the paper quoted in footnote (2). Additional information on the experi- i 
mental apparatus may be found in the following papers: P. G. Borpon1: Nuovo Cimento, 
4, 177 (1947); Ric. Scient., 18, 103 (1948); P. G. Borponr and M. Nuovo: Acustica, 
7, 1 (1957); Ric. Scient., 27, 695 (1957); Acustica, 8, 351 (1958). 
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In platinum an intermediate vibration mode has also been employed for 
dissipation measurements. 


\ TaBLE I. — Characteristics of the specimens. 


Diam- | Thick- 
Thermal and mechanical treatments eter ness 
(mm) | (mm) 


Spec- 


Metal | Purity LARA 


1 (a) Cross rolled after casting from 12 mm | 36.0 10.00 
to 10 mm; machined. | 


2 (a) Cross rolled after casting from 12 mm | 29.0 5.00 
| to 10 mm; machined. | 


Silver 0.9980 (b) Cold rolled from 10 mm to 7.5 mm; 
annealed 1h at 900 °K; cold rolled 
from 7.5 to 5.0 mm; machined. 
3 (c) Annealed 3 h at 460 °K. 36.0 4.55 
| Annealed 4h at 493 °K. (È) 
3 | (e) Permanently strained 2.7% at room 
| temperature. 


& 


a) Machined from the cast. 

b) Annealed 7h at 502 °K. 25.0 
c) Permanently strained 3.8% at room 

| temperature. 


| Gold 0.9988| — | ( 


I 
=) 
bo 


(a) Machined from the cast and perma- 
nently strained 1% at room tem- 

Palladium | 0.9991) — | perature. 

Annealed 6 h at 423°K. 

| (c) Permanently strained 6.9% at room | 

temperature. 


bo 
ou 
bo 
DO) 
w 
ur 


— 
im) 
Se 


(a) Machined from the cast. 
| Annealed 4h at 500°K. 
Platinum |0.9989| — | (c) Annealed 11h at 693 °K. | 36.0 3.00 

| (A) Annealed 1h at 793 °K. | 


— 
a 
— 


(e) Permanently strained 3.2% at room | 
temperature. | 
| 


(*) With a ring of 44.0 mm external diameter and 0.6 mm thickness. 


à 


As it is shown by the lower diagram of Fig. 1, the dissipation peak in silver 
is quite well marked, whilst a much lower and flatter peak was found in a 
previous investigation on a not very pure specimen of the same mate- 


è 


1045. 


DS P. G. BORDONI, M. NUOVO and L. VERDINI 
rial (°). As it was observed in lead, copper and aluminium, the temperature 7, 


corresponding to the maximum dissipation increases with the vibration frequency 
fn of the specimen (5). The same peak is also clearly shown by the dissipation 


600-10- 


400 


a) 200 


50 75 100 125 150 
Temperature (°K) 


I 
50 75 100 125 150 


Temperature (°K) 


Fig. 1. — Dissipation peak and frequency relaxation in silver; curve A, f,,=18.58 kHz; 
B, jm=—50.55 kHz; ©) fn=:421.4 kHz; D, f,,—820.4 kHz. 


(*) P. G. Borpoxni: Journ. Acoust. Soc. Am., 26, 495 (1954); Ric. Scient., 23, 
1193 (1953). 

(*) As the resonant frequency changes with temperature, f,, is the value measured 
at the temperature 7,,. 
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measurements in gold (Fig. 2, lower diagram), whilst in palladium and plati- 
num the peak is somewhat smaller and flatter but still very evident (Fig. 3, 4 
lower diagrams). As in the case of silver, the temperature 7, is an increasing 
function of the frequency f,,. 


a) 


50 - 100 150 200 250 
dh Temperature (°K) 
300-105 
> Ai IRA Au 
si È; O 
Q J 
@) ™~ O 
200 we Sa 
Wy ù 
xf p O 
» @ 
100 
b) (A) 
| I 
0 
50 100 150 200 250 


Temperature (°K) 
Fig. 2. — Dissipation peak and frequency relaxation in gold; curve A, f= 15.99 kHz? 


B, f,,=55.60 kHz; CO, fn=561.0 kHz. 


The upper diagrams of Fig. 1-4 give the relative difference (f —f,)/f, be- 
tween the resonant frequency f measured at the temperature 7, and the re- 
lared value f, of the resonant frequency computed for the same temperature 


È 
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by means of a linear extrapolation of the room temperature slope (’). These 
differences have been computed only for the flexural modes, as the experimental 


60010 


40 


0 ; 100 200 300 400 
Temperature (°K) 


100 200 300 400 
Tem perature (°K) 


Fig. 3. — Dissipation peak and frequency relaxation in palladium; curve A, f,,= 24.92 kHz; 
B, fm=88.35 kHz; ©, f,,=—387.9 kHz. 


difficulties due to the grain structure of the specimens make the frequency 
measurements less accurate for longitudinal vibrations. In the upper diagrams 
of Fig. 1, 2 (silver and gold) an inflexion point is found at a temperature very 


(") See the paper quoted in footnote (2), p. 307. 
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near to that of the dissipation peak for the same vibration mode. This coinci- 
dence shows that the peaks are due to a relaxation effect. 


600 + 10- 
400 
ESTA 
fp 
a) 200 
0 a ni 
0 100 200 300 400 


Temperature (°K) 


=5 
15010 = 


fom 


100 


b) 


0 


0 100 200 300 400 
Temperature (°K) 


Fig. 4. — Dissipation peak and frequency relaxation in platinum; curve A, f,=9.6 kHz; 
B, fm=42.41 kHz; O, fm=111.8 kHz. 


The inflexion is not evident in the upper diagrams of Fig. 3, 4 (palladium 
and platinum) being masked by the pronounced curvature of the elastic mo- 
dulus vs. temperature diagram exhibited by the transition metals in this tempe- 
rature range (8). This curvature has an anomalous character as the concavity 


(8) J. A. Rayne: Phys. Rev., 118, 1545 (1960). 
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is turned towards the high values of the modulus, and it is not related to any 
relaxation effect but to the conduction electrons. However there can be no 
doubt that a frequency re- 
laxation takes place also for 
palladium and platinum in 
the same temperature inter- 
val in which the dissipation 
peak is observed. 

The experimental values 
of f, and 7, are listed in 
Table II together with the 
maximum energy dissipa- 
tion Q;*. When the loga- 
rithm of f, is plotted 
against 7, as it is done 
in Fig. 5, the points for 
each metal lie on a straight 
line; hence the relaxation 
effect is thermally activated, 
as it has already been found 
in the case of lead, copper 


and aluminium. 
Fig. 5. — Frequency dependence of the temper- stablish an exper- 
a q Pp p o establis exp 
ature 7,, of maximum dissipation. 


25-10 


be 
° 


k loge fm (eV degree ) 
an 


imental relation between 
the relaxation effect and 
the motion of dislocations, the influence of anneal and cold work on the dis- 
sipation peaks and on-frequency relaxation has been systematically investigated. 

It has been found that a moderate anneal decreases the height of the 
peak, as is shown for silver and platinum (Fig. 6, 9, lower diagrams) by 
the comparison of the curves A, obtained when the material was in the 
original state (°), and the curves B, taken after an annealing treatment (!°). 

If the temperature of the treatment is high enough, the peak may even- 
tually be cancelled as it is shown by the curve C of Fig. 6 (lower diagram) 
for silver and by the curves B of Fig. 7, 8 (lower diagrams) for gold and pal- 
ladium. In the case of platinum the height of the peak is also dependent upon 
the temperature of anneal, but some trace of it is still present after the spec- 
imen has been kept for 1 h at 793 °K (Fig. 9, lower diagram, curves C and D). 
It may be added that the effect of an high-temperature anneal is somewhat 


(9) That is immediately after the mechanical treatments listed in Table I under 
(a) for gold, palladium and platinum, and after (a) and (b) for silver. 
(!°) For the details of these treatments see Table I and the captions of the figures. 
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TaBLe II. — Experimental data on the relaration peak. 
| 
} 3 ital Mppolmen a nella | CE) mi 
ye | | 
1 a 28.61 | 73.3 129-10-6 | 
2 7 a 79.99 Dic 150 si 
Pu b 18.58 75.8 228 | 
ee i a ee 
b 421.4 91.0 288 | 
Silver 3 b 820.4 ~ "93.0 255 i 
€ o 55.28 nr 72.5 235 
| d RASTA ae Oe 
e | 55.35 | 80.0 184 | 
t a 15.99 di | 202 da 
+ a a 55.60 .| 127.7 i 291 
Gold — a 561.0 153.0 284 È 
DIG fee SRN A va 
tufo | 15.36 Tigo |}, 160 
È Bs a 102492 | 1890.0 30. 
a "88.35 | 186.5 59 7 
Palladium — i a 387.9 200.0 75 
b 9 as A fas 
‘ 21.69 | 187.0 125 > 
a È ay ere 9.600 id TO 
a ù 42.4] / 201.0 : 116 Li 
: a ia; 213.0 138 
Platinum — î b ; if 9.641 166.0 21.5 
e 9.652 166.0 10.0 ; 
d y 9.653 166.0 9.0 
e vo: 8.778 180.0 13% 


ri (*) For the details of thermal and mechanical treatments, see Table I. 
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masked in silver by the large increase of the dissipation measured at room 1 
temperature which has already been observed in copper after a thermal 


570 


550 


53.0 


ga 


Frequency (kHz) 


51.0 


0 50 100 150 200 250 300 
Temperature (°K) 


200 


0 50 100 150 200 250 300 
Temperature (°K) 


Fig. 6. — Influence of anneal on the height of the dissipation peak and on the frequency 
relaxation in silver. Curve A, after the treatments (a) and (5) of Table I; B, after 
3h at 460 °K; €, after 4h at 493 °K. 
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Fig. 7. — Influence of anneal on the height of the dissipation peak and on the fre- 
eroney relaxation in gold. Curve A, after the treatment (a) of Table I; B, after 
7h at 502 °K. 
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Fig. 8. — Influence of anneal on the height of the dissipation peak and on the fre- 
quency relaxation in palladium. Curve A, after the treatment (a) of Table I; B, after 
6h at 423 °K. 
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treatment (1) and which is not related to the relaxation effect. No such 
increase is noticed in the other metals, which seem therefore to be better 
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Fig. 9. — Influence of anneal on the height of the dissipation peak and on the fre- 
quency relaxation in platinum. Curve A, after the treatment (a) on Table I; B, after 
4h at 500°K; C, after 1l h at 693 °K; D, after 1h at 793 °K. 


(4) See the paper quoted in footnote (2), p. 292, or D. H. NisLert and J. Winks: 
Phil. Mag., 1, 415 (1956); 2, 1427 (1957); Adv. in Phys., 9, 1 (1960). 
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suited for an investigation on the influence of anneal upon the dissipation 
peak. 

All the annealing treatments were made in an evacuated container and 
they lasted long enough to allow an equilibrium state to be reached by the 
material. One of the resonant frequencies of the specimen was measured during 
the treatment, and the anneal was not stopped until the frequency was ap- 
proximately constant with time, as it is shown for silver and gold, in Fig. 10. 


1.018 


1.012 


End of anneal 


1.006 


0 1 2 3 
Time (Hours) 
Fig. 10. — Increase of frequency f(t) with time ¢ at the constant temperature of 
493 °K for silver and of 502 °K for gold. 


The upper diagrams of Fig. 6-9 show the influence of anneal upon the 
frequency of the lowest flexural mode of the specimens. For silver and gold 
(Fig. 6, 7, upper diagrams) this influence is similar to that already observed 
in copper (12). The frequency vs. temperature curve is shifted after anneal 
towards higher values of the frequency. Its inflexion decreases with the height. 
of the dissipation peak, as it is required by the relaxation theory, whilst its 
slope is little or not changed at temperatures higher than Tn. After an high- 
temperature anneal, when the dissipation peak has been cancelled, the fre- 


quency curve no longer exhibits any trace of inflexion. 


(12) See the paper quoted in footnote (?), p. 291. 
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\ _ In palladium and platinum (Fig. 8, 9, upper diagrams) the effect. of anneal 
on the resonant frequency is less obvious, as the frequency vs. temperature 
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Fig. 11. — Influence of cold work on the height of the dissipation peak and on the 
frequency relaxation in silver. Curve A, after an anneal of 4h at 493°K; B, after 
d 2.7% permanent axial strain. 


curve, even after a complete anneal, is a curved line turning its concavity 
towards the high frequencies, in agreement with the previous remark on the 
elastie modulus of the transition elements. However a slight reduction of the 
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curvature is observed after each annealing treatment and this is what one é 


would expect from the relaxation theory. i 
175 — eq 
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Fig. 12. — Influence of cold work on the height of the dissipation peak and on the 
frequency relaxation in gold. Curve A, after an anneal of 7h at 502°K; B, after 
3.8% permanent axial strain. | 
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Fig. 13. — Influence of cold work on the height of the dissipation peak and on the 
frequency relaxation in palladium. Curve A, after an anneal of 6 h at 423 °K; B, after 
| 6.9% permanent axial strain. 


7 


10. 


TONE i P. G. BORDONI, M. NUOVO and L. VERDINI 


The influence of cold work is the opposite of that of anneal. If the peal 
has been cancelled by a high temperature treatment, as it is the case wi 
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Fig. 14. — Influence of cold work on the height of the dissipation peak and on the fre- 
quency relaxation in platinum. Curve A, after an anneal of 1 h at 793 °K; B, after 
3.2% permanent axial strain. 


silver after the treatment (d) of Table I, a permanent strain of 2.7% due to 
an axial compression of the plate at room temperature reintroduces the dis 
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sipation maximum (Fig. 11, lower diagram, curve B). At the same time an 
inflexion appears in the frequency vs. temperature curve (Fig. 11, upper dia- 
gram, curve B). The same influence is observed in gold, when it is permanently 


strained 3.8% after a complete anneal (Fig. 12) and in palladium after a per- 


manent strain of 6.9% (Fig. 13). In the last case the peak is even higher than 
in the original state of the material corresponding to a permanent strain of 
only 1% (Table I). 

As it has already been said, even an annealing treatment at 793 °K does. 
not remove completely the peak in platinum, but reduces it to about one tenth 


of its original value. This small peak is increased again to nearly twice its 
‘original value by a permanent strain of 3.2% (Fig. 14). Hence the influence 
‘of cold work is the same as in the other metals. 


All the previous measurements have been made on polycrystalline speci- 
mens. However the dissipation peak cannot be due to the grain structure, 
owing to its low temperature and to the fact that in copper the same peak 


has been observed by various experimenters in single crystals and in the poly- 


crystalline material (1). 

The opposite influences of anneal and cold work point out that the relax- 
ation effect is due to some geometrical imperfections of the crystal lattice. 
As it has been suggested by SEEGER (!*) the comparatively high value of the 
annealing temperature required to cancel the peak excludes that the relaxation 
effect can be due to the reorientation of divacancies or similar lattice defects 
created by the plastic deformations as these defects anneal out rapidly at room 
temperature, and indicates that the peaks are produced by the motion of dis- 
location segments near their equilibrium positions. 


3. — Discussion. 


3°1. Activation energy. — The straight lines of Fig. 5 show that the vibration 
frequency fm is related to the temperature T,, at which the dissipation peak 
is observed by an Arrhenius equation 


i 


È [= foo 0 |= gr 


where W is obtained from the slope of the lines, and characterizes the activ- 
ation energy associated with the displacement of the peak; fm. is obtained from 


(13) See P. G. BorponI: Journ. Acoust. Soc. Am., 26, 495 (1954); Ric. Scient., 23, 
1193 (1953); H. L. CasweLL: Journ. Appl. Phys., 29, 1210 (1958); V. K. Part: Thesis 


Cornell University, AFOSR-TR 58-92. 


(14) A. SerGER: Phil. Mag., 1, 651 (1956). 
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the intercept of the straight lines with the frequency axis and gives the limiting 
value of the frequency at which the peak would be observed for an infinite 
value of the temperature T,. 

From a purely formal stand point a new parameter, the characteristic 
time 7, may be substituted to f,, by means of the relation 


The temperature dependence of 7 is obtained from (1) and (2) 
7] 
7 


(3) T= di, exp KT, 


where 7, is the value given by (2) when fmo is substituted to f.,,. 

The values of W and 7, obtained from the experimental data by means 
of (2) and (3) are listed in Table III. For comparison purposes the values 
obtained for the same parameters in previous investigations on copper (1) 


TaBLe III. — Fundamental paramelers of the relaxation effect. 


î = = Debye’s = = 
* -1 Tot, W 
Metal Lt) W To frequency vp| o8/u = 
treatment (eV) (gs) (81) btu 
Silver (**) (b) 0.124 40-10" 45-10 6.0:10-4 0.04 
(+0.005) 
Gold (a) 0.158 0.7 30 8.6 Due 
| (+0.002) 
Palladium (a) 0.260 12.0 DINI 10.3 0.3 
(+0.013) 
Platinum (a) 0.192 0.06 47 SI 23.0 
(+0.006) 
Copper (*+*) — 0.122 3.8 67 4.2 0.6 
| (+0.C05) 
| Aluminium (+*s) — 0.11 0.13 82 5.0 24.0 
(*) For the details of thermal and mechanical treatments, see Table I. 
(**) Measurements made on the specimen no. 3. 
(*,*) See the paper quoted in footnote (*), p. 296. 
(x*x) See the paper quoted in footnote (2°), 


(7°) See the paper quoted in footnote (2), p. 297. 
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and aluminium (?*) are also listed at the end of the Table. It is worth of notice 
that the activation energies W for the different metals are rather near each 
other, the largest value being that measured in palladium (0.260 eV), and the 
smallest that found in aluminium (!) (0.11 eV). This similarity shows that W 
is associated with some intrinsic property of the dislocations and not with their 
density or with the impurity content. 


The inverse 7, of the limiting value of the characteristic time is rather 


close to the Debye’s frequency », computed by means of the simple equation 


() y=», 


kOy 
h 


and given in the fourth column of Table III (*). 
The experimental values of W and T,* may be employed to test the validity 
of the present theories of dislocation motion. As it has been shown in a recent 


‘analysis (19) the fundamental hypotheses adopted by WEERTMAN (?°), SEEGER (21) 


and by DONTH (??) are rather different, but the final formulae for the numerical 
eomputation of the potential barrier H which opposes the motion of dislo- 
cations and of the frequency v associated with this motion (?) are surprisingly 
similar. 

In Weertman’s theory the value H,, of the potential barrier is given by 


9 0\ 1 


5 Hw=- pa (Le) 
Re, Tye: | 
(5) Ww he ) 


where b is the length of Burger’s vector; u is the torsional modulus; a, is the 
Peierls’ stress for the limiting case of a vanishing thermal agitation. 
In Seeger’s theory the potential barrier H, depends upon the stress pro- 


(15) N. G. ErnsprucH and R. TRUELL: Phys. Rev., 109, 652 (1957). 
(1?) An even smaller value of W seems to be characteristic of lead. See W. P. MASON: 
Physical Acoustics and the Properties of Solids (New York, 1958), p. 267. However 


| the experimental evidence for this metal is still too limited to allow any conclusions 


to be drawn safely on the value of W. 

(18) A more refined computation of the Debye’s frequencies associated with oe 
different types of waves is obviously not necessary as the comparison between T, * 
and » has only a qualitative meaning. 

(1°) P. G. BorpoNI: Dislocation Relaxation at High Frequencies, in Suppl. Nuovo 
Cimento, vol. XVII, pag. 43 (1960). 

(59) J. WEERTMAN: Journ. Appl. Phys., 26, 202 (1955). 

(21) See the paper quoted in footnote (!*). 

(22) H. Donrn: Zeits. f. Phys., 149, 111 (1957). See also A. SEEGER, H. DonTH 
and P. Prarr: Discussions Faraday Society, 23, 19 (1957). 

(23) In this case the word « frequency » has only a statistical meaning. 
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duced by the vibrations. However, when the strain is from 10-8 to 10° as i 
the present measurements, a satisfactory average of H, is given by the very 
simpie relation 


(6) He= 1.95 Ho. 


In the more recent theory due to SEEGER, DonTH and PFAFF (”), the relation 
between » and the temperature is not exactly of the exponential type. As the 
temperature interval in which the peaks are observed is rather limited, the 
temperature dependence of the dislocation frequency v can be represented with 
a satisfactory approximation by an Arrhenius equation, with an activation 
energy Hspp given by 


(6°) Asp = 14H, . 


The limiting value 9, assumed by the frequency associated with the dis- 
location motion for very high values of 7 has the following expression ac- 
cording to SEEGER 


én (o 
(7) - (%)s= by2a (©) ’ 


where c, is the velocity of transversal waves. In the S.D.P. theory the coeffi- 
cients (¥%)spp of the exponential depends upon the temperature according to 
the equation 


v5 , CkT 


ay) i, (%)spp =-73 16 ar a 


where e is the base of logarithms and % is the Boltzmann’s constant. However 
it must be remembered that the parameter given by the experiment is not 
the frequency fn» but its logarithm. Hence in the narrow interval of tempe- 
ratures in which the peaks are observed the temperature dependence of CAR 
can be neglected, taking for 7 an intermediate value between those listed in 
Table II for 7,, (24). When this is done the approximate values of (Yo)gpp are 
near to those given by (7). 

A first check of the theory may be obtained assuming with SEEGER that 
the value H of the potential barrier equals the activation energy, that is 


(3) Haye 


(#4) The relative error introduced in this way in the computation of In (Yo)app 
does not exceed +6%. For the derivation of (7') see the paper quoted in footnote (19). 
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and employing the experimental data for W given in Table III to compute 
the ratio of the Peierls’ stress 0° to the torsional modulus 4 by means of (5) 
and (6). 

The values obtained are listed in the sixth column of Table III and are 
of the order of a few units in 10-4, that is somewhat larger than those obtained 
in static experiments. This has been explained by SEEGER (?°) observing that 
a great portion of dislocations in a crystal does not lie along close-packed di- 
rections. Therefore their effective Peierls’ stress is several orders of magnitude 
smaller than 0°, and, even at low temperatures, the crystal deforms at stresses 


‘smaller than o°. Only if all dislocations involved in an experiment lie along 


a crystallographic direction it is possible to observe the Peierls’ stress in static 
experiments. 

Another significant test of the theory can be made by eliminating (05/1) 
between (5), (6) and (7) and observing that the ratio (¢,H,)/[27(%)sb*u] is an 
invariant for all metals. If we make the reasonable ua that the 


‘unknown ratio between the theoretical value of (v)s and the experimental 


value of 7,1 has the same value for all metals having the same structure (5 
the above expression must also be invariant when the experimental values 
of W and 7," are substituted to H and 27(vo)s. The last column of Table III 
shows that ‘this is not the case. As the hypothesis (8) cannot be reasonably 
questioned, it seems that the lack of agreement between theory and exper- 
iment is essentially due to the evaluation of (»,), which is the less satisfactory 
part of the present theories of dislocation motion. 


3:2. Relaxation spectra. — It has been shown elsewhere (””) that all the 
experimental points obtained in energy dissipation measurements made at | 
different frequencies must fall on a single line when the ratio Q-1/Q,,* is 


plotted against Wk-(T-1 — T-), provided the relaxation effect is due to the 


sum of many elementary effects having the same activation energy W. 
Fig. 15-18 show that this is the case for the four metals; it is also evident 
that the experimental peaks are much broader than the theoretical curve for 
a single relaxation time (dotted line), and that the low-temperature branch 
is in every case higher than the high-temperature one. These results are quali- 
tatively the same which have been obtained for copper, and they show that 
the dissipation peak is associated with a spectrum of relaxation times. The 


(25) See the paper quoted in footnote (™4). The value of o)/u would essentially be 
the same if Hw or H,,p were taken instead of Hg. 

(28) SEEGER takes (v%)s3 = To 1/(27). 

(27) See the paper quoted in footnote (*), p. 300, or P. G. BorDONI: Theory of rela- 
ration effects with a continuous spectrum, in Proc. of the Third International Congress 
of Acousties (Stuttgart, September 1959), in press. 
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activation energy W associated with each time t has essentially the same 
value W which characterizes the frequency dependence of the temperature 
of the dissipation peak. 


WKN T_T) 


Fig. 15. — Values of the ratio @-1/Q7) as a function of Wk1(T} — T-) for silver, after 

the treatments (a) and (b). Measurements at: (@) 19 kHz; (x) 51 kHz; (©) 421 kHz; 

(A) 820 kHz. Heavy line, sech [0.27 Wk-1(77} — T-1)]; dotted line, Zener’s dissipation 
curve for a single relaxation time. 


SR The symmetrical part of the dissipation curve, plotted as a function of 
Wk-(T71— T-), is due to a symmetrical spectrum of relaxation times, whose 


-10 Sea) 0 +5 +10 

Wk' (Tg 7) 
Fig. 16. — Values of the ratio Q-1/Q71 as a function of Wk (T,, — T+) for gold, after 
the treatment (a). Measurements at: (@) 16 KHz; (x) 56 KHz; (©) 561 kHz. Heavy 


line, sech [0.35 Wk(T}— T-)]; dotted line, Zener’s dissipation curve for a single 
relaxation time. 


center coincides with the characteristic time 7 which is directly related through 
eq. (1) and (2) to the frequency dependence of the dissipation peak. 
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The width of the spectrum can be evaluated with good approximation by 
applying to the high temperature branch of the experimental curves, which 


0 
-10 Sb 0 +5 + 10 
4021 I 
Whol T= 3) 
Fig. 17. — Values of the ratio Q-1/Q;,' as a function of Wk(T;,,. — T-) for palladium, 
after the treatment (a). Measurements at: (®) 25 kHz; (x) 88 kHz; (o) 388 kHz. 


Heavy line, sech [0.22 Wk-(7;,1— T)]; dotted line, Zener’s dissipation curve for a 
single relaxation time. 


is less affected by the asymmetry of the spectrum, the Fuoss-Kirkwood equa- 
tion (8) 


(9) ee = sech [yWk-(T-1— T-)], 


m 


+10 


Wik Ki na) 


Fig. 18. — Values of the ratio Q-/Q;,' as a function of Wk-(T,} — T-) for platinum, 
after the treatment (a). Measurements at: (@) 10 kHz; (x) 42 kHz; (0) 112 KHz. 
Heavy line, sech [0.35Wk-1(7;,'— T-1)]; dotted line, Zener’s dissipation curve for a 
single relaxation time. 


(#8) R. M. Fuoss and J. G. Kirkwoop: Journ. Chem. Phys., 63, 385 (1941). For 
additional information on this approximation of the spectrum see also” the paper 
quoted in footnote (2), p. 305. 
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where y is a parameter which characterizes the spectrum width; its value is 
unity for a single relaxation time, whilst it vanishes for an infinite spectrum 
width. According to the theory developed by Fuoss and Kirkwood, the den- 
sity d(7) of the relaxation spectrum is a bell-shaped curve, having its maximum 


TABLE IV. — Relaxation spectra. 


| si Moti | Last treatment (*) y |: Ing 
Silver (**) (a) 0.28 5.30 
| ’ (b) 0.27 5.60 “4 
| Silver (c) 0.30 4.87 
| (e) 0.18 (**") 4 9.05 
(a) | 0.35 3.91 n 
Gold 
| (e) 0.27 (**) 5.60 
(a) 0.22 7.38 
Palladium aa 
| (0) 0.24 ("+") | 6.55 
(a) 0.35 © E, 
(5) 0.35 ae 
Platinum | (c) 0.35 | 3.91 
(d) 0.35 i | | 3.91 
(e) 0.24 (*»*) | 6.55 
Copper | a 0.39 (+**) 2.89 
> =f 


(*) For the details of thermal and mechanical treatments, see Table I. 

(**) Measurements made on specimens no. 1 and 2. 
(*,*) Computed with the value of W given by Table III and measured before cold work. 
(«*x) See the paper quoted in footnote (*), p. 209. 


at 7; the times 7, and 7; for which the height of the spectrum is 1 /V2 of the . 
maximum height may be taken as upper and lower limits of the spectrum; 
the logarithm of their ratio gives a physically significant measure of the 
spectrum width. 

The heavy lines of Fig. 15-18, which have been computed by means of (9), 
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fit well the experimental points for 7> T,, showing that the Fuoss-Kirk- 
wood theory gives a good approximation of the symmetrical part of the spec- 
trum. The values of y and of In (t,/t,) have been computed for all the meas- 
urements made, and the values obtained are listed in Table IV. 

The widths of the spectra are somewhat larger than the value previously 
found in copper which is given for comparison purposes at the end of the same 
Table. It may be noticed that anneal decreases the width, whilst the opposite 
effect is produced by cold work. 


3°3. Frequency relaxation. — As it has been said, the experimental values 
of (f—f,)/f, are available only for the lowest vibration modes (Fig. 1-4, upper 
diagrams). 

For silver and gold the agreement between theory and experiment is satis- 
factory. According to the Fuoss-Kirkwood’s computation, the maximum 
dissipation @,,! is related to the total relaxation of the resonant frequency by 
the following equation 
(10) VA Mice Sa 

fr 2 

‘where S is the relaration strength, and (f,—f,) is the value taken by the 
‘difference ( (f—f,) for T<T,,, that is when 7 is much longer than the vibra- 
tion period (?°). 

| The ratio (f,—fx)/f, can be obtained from the upper diagrams of Fig. 1 
and 2, extrapolated towards the low temperatures. The values given by this 
extrapolation are in satisfactory agreement with those computed by means 
of equation (10) taking for Q_} and y the values given respectively by Tables II 
and IV. It must however be noticed that the inflexion point of the curve 
which gives (f—f,)/f, as a function of 7, is not found at the temperature 
W.as theory requires, but at a somewhat lower temperature. 

The above comparison between the maximum dissipation and the relax- 
ation strength cannot be made for palladium and platinum, as the compu- 
tation of f, cannot be made by a linear extrapolation owing to the anomalous 
behaviour of the modulus vs. temperature curve which has already been 
pointed out. 


34. Influence of anneal and cold work. - The influence of anneal on the 
relaxation effect, that is the reduction of the dissipation peak and of the 
relaxation observed in the resonant frequency (Fig. 6-9) is readily explained 


(5°) The difference (f, —f,) must be computed at temperatures lower than 7, 
but not too near to the absolute zero, otherwise the extrapolation which gives f, as 
4 linear function of 7 is not valid. 
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by a decrease of the number of dislocations which are active in the relaxation 
process. The results obtained throw also -some light on another effect due to 
anneal, that is the temperature-independent increase of the frequency of the 
specimen, which was first observed in copper (?°). This increase is shown by 
the upper diagrams of Fig. 6-9. In silver, as in copper, a contemporary in- 
crease of the room temperature dissipation is observed (Fig. 6, lower diagram). 
No change in dissipation is found in gold (Fig. 7, lower diagram) whilst in 
palladium and platinum anneal gives rise to a decrease of the room tempera- 
ture dissipation (Fig. 8, 9, lower diagrams). « 

The relative values of the permanent changes observed in frequency and 
dissipation after anneal are listed in Table V, together with those obtained for 


Taste V. — Permanent changes in resonant frequency and energy dissipation at room 
temperature after anneal. 
Changes in 
Metal Speci- Last ewes 
men treatment (*) ci Dissipation 
requency (Af/f) AQ= 
(AQ) 
(c) + 9000-107 + 80-10 
Silver 3 | | 
(d) +11 600 +130 
Gold — (b) + 6800 ae 
Palladium — (b) + 320 — 10 
(b) + 460 — 25 
Platinum sea 
(c) + 630 — 27 
3 I in ee Log + 980 + 29 
Copper (**) 2 6h at 516 °K + 2370 +130 
1 (b) 1h at 873 °K + 2920 +120 
(*) For the details of thermal and mechanical treatments, see Table I. 


(**) See the paper quoted in footnote (*), p. 291. 


copper in the previous work. In every case the changes in frequency are at 
least one order of magnitude larger than those measured for Q-1. Moreover 


(3°) See the paper quoted in footnote (2), p. 291. 
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the former changes are always positive, whilst the dissipation may either in- 
crease or decrease after anneal. Hence no direct relation seems to exist be- 
tween the influence of anneal upon the room temperature frequency and upon 
the dissipation at the same temperature. ; 

The annealing temperatures at which the permanent increase of frequency 
begins to be observed show that this effect in some cases may be related to 
the recrystallisation of the specimen. On the other hand the changes in dis- 
sipation, which have also been observed in copper single crystals, may be 
associated with some changes in the mobility of the dislocations. 

For platinum the values of Q;,1 have been plotted in Fig. 19 in a logarithmic 


0 01 02 03 0.4: 10° 


-1 0,7! 
SR 


Fig. 19. — Exponential relation between the height of the dissipation peak and the 
inverse of the annealing temperature for platinum. 


‘scale against the inverse of the annealing temperature 7,. As the experimental 
points fall on a straight line, the number of residual active dislocations is re- 
lated to the annealing temperature by an Arrhenius equation, with an activ- 
ation energy of about 0.09 eV. This energy is nearly one half of the acti- 
vation energy associated with the relaxation effect (Table III). The reason 
why at room temperature dislocations move near their equilibrium position 
without annealing, is that the time required to form a kink is of the order 
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of %, that is ~10-"s, and the probability of an energy fluctuation higher 
than W is large for such a short time. To make the dislocations inactive with 
respect to the relaxation effect, the energy fluctuations required must last a 
much longer time, and the corresponding probability is small at room tempe- 
rature, notwithstanding the fact that the energy barrier seems to be lower 


than W. 
The influence of cold work on the relaxation effect (Fig. 11-14) is explained 


as easily as that of anneal, by an increase of the number of active dislocations. 
A temperature-independent decrease of vibration frequency after cold work is 
very evident in gold (Fig. 12, upper diagram) and in platinum (Fig. 14, upper 
diagram). This decrease is due to an effective change of the elastic coefficient, 
as the values of frequency given in Fig. 11-14, (upper diagrams, curves B) 
have been corrected for the thickness reduction of the plate and for the increase 
in diameter, due to the permanent strain. 

For a given vibration frequency a slight dependence of the temperature T,, 
upon the annealing temperature 7, has already been observed in previous 
experimental work (?). As 7, decreases together with Q,,! when 7, increases, 
it was concluded that the dislocations more easily removed in copper are 
those associated with the largest values of the relaxation time, as the ave- 
rage value 7 is decreased by anneal. The same dependence is observed in 
the four metals, being very evident in silver (Fig. 6, lower diagram) and in 
platinum (Fig. 9, lower diagram). It may be added that cold-work has ap- 
parently the opposite effect, as 7° for a given frequency increases slightly 
with the amount of cold-work and therefore with the height of the dissipation 
peak (Fig. 11 and 14, lower diagram). 

The most significant parameter to which the observed changes of T,, may 
be related is the total number of dislocations which are active in the relax- 
ation effect. In the Fuoss-Kirkwood approximation this number is propor- 
tional to the ratio Q;1/y, which may be computed from the data of Tables IL 
and IV. The diagrams of Fig. 20 show that the temperature of the peak in- 
creases regularly with the number of dislocations; in the narrow range of values 
considered this dependence may be satisfactorily approximated by a straight 
line. The changes of T corresponding to a given change of the total number 
of active dislocations may be easily computed from the diagrams of Fig. 20 
and from the values of the activation energy given in Table III. For silver 
when Q="/y goes from 1500:10-> to 750:10-* the value of 7 is decreased by 
a factor of about 5. In platinum, when Q;1/y goes from 500-10~° to 250°10-% 
the characteristic time 7 is reduced by a factor of about 2 (*). 


(21) See the papers quoted in footnotes (2), (5) and ("*). 

(32) In the paper quoted in footnote (2), p. 312, a factor of 1.5 was computed fo 
the change of 7 in copper. This value is correct, but one must read: « The change i 
the value of (7,1 — T-1) is of the order of 0.3-107-8 °K ». 
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The above results show that anneal and cold work are effective in 
changing both the number of active dislocations and their average relaxation 
time. In fact when the number is reduced by a factor of 2, the characteristic 
time may be reduced by a larger factor. 
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Fig. 20. — Dependence of the temperature of the peak upon the total number of active 


dislocations. Silver, experimental data reduced to the vibration frequency of 50 kHz; 
| platinum, vibration frequency of 9.6 KHz. 


The experimental evidence which is presently available is not sufficient to 
‘decide if the changes of 7, which are observed at low temperatures, are due 
to corresponding changes of the limiting value 7, of the characteristic time or, 
jas it is more probable, to a slight variation of the activation energy. To give 
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a satisfactory answer to this question the activation energy must be measured 
after different anneals with an accuracy of the order of a few units percent. 


3'5. Subsidiary peaks. — The dissipation curves for palladium (Fig. 3, 8 
and 13, lower diagrams) and those for platinum (Fig. 4, 9 and 14, lower dia- 
grams) clearly show a subsidiary peak below the temperature of the main max- 
imum, of the same type as that found by NIBLETT and WILKS in high purity 
(99.999%) polyerystalline copper, and successively investigated in the same 
material by PARÈ and by others (**). 

In silver and gold the existence of a N.W. peak is less obvious (Fig. 1, 2;. 
6, 7, 11, 12, lower diagrams). This is partly due to the fact that the exper- 
imental values are plotted against temperature. As it has been shown (84) 
the dissipation curves due to a symmetrical relaxation spectrum are not sym- 
metrical with respect to 7 but to 7-1. When the values of the dissipation are 
plotted against the « symmetrical » variable Wk-(T;1 — T-), a flat «bump » is 
observed in the low temperature branch of the curves for silver (Fig. 15) and 
gold (Fig. 16). : 

The N.W. peak is clearly related to the dislocation motion, like the main 
one, as it is reduced or cancelled by anneal (Fig. 6-9, lower diagram) and 
increased or reintroduced by cold work (Fig. 11-14). 

To separate the N.W. peak from the main maximum, the experimental 
values of dissipation have been plotted against 7-1 as it is shown in Fig. 21 
(heavy line). A Fuoss-Kirkwood’s curve (eq. (9)) has been drawn to fit the 
experimental points, and the values of T,,, W, 7 and y given in the Tables Ty 
III and IV have been employed to this purpose (Fig. 21, dotted line A). The 
difference between the experimental points of the low-temperature branch and 
the A curve has been computed, and two maxima have been found for 
T-1= 0.99-10-2 °K-1 and 7-1=1.37-10-2°K-! A second Fuoss-Kirkwood’s 
curve has been drawn choosing the temperature of the first subsidiary max- 
imum 7’ and the new value of the product y'W'%-! by a trial and error method, 
to fit the computed differences in the best possible way (Fig. 21, dotted line B). 
The procedure has been repeated for the second subsidiary maximum, and a 
third curve of the type (9) has been drawn through the new differences (Fig. 21, 
dotted line C). In this way the experimental dissipation curve has been ap- 
proximated by the sum of three « symmetrical » curves, corresponding to sym- 
metrical relaxation spectra. 

The same procedure has been followed whenever possible for the other 
metals, and for the different vibration frequencies. The results obtained, that 
is the temperatures 7’, and 7" of the two subsidiary peaks, their heights 


(33) See the papers quoted in footnote (11) and (18). 
(34) See the papers quoted in footnote (?), p. 304 and in footnote (+7). 
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Fig. 21. — Separation of the N.W. peaks from the main dissipation maximum. Heavy 
line, experimental data; dotted lines A, B, C, computed Fuoss-Kirkwood’s curves. 


TaBLE VI. — Experimental data for the subsidiary peaks. 


Last| Fre- Main peak | First subsidiary peak |Second subsidiary peak 
treat-| quency fe he i 
Metal | nent|'at La RON SA 8 yWk| Tm (Qh) v Wiki Th (ny Wk 
GHZ): (CK) 0" 6 (SK) | CR) | [CEI CK) (°K) 
| ie a È Mei | È | 
| | ines. 4 
Silver (**) | (6) |421.4 | 91.0288-10-5 430 | 61 |120-10-) 460 | 
: (a) 15.99 |117.7\202 63a 090) p90 ; 
Gold — - -—- = —|— 
(a) | 55.60 |127.7/221 (635) 74 Soe 110 | | — 
(a) | 24.92 |180.0| 30 | 650 | 107 20) 1650) 10 )5 
| — — - n 
alladium | (a) | 88.35 |186.5| 59 650 | 111.5 36 650 |<80| > 20 = 
a Ai er da daga Sa a) 
(a) |387.9 |200.0| 75 650 | — |>.32| = 
(a) 9.60 |173.2| 90 780 | 101 49 630 73 32 630 
Platinum - 
(a) | 42.41 |201.0/116 780 | 111 65 630 


(*) For the details of thermal and mechanical treatments, see Table I. 
(**) Measurements made on specimen no. 3. 
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(Qi) and (Q,,)-4 and the values of the product y'W'%-1 for the first subsidiary — 
peak are listed in Table VI. 

The experimental data are too limited for a systematic evaluation of the 
wetivation energy W' associated with the first N.W. peak. However, the energy ¢ | 
W'= 0.16 eV which may be computed for palladium from the values of fre- 
quency and temperature given in Table VI, and the energy W'= 0.11 eV 
which is obtained for platinum, are both consistent with the hypothesis that 
at high temperatures the corresponding values of the characteristic time are of — 
the same order as the period of the Debye’s frequencies. As it was found for _ 
copper, these energies are considerably smaller than the energies of the main 
peak. In silver and gold the N.W. peaks fall below the experimental range 
of temperatures and no information can be derived from Table VI on their 5 
activation energies. : 
4 


Li 
4. — Conclusions. È 

The fundamental parameters W and 7, of the relaxation effect due to dis- i 
locations in silver, gold, palladium and platinum have values near to those è 
found for the same parameters in other f.c.c. metals. This seems to point — 
out that the motion of dislocations has an « intrinsic » character, essentially in- — 
dependent of the impurity content, in agreement with the present theoretical : 
treatments of the same problem . 

From a quantitative standpoint, the agreement betwen experiment and the- 
ory is satisfactory as far as the activation energy is concerned. The evaluation | 
of the characteristic time must be improved, as the present values hardly : 
give the right order of magnitude. In this respect the theory must also be 
extended to account for the wide spectrum of characteristic times which is 
always associated with the effect, and for the relation between the total number 
of dislocations and the value of 7 at low temperatures. 

The asymmetry of the dissipation vs. temperature curves is a general pro- 
perty of the effect, and is possibly due to subsidiary peaks on the low temper- 
ature side of the main dissipation maximum. These peaks are of the same type 
as that found by NIBLETT and WILKS in copper. The dislocations responsible 
for this asymmetry seem to have a lower activation energy than those asso- 
ciated with the main peak. 

The influence of anneal and cold work on the relaxation effect is now satis- 
factorily understood. The new results on silver, gold, palladium and platinum 
have removed an apparent contradiction with theory, showing that the 
temperature-independent increase of frequency which is produced by anneal, 
and the opposite effect due to cold work are not related to the change observed 
in dissipation but are probably due to recrystallization effects. 
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or! J.-P. AMIET (**) 


Instituts de Physique des Universités - Genève et Neuchatel 


(ricevuto il 4 Febbraio 1961) 


Contents. — 1. Introduction. — 2. Elimination partielle des mésons. — 
Bo 8. Introduction d’un potentiel scalaire. — 4. Autre méthode d’élimination 
4 des mésons. — 5. Résolution formelle de l'eq. (3.1); construction de la 


fonction d’état générale. — 6. Calcul de la source du champ Vy. — 7. Caleul 
d’un exemple. — 8. Niveaux d’énergie. — 9. Parité. — 10. Utilisation des 


fonctions d’onde; moments magnétiques. — APPENDICE A) Unités et 
symboles. — B) Champs d’opérateurs mésoniques Vx). — 0) Parité M. 


D) Quelques proprietés utilisées des éléments de matrice de V,(x). 


4. — Introduction. 


i Pour autant que la théorie des champs soit aussi fructueuse dans le cas 
des nucléons que dans celui des électrons, les solutions stationnaires de Véquation 


(1.1) A E I d*x(VV,-VV=4 I-II, + 2V,V%) + 


i J f 
al 3 * / ‘ x 
fre | 3 road VV, Li 


2m 


apporteraient la contribution essentielle a la description des noyaux. (1.1) est 
la première approximation non-relativiste communément admise de l’équation 


Ss (*) Recherche subventionnée par le Fond National Suisse pour la Recherche Scien- 
tifique. 
(**) Ce travail constitue une thèse présentée à l’ Université de Genève. le 21 Mars 196), 


pour l’obtention du grade de Docteur ès Sciences. 
y (***) Actuellement a l'Institut fiir theoretische Kernphysik, Universitat Bonn, 


Deutschland. 
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fournie par la théorie du couplage P.S.-P.V. et devrait suffire pour les basses 
énergies. Malheuretsement, ses solutions ne sont pas encore connues, et l’on 
ignore méme s’il en existe. Pourtant les mésons x jouent un réle prépondérant 
dans la liaison des nucléons au sein des noyaux. Les expériences de diffusion — 
pions-nucléons révèlent une forte interaction et parmi toutes les particules ayant 
une forte interaction avec les nucléons, les mésons z sont celles dont la masse 
est de loin la plus petite, ce qui entraine qu’ils engendrent les forces de rayon 
d’action le plus grand. D’autre part, les bons résultats obtenus en introduisant 
un potentiel moyen dans les modèles phénoménologiques 4 particules indé- 
pendantes (1) deviennent incompréhensibles si ’on considère des forces de rayon 
d’action beaucoup plus petit que celui des forces 7-mésoniques. Par là, seuls 
les mésons 7 restent en cause. 

Un essai de dériver de l’interaction méson z-nucléon un modeéle simple 
décrivant le comportement des noyaux sphériques moyens et lourds aussi bien 
que le modèle phénoménologique en couches, est par conséquent justifiable. 
En termes plus précis, il s’agit de trouver l’origine mésonique des deux termes 
essentiels du modèle en couches: le potentiel moyen et le couplage « spin- 
orbite »; c’est le sujet que nous abordons ici. 

La méthode des forces statiques à deux corps, utilisée fructueusement pour — 
le deutéron (?), conduit pour les noyaux lourds 4 un probleme d’une extréme ~ 
complexité (?), et nécessite l’introduction d’un cut-off; aussi l’éviterons-nous 
en utilisant une méthode simple qui conduit rapidement è un modele en 
couches, contenant un terme jouant le méme réle que le couplage L-S. Le 
potentiel scalaire ne vient pas directement, nous en verrons la raison à la ~ 
Section 3, et doit étre introduit ad hoc. 

Si Pon considére un système de A nucléons dans un état stationnaire, on 


peut transposer (1.1) dans l’espace de configuration et l’écrire directement: — 


HY = (H+ HS + HY = BY, i 
Hy = d?a(VV,(x)-VV%(x) + TI*(x)IL,(x) + yw? V,(x) V=(x)), 3 


A if A 
Hy = 3-24, Prat Se ev Vie 
a=1 


{ 
q 
Ma=1 j 
“ 
5 


x 
PSP. SMS ORI ee 


(') K. BLEULER et CH. TERREAUX: Helv. Phys. Acta, 28, fasc. 2-3 (1955); 30, 
fase. 2-3 (1957). M. GoreppeRT MayER et J. H. D. Jensen: Elementary Theory of ki 
Nuclear Shell Structure. È 

(?) L. HuLtHÉN et M. Sucawara: The two-nucleon problem, in Encyclopedia of — 
Physics, vol. 39. 

(*) A. M. FeIxGOLD: Phys. Rev., 101, 258:(1956); T. Terasawa: Progr. Theor. Phys. 

22, 150 (1960) et 28, 87 (1960). 
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MODÈLE SEMI-PHÉNOMÉNOLOGIQUE DU NOYAU ETC. |, =O 
Xay Say ta = variables d’espace, de spin et d’isospin du a“ 
pres d’occupation des mésons (*). ; 

On peut éliminer partiellement de (1.2) les variables mésoniques, sans passer 
| par les forces à deux corps. Il suffit de postuler un Y de forme particulière, 
puis de prendre la moyenne de H relativement aux mésons. On obtient par 
ce procédé un système d’équations différentielles régissant les fonctions qui 
décrivent le mouvement des nucléons. 


nucléon; M*,, = nom- 
m 


2. - Elimination partielle des mésons. 


Considérons une fonction d’état dans laquelle les variables nucléoniques et 
| mésoniques sont séparées: Y = g(X,n)|M}); X= (1, ...; x,) symbolise les vari- 
ables d’espace des A nucléons et 7 leurs variables de spin et d’iso-spin; |M} 

est une combinaison linéaire de vecteurs 4 nombre fixe de mésons ("*); les 
normes des trois fonctions valent (2, ¥)=(y, y)=<M|M>=1. Introdui- 
sons 7 dans Véq. (1.2) et multiplions scalairement chaque membre par <M |; 
on obtient l’équation différentielle: 


A 1 eee N ) CORE ar 
(2.1) D- - Im dat on Vi J ete p(X, n) = Exg(X, N) 
ou 

V.(x,) = <(M|V.(x)| ©), Ey=E- M|\H,|M}. 


a 


Les fonctions V,(x) = <M|V,|M> sont inconnues car on ne connait pas 
|M), mais on peut établir une équation différentielle qu’elles doivent satisfaire. 
L’éq. (1.1) est une première approximation non-relativiste d’une équation hamil- 

‘tonienne équivalente au système d’équations de mouvement: 


(2.2) | (79,4 m + I 7° y%VA) Pa) = 0, 


(2.3) (O~ Wala) = Late). 
En prenant la méme approximation non-relativiste de (2.3) et transposant 
la source dans l’espace de configuration des nucléons: 


A A 3 È 5 la) Ha 
(A a pu) V(x) Sa I Lo Xo) Ta 2 Vo (È <a cd = 0x (*#") . 


(*) Dans l’appendice sont rassemblées les notations particuliéres, des définitions, 
et quelques propriétés des opérateurs utilisés. Etant donné sa briéveté et son carac- 
tère de memento, nous l’indiquons en bloc dans les renvois. 

(**) Appendice. 

(*»*) 6/dai=dérivée è gauche et 0/0a*= dérivée a droite. 
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La valeur moyenne de cette équation pour la fonction Y donne l’équation 


cherchée: 
(2, (A—w)V (x) ¥) = (gp, pI(A — 29M V(x) | MD , 
(Y, 0.7) = <M| MY 2 g*(X, 0) p(X, n), 
d’où i 
(24) (4 — 22)Va(x)= -_ £v i X, 1) t102P(X; M)Ò(x — Xo) = 0a(x) - 


Les éq. (2.1) et (2.4) constituent un système d’équations différentielles 
couplées dont les solutions présentent une particularité: prenons en une, 
{g(X, n), V,(x)}, correspondant a l’énergie E,. La fonction P(X, n)=g(-4, n) 
satisfait naturellement . 


4 1 ae by 
(2.5) DI (- 4 7 VARI x,)x4) P(X, 7) = Hy P(X, n) , 


(o Pel 2 


pour la méme valeur £,; la source engendrée par @ est: 


si g donne g,(x); 6, est source d’un nouveau potentiel V(x) =—V,(— x). En 
introduisant UA dans (2.5) on retrouve la forme (2.1). Ainsi, le couple 
{P(X, n), 7. (x)} est aussi solution du système (2.1-4) pour la valeur H, de 
l’énergie; c’est une dégénérescence, d’un type particulier. 

En supposant que Y= g|M) représente un état d’un noyau, ou inversé- 
ment que les états des noyaux sont représentables par des fonctions où les 
variables mésoniques et nucléoniques son séparées, le système d’équations (2.1-4) 
constituerait un modéle permettant de décrire le mouvement des nucléons. 
Les équations de ce modéle seraient exactement, d’autre part, celles que fournit 
la théorie quantique classique du couplage pseudo-vectoriel, qui consiste è 
envisager classiquement les équations de mouvement (2.2-3) des opérateurs. 
Selon ce point de vue, après une approximation non-relativiste et statique, 
(2.2-3) ménent 4: 


(2.6) i Sah A vf o VV, (e) Ya) = Ba zal) ; 
(2.7) (i be 2)Va (x) = — “ los V-(YiG0t 2 %a(X)) FI 


où a numérote les nucléons du système. 
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Les variables x, étant séparables directement dans (2.1), on trouve (2.6) 
en posant g(X, 7) = gi(x1) --- X,(x,); du méme coup, 0, de (2.4) prend la forme 


de celui de (2.7) si on norme les ya. 


Un tel modéle, bien sùr, n’a pas-de sens physique car les fonctions d’onde 
des nucléons n’ont pas une parité définie, et par conséquent g et Y non plus; 


ce fait provient de ce que [z, o: V V,(x)] #0 (x étant l’opérateur de parité 


défini par 7 y(x) = x(- x)) (*), et il est lié au caractère pseudo-scalaire du champ 
des mésons x. On serait tenté d’utiliser la dégénérescence vue avant, et de la 
lever en introduisant un terme supplémentaire après diagonalisation des solu- 


tions selon 7. En fait c'est impossible car on a affaire a une dégénérescence 


des solutions du système (2.1-4), et non pas de (2.1) seule, et on ne peut com- 


‘piner linéairemente deux solutions {p, V,} et (9, 7} car ¢ intervient non linéai- 


rement dans la source 0,. A titre illustratif on pourrait dire que les deux 
solutions représentent deux « noyaux » fictifs, « miroirs » l’un de l’autre, et seul 
un mélange des deux peut correspondre à un noyau réel. 

Ce mélange, impossible classiquement, peut étre essayé en conservant les 
mésons: dans la fonction .d’état. Remarquons tout d’abord que la seconde 


| solution {7, V,} peut étre obtenue par le méme procédé que la première, en 


introduisant dans (1.2) une fonction =F | M> puis en prenant la moyenne 
relativement è <M |; si | I est tel que <M|V,(x) | =<M|-V,(-x)|M>, alors 


| 7 satisfait (2.5) avec le méme spectre de valeurs propres que celui de (2.1). 
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Le vecteur mésonique | > peut étre facilement défini si l’on fait une sup- 
position, nécessaire par ailleurs. Nous postulons que les vecteurs |M} ne con- 
tiennent que des mésons S (moment angulaire {= 0 (*)). Les équations devant 
s’appliquer aux noyaux, il faut qu’elles aient la méme symétrie spatiale qu’eux, 
cest-A-dire la symétrie sphérique (nous ne considérons pas, ici, les noyaux 
déformés). Cela impose que les potentiels pseudo-scalaires ne dépendent que 


‘ de la coordonnée r=|x|: <M|V,|M>=V,(r), done que | M> ne contienne que 


des mésons S (*). Cette hypothese est compatible avec l’éq. (2.4) lorsque les 


| niveaux J sont complets (2J+1 particules), car dans ce cas 0,= 0,(r) (*); 


on peut, sans perdre en généralité, supposer que tous les niveaux sont pleins, 
car les écarts sont faibles dans un modele en couches. Décomposons main- 
tenant | M)> en parties paire et impaire: 


SA A 
HO = Fall +19), 
TI, |=|+>, H,| 2 =—|-, Il, = opérateur de parité dans Vespace des 


mésons; |-+) est une somme de vecteurs contenant un nombre pair de mé- 


(*) Appendice. 
(**) Voir Sect. 6. 
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sons S, et |—> un nombre impair. Substituons cette expressions dans la défi- 
nition de V,: 


V(r) = CMV (8) | AD = 3, |A 


Il est alors clair que: 


10) = Th, |) = (15). 


Ainsi: 7=II9, si II est Vopérateur de parité dans l'espace de Hilbert 
nucléons-mésons (*). Les fonctions Y% et Y conduisent & des « modéles » dif- 
férents ayant méme spectre d’énergie, mais sans sens physique, pas plus que 
les fonctions elles-mémes qui n’ont pas de parité. Malgré cette difficulté, il 
n’en reste pas moins que la séparation des variables n’est pas dénuée de sens 
puisqu’elle correspond à l’hypothèse que les noyaux sont approximativement 
représentables par un modéle en couches, ce qui est vérifié (4). Aussi tentons 
nous l’essai suivant: Prenons Y et Y comme fonctions de base pour un calcul 
de variation, et posons D(x, 8) = a ?+Y (a*x+f*B=1); D(a, B) correspond 
à un mélange des « modèles miroirs ». Puis, calculons les valeurs de « et f 
qui rendent minimum la valeur moyenne <@|H|®) de Vhamiltonien H: 


E(a, 8) = ‘D|H|D) = a*xP|H|P) + *RKP|H| Py) + 
+ ap*<P|H|P) + B*BP|H|P). 


On peut décomposer ¥ et 7 en toute généralité: 


= 9(£, 9): 1) = (p(X, 0) + E Fe = (+= +, 


Si eel 
P= WX, 10 = WF = (VO), 


ou 


gE(X, n) = 3(9(X, mM) +9(-X,m)) YF =yt|t>+oe|F>- 
Posons E* = <P*#|H|Y#); comme <P*|H|P)=0 (*), on a 
E(a, B) = }(a*x + B*B)(E* + E°) + 4(a*B + af*)(E*— E"). 


(*) Appendice. 
(4) K. BLEULER et Cn. TERREAUX: Helv. Phys. Acta, 30, 183 (1957). 
(*) Car LH, II}=0. 
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Soient E,=}E*+F), AE=E-—E*, a=|alexp[iw], B=|Blexp[ix]. 
La propriété «*«+/*6 =1 entraine finalement: 


E(a, B) = E, —|\a|V-1—|a|? cos (0 — x) AE. 


E(«, B) est extrémale pour |«|= 1/V2 et x=, 0+7. Comme © est définie 
à une phase près, on peut poser © = 0. D’où: 


1 
| X= lo) => J? si AE > 0 ’ 

\ E(x, f) est minimum pour 1 
| xa= to) = WD) Sl AH< 0 ; 


et 


1 1 ib ~ < 
o) Iieugg dea, 


cay dee a MM re A 2 rs 
y =2( aa Dean ae di 


Ainsi, les fonction D(«, 8) rendant minimum l’énergie ont une parité déter- 

minée; elles ont un sens physique et pourront étre utilisées pour représenter 
des états de noyaux. La parité de l’état fondamental est déterminée par le 
signe de AH: 


Etat fondamental Y* si AH>0, 
Etat fondamental YW si AH<0. 


Cette question est traitée a la Section 9. 
En appliquant [H°,, I)=[H3,11}=[H,I]=0, ce qui entraine 


(yt, Hye) = Hi (+, <+|H’|—¢*) = 0, 


et en utilisant (p+, )+(g-,@)=<+|t>=1, on obtient 


(2.8) AE = Ko", p*) — (97; pH +2, Be) — (9 Hm}, 
ou 
A 
Pe eh Fh oo Ch | 4) SE i ¥ o,-V. Vite. 
a=1 


Le potentiel V,(r)= <M|V,| M> est évidemment réel, mais pas forcément, 
<+|V,|—>, qui n’est pas un élément de matrice diagonal; il reste VI = partie 
imaginaire de <+|V,\|—> dans AK. 
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Il est possible, par un raisonnement indépendant du précédent, de prévoir 
la forme des fonctions Y les plus simples qu’il convient d’utiliser. En intro- 
duisant ces fonctions dans (1.2) et en prenant des moyennes relativement aux 
mésons, on tombe sur des systèmes d’équations différentielles ressemblant è 
(2.1) mais dans lesquels on ne peut séparer les variables x,. Il est alors pratique 
d’utiliser les solutions de (2.1) qui sont une approximation des solutions de 
ces systèmes. Ce sujet est abordé a la Section 4. 


3. — Introduction d’un potentiel scalaire. 


Un premier test révèle qu’un modéle de noyau basé sur les éq. (2.1-4) ne 
donne pas des résultats satisfaisants. Pour avoir une idée du spectre d’énergie 
nous avons résolu le problème 4 une particule (2.6) en postulant différents 
V.(r). Il s’avére qu’en. aucun cas on obtient un nombre correct de niveaux 
et une bonne succession; il manque visiblement dans l’hamiltonien un terme 
scalaire U(r). 

Un tel terme apparait pourtant dans les théories mésoniques. Si l’on aborde 
le problème du couplage entre nucléons en résolvant (1.2) au moyen de l’ap- 
proximation des forces statiques 4 deux corps et d’un calcul de perturbation, 
on obtient un hamiltonien pour les nucléons seuls composé des termes d’ordre 
0, 2, 4, ete., du calcul de perturbation (les ordres impairs sont nuls) (?). Les 
deux premiers termes peuvent également étre obtenus par la théorie des champs 
classiques (non-quantifiés), mais pas ceux d’ordre supérieur. Le 4” ordre 
contient un terme scalaire d’interaction 4 deux corps, et quelques autres dont 
Vapport a l’énergie totale est pratiquement nul si on prend la moyenne pour 
toutes les particules. Dans son travail sur le deutéron, L. HULTHEN (?) montre 
qu’il faut tenir compte de ce terme scalaire pour obtenir des résultats satis- 
faisants. Il n’est pas étonnant que nous ayons besoin de ce, potentiel U(r); 
on peut montrer que l’équation hamiltonienne pour A nucléons ne contenant 
que les termes d’ordre 0 et 2 correspond a l’éq. (2.6) pour un nucléon plus 
une équation pour V,. En outre, par la méme opération, on peut évaluer 
que le terme scalaire du 4"° ordre (avec cut-off) donne un potentiel scalaire 
moyen d’ordre de grandeur semblable à celui des potentiels nucléaires me- 
surés (9). 

Il se pourrait, soit dit à titre de remarque, qu’un perfectionnement de 
notre procédé donne le terme désiré, par exemple en introduisant dans (1.2) 
une fonction d’essai ‘7 de plus de deux termes. Nous ne nous hasarderons pas 
dans une telle recherche; étant donné, d’une part, les considérations ci-dessus, 
d’autre part la nécessité où nous nous trouvons, nous introduisons a priori 


(°) P. HUGUENIN: communication privée. 
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un potentiel U(r) dans l’éq. (2.6); avec cette adjonction phénoménologique 
notre modèle pseudo-scalaire est, par définition, basé sur le système d’équa- 
tions 


(3.1) (1, + È o-VV,1°) Xa(x) = EaXa(x) , 
(82) (A= JT, = — 4 $ Wat eer (x), 


ou H,=—(1/2m)A+ U(r), U(r) a choisir en fonction des résultats expéri- 
I mentaux. 

Les fonctions d’état d’un noyau sont construites 4 partir des solutions 
‘de (3.1): 
PE = g*(X,7)|+>+ pn); 
ou 


p(X, n) = (p(X, 7) + o(—X, n), — P(X, N) = (a) Yulx,) - 
> Les fonctions gt(X, n) et p_(X, n) servent a décrire le mouvement des nu- 
cléons, mais elles ne sont pas des fonctions d’onde au sens courant du terme. 
i Nous voyons à la Section 10 comment les utiliser. Le fait qu’on ne puisse 
représenter un état de noyau par une fonction dont les mésons sont absents 
provient du caractère pseudo-scalaire du champ mésonique. 


Remarque 1. Le potentiel U(r) ne change en rien la structure du modele, 
dérivée par le procédé de la Section 2, qui a l’avantage de permettre l’utili- 
sation du terme o-V V,t* dans un modéèle a particules indépendantes. 


Remarque 2. U(r) englobe des effets qui ne semblent guère explicables au 
moyen des théories mésoniques actuelles, comme par exemple le «coeur ré- 
pulsif ». 


Remarque 3. On peut remplacer (3.1) dans la définition par 


6.1) (Ho-+ 4S o4-V.Va(.)22) AL M= BX, 1) 
ou 


A il 
H, sà. 2 a A,+ U(r,)). 


4. — Autre méthode d’élimination des mésons. 


Au lieu d’essayer au départ une solution approximative de (1.2) revétant 
la forme la plus simple, soit le produit g(X,n)|M}), on peut partir du point de 
vue opposé, c’est-a-dire de la solution exacte formellement développée en 


7 - Supplemento al Nuovo Cimento. 


an ge» 
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série: 


Y= > prlX, n). 


Les fonctions p(X, 7) dépendent des variables nucléoniques, et les vecteurs 
n) des variables mésoniques. Nous choisissons, pour les mésons, la représen- 
tation des nombres d’occupation en coordinnées sphériques. Les vecteurs |n> 
sont des combinaisons linéaires de vecteurs |... M%,,>, Mim = nombre d'occu-. 


pation (*). 

La solution exacte Y est donnée par un développement infini; l’idée d 
Vapproximation consiste à ne retenir qu’un nombre fini de termes, mais sans) 
choisir d’avance les fonctions de base g,, |n>; & chaque nombre de termes” 
retenus correspond un modéle, dans le cadre duquel on calcule ces fonction 
Nous choisirons un nombre de termes tel que le modéle obtenu soit le plus. 
près possible d’un modèle en couches, l’existence de couches dans la structure. 
des noyaux étant un fait acquit (). : 

Le point de départ utilisé à la Section 2, soit la séparation des variables, | 
équivaut è dire, dans la présente formulation, qu’on ne retient qu’un terme. 
de la série. On est alors conduit à un modéèle en couches pur, mais qui ne 
permet pas de concilier deux exigences physiquement nécessaires: parité déter- 
minée des états et sphéricité. Le modèle le plus voisin permettant de satis- 
faire ces deux conditions est obtenu en gardant deux termes: | 


P= p(X, n): |D + p(X, m)|2). 


Les fonctions gn, |”>,-,, ne sont pas absolument arbitraires: la relation 
TY =1I,9,:I1,|1) +11,9,:1I,|2>= + 7 west satisfaite que si Ip, = pny 
II,,|7> =-+|n>; de plus, pour que les équations aient la symétrie sphérique, 
il faut que les éléments de matrice <n|V,(x)|n’> soient des fonctions de r= o 
exclusivement; nous exigeons, en conséquence (*), que |n)= combinaison li- 
néaire de vecteurs |... M%,...> où My 0 si 740 (on ne tient compte que 
des mésons SS). 
En remplacant les indices 1 et 2 par (+) et (—), qui indiquent la parité 

de la fonction, les deux vecteurs d’état approximatifs s’écrivent: 


a 
klm * kim T 


(4.1) Ut ze pt (Xa) | ee eC, I> 
(4.2) Y- = ot (X, n) | + @ (4, )|+>.- 


Introduisons (4.1) dans Véq. (1.2), et multiplions scalairement chaqu 


(*) Appendice. 
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membre par <-++|, puis par <—|; il vient: 


pi A.) 94,0) + (L$ 04-0, ERO 2) gtx, m) = 


eles 
EEN | 
= jie cia XG y 
| 4 |-+> Pp ( ’ n) , 
SA È PE ig Kala] no ; 

La <—|Hi|- : 
a =(z dit 2\ (Xm) 
car <—1V,|+>=0 et <+[H},|>=0(*): 

On peut éliminer les valeurs <+|+) et <—|—) des équations au moyen 
de la transformation: 
PO, TAI NES ~ 1 
t= V+|+)gt, ee 


qui laisse la norme de Y invariante. Cette liberté est naturelle, car les normes 
individuelles de g* et | +) n’interviennent jamais isolément, et seules les normes 
des produits y|> ont un sens physique de probabilité. 

En effectuant cette transformation et posant 


Aa,  Ex=<&|Hx|tb, TM=+]kl®]=>, 


| Hog, + n) +( yo VV) e (Km) = (B— PA, 
(4.3) ; | 
| | Hose n+(Z DIS VAL: ra) ts) px, n) = (B— By) (Xn). 


En introduisant (4.2) dans (1.2), on obtient un système d’équations ana- 
logue pour les états de parité (—): 


A 
Hog*(X, n) + (2 Do, VV) XX, mM) = (B— Ew p(X,) 
a=1 


(4.4) Pa 
[term + (£ È ov) 1) = (8 Bile 9) 


(*) Appendice. 
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Il est essentiel que les systémes soient différents (on passe d’un à autre © 
en intervertissant H* et HZ, ainsi que V, et V*). S'ils étaient identiques, il — 
y aurait dégénérescence par rapport à la parité, car on pourrait construire 
une solution Y*+ et une autre ‘77 a partir des mémes fonctions yt et g-, et des 
mémes kets |+) et |—>. : 

Les systémes d’éq. (4.3) et (4.4) sont trés semblables à l’éq. (2.1); il y a 
tout lieu de croire que leurs spectres d’énergie ne seront guere différents, et : 
que le méme manque de potentiel scalaire se fera sentir. En conséquence nous 
introduisons ici, également, un potentiel phénoménologique U(r); il n’y a pas 
d’équivoque quant a sa place, étant donné ses caractères réels et scalaires; — 
il suffit de l’ajouter au laplacien, et redéfinir î 


(4.5) ato G A Ut) sia 


a=1 


Les systèmes (4.3) et (4.4) présentent l’inconvénient qu’on ne peut les ra- | 
mener 4 une équation pour une particule en séparant les variables; si on pose 


+ 
V(X, n) = Di 
YP 


) fonction & deux composantes, (4.3), par exemple, s’écrit 4 


H, Hh | == Wil. Hi 
| Jola, n+ ‘ eta n)= Eyp(X, n) ’ 


He lhe tH, = Bie 
ou 
> 
; f A R a 
Hy =— Da VaValta)Ta; h= By— Ey, di 
I Ma=1 ‘ 
Vi= Re Var), Vi=ImV,(r), EE HE+E. | 
i = hj)? 24H; e 
L’opérateur Bea ne peut se mettre sous forme d’une somme — 
— iH, h/2 


d’opérateurs partiels relatifs 4 un nucléon, d’une part, et d’autre part empéche | 
de sommer les deux équations d’un système. 

Il faut faire une approximation pour ramener le tout à un modèle a par- 
ticules indépendantes. Le terme génant est probablement petit devant l’autre, 
qui contient le U(r) assurant l’essentiel de la liaison et la partie réelle de V,. 
Négligeons ce terme pour le calcul des fonctions g, et tenons-en compte onesti 
par un calcul de perturbation du premier ordre. Les systèmes (4.3) et (4.4) 
deviennent identiques et on obtient, en sommant les équations 


(4.6) (H, + H,) p(X, n) = Erp, n); 


où Po Pi + Pr 
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On retrouve en (4.6) l’éq. (3.1'); done la fonction q utilisée a la Section 2 
est une première approximation des solutions de (4.3) et (4.4), identique à q. 
En ce qui concerne Vénergie, le calcul de perturbation donne: 


h sr 2 isa =, Core 
AE+ = 5 (Pos Yo) — (wi Go )} + (wi, Hive) — (go, Hige)} ; 
NESTA, 


AE = HE, — Et = AE — AE* = écart entre niveaux 4 parité opposte. Ce 
AE est indentique è celui qui découle du calcul de variation (voir (2.8)). 


5. — Résolution formelle de l’eq. (3.1); construction de la fonction d’état générale. 
Nous résolvons l’éq. (3.1) au moyen d’un calcul de perturbation: 
(3.1) H x(x)= (Ay + H') x(x) = Ey(x), 


H,=—(1/2m)A+ U(r) est Vhamiltonien sans perturbation, et H' = (f/u)o- 
‘VV, (r)t* est le terme de perturbation. 

Nous supposons en cela <H')><<H,>, ce qui se confirme dans la suite. 
Il est commode d’établir d’abord la forme des solutions. L’hamiltonien H 
commute avec les opérateurs de rotation de l’espace autour de l’origine, car 4, 
U(r) et o:VV,(r) ont la symétrie sphérique; par conséquent [H, J?]=[H, J.]=0, 
ce qui permet aux solutions y d’étre fonctions propres des opérateurs de mo- 
ment cinétique total J? et J,. Par contre, H ne commute pas avec l’opéra- 
teur 7 de réflexion de l’expace (parité): [H,, 7]=0, car [U(r), 7] =[A, x] = 9, 
mais [H’, x] #0. Posons 


V i 
7 veto e e 


x 
r dr 


o= 6° 
La propriété: 

(5.1) on + mo = 0 

entraine que oV.(r), done H’, ne commutent pas avec x. Les solutions y 

n’ont par conséquent pas de parité déterminée. L’ensemble des fonctions des 


deux variables 0 et y (*) qui sont fonctions propres de J? et J,, pour un couple 
de valeurs propres déterminées {J, M}, constituent un sous-espace vectoriel a 


(*) x= (r, 6, p). 


1087 


102 J.-P. AMIET 


deux dimensions; puisque [J?, 7]=[J.,a]= 9, ce sous-espace est sous-tendug 
par les deux fonctions propres de x: Di,,(0, g) et ®,,,(0, 9); on peut ainsi écrire i 
directement les solutions sous forme de combinaisons linéaires de ces fonctions | 
au moyen de fonctions radiales: : 


(5.2) Yon(X) = = 20 ale D;,,(9, 9) pato ~ fit ) OF, (9, P) . 


Les fonctions angulaires + on deux composantes de spin, et les fonetions 
radiales f= en ont deux d’iso-spin. 

On peut remplacer les indices + et — par le nombre quantique J, car di, 
sont fonctions propres de L?; nous notons dorénavant les fonctions angulaino i 
D,in-5-1(0, 9) (€E= +1), ou encore D,,,,_131 (€ = +1), M étant-valeur propre | 
de J.. Les fonctions (5.2) deviennent: 


ah 1 1 

(5.3) Ywou(X) = 7 Iva") Dyuz(0, p) + x frare(?)* Dymite(0, P) , | 
| 
où l=J—te, e= + 1, et N est introduit pour compléter la numérotation des | 
solutions. 1 LI 
Les fonction ®,,, ont les propriétés: : 
a) AP su 1) Dora ’ b) so Dns ’ 
ui 
(5.4) e)- PO,,,=IT+ ICE d) CID, = MODE bi 
e) L°D, = (W+1)D,1= (J —te)(J —32e4+1)@,,,, 
f) (O Dry) = dg O55" Bu Ow % 
J 

Calcul de perturbation. — Par un calcul de perturbation, nous allons déve- 
lopper les solutions de (3.1) dans le système des solutions de l’équation A 
5 


(5.5) HyC(x) = E°C(x). 


A cause de (5.3), nous choisissons directement des ¢ qui sont fonctions 
propres de J? et J,: 


ar 1 
(5.6) ; Cnims(X) = n Ini?) i Piys(0, g) ’ J=1 +3 


I,;(7) devant satisfaire 


“i 1a 
(5.7) (Tsi D0)) dul) = Bi gual?) 


} Ss 


2m dr? 2m Tai 
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Du fait qu’on néglige l’interaction électromagnétique, protons et neutrons 
sont traîtés symétriquement, et sont plongés dans le potentiel commun U(r); 
ils ont done mémes fonctions d’onde avant l’introduction de la perturbation 
b) 


Ce n’est plus le cas après, aussi introduisons-nous l’indice ¢ d’iso-spin; par 
définition : 


1 fonction radiale à deux composantes 
Gnu=1(") Ta Ini(?) , x È 
0 attachée è un proton. 


7 (r) La fonction radiale è deux composantes 
u=-1\7) = iP Pe 
‘i ach attachée & un neutron 


Let g,, sont fonctions propres der 7 — i n") T*Qnu= (nun Lies rela- 
tions [H', #] 40 et [H', L?] #0 font que les y ne peuvent étre fonctions 
propres de 7° et L?; il est toutefois pratique de conserver la numérotation 
(nLMJt) pour ces iii 

L’iso-spin et la charge ne sont pas conservés, mais tout comme la parité, 
‘ils sont «répartis » entre les nucléons et. le champ mésonique par le terme 
EVI. 


Notons pour fixer la notation les résultats nécessaires du calcul de pertur- 
bation; en abrégeant (nlmJt) par v: 


1° ordre: 
NES = <y|-H?|y> ; 
(5.8) ial <y' |H'|w (#7), 
Dylyi. TT EX > pi ’ 
2°" ordre: 
7 H' STE r H' : 
po Amey SEMO Es ee BD. 
v'#v v' 4y v'#V 


I Les éléments de matrice de H' sont aisément calculables 


(5.10) 


v|H' |) © Gm | Vama" gni <Pins|o| Pra? = 


[co] 
| 


= 657 Ou'+e'O mm’ thar Gru’) V, alt) U Iner(T) - 
0 


en faisant usage de (5.4 b) puis de (5.4 f). 
Les seuls coefficients non nuls sont done: 


» (5.11) - pi (1, J) = — (E, E MMI H'|n'1+eMJt}. 
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Ils ont la forme explicite: 


1 
sl EI da E, farsa Vara 
rae n'i+e 
0 
(5.12) Brn'(l, J) = Ban *(1, J) = 
SARÒ ol. / . r 
une ee AP Jnil?) Inte?) (Var) — 1 V(r) . 


0 


Les gni(r) sont supposées réelles car nous utilisons dans la suite des fonctions 
radiales réelles. 

La correction du premier ordre AE” est indentiquement nulle 4 cause de 
‘v|H'|v = 0; nous gardons par conséquent le second ordre AE‘? pour l’énergie. 

En résumé: 


1 
(5.13) Ynust(X) = n {Cnn (x) + 23 bee (2, T)Cn'ivemar(X)} (J=14 fe). 


La constante de normalisation N 


inez (1 DARE 
n't! 


| AE us = DB (J, J) ) |? (# wo 
(5.14) 
| Boo ia Ei af AEs 

Remarque: En examinant (5.13) on voit que pour n et J donnés il existe 
deux solutions, correspondant 4 1=J-+4 (abstraction faite de M et t); ces 
deux solutions contiennent chacune les deux mémes fonctions angulaires 
D,.3 y,,(9, y). Ces paires de solutions jouent un ròle important dans les ap- 
proximations futures, et correspondent à des niveaux d’énergie que nous ap- 
pelons « associés ». Le caractère accessoire de numérotation de 1 apparaît si on 
compare (5.3) et (5.13): (n, 1, t) > N. 


Construction de la fonction d’état ‘*Y. — Les fonctions d’onde approximatives 
du systéme nucléons-mésons ont la forme: 


P? = p(X, n): |p) + p(X) ID, 
où |p> et |g) sont des vecteurs mésoniques de parité déterminée (p= + et 


q=+); YW” a la parité p. 
Les fonctions g* sont égales a: 


G(X, n) =3{(X,mM+9o(-X,m}, w(X,n) = 4 {o(X, n) — o(—X, n}, 
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ou p(X, N) = yi(%1 8141)... X,(x,8,t,); les xa sont les fonctions de quatre compo- 
santes, deux de spins et deux d’iso-spin, solutions de (3.1); notons-les 4, (Xa8ata)- 
La fonction ‘7? doit étre complètement antisymétrique relativement 4 l’échange 
simultané des variables d’espace, de spin et d’iso-spin de deux quelconques 
nucléons. Cela impose que g* et-q le soient également, et elles le seront auto- i 
matiquement si on antisymétrise g. On redéfini done: 


g(X,n) = antisym. (4, sia Ga : 


Si les y, sont orthonormées: (y,, x) = 6 
ment la norme 1: 


la fonction Y” a automatique- 


vy’) 


(P,P?) = <p |p>(¢*, o*) + <a> gd) = x [ov Mont 
n 


. Car Vantisymétrisation ne change rien, et <p|p) = <q|g>=1 par hypothése. 


6. - Calcul de la source du champ V,. 


Introduisons dans o* de (2.7) les solutions (5.3), connues aux fonctions 
radiales près. Les sources étant additives, commencons par calculer le 0° dù 
à une couche complète de nucléons (N, J): 


+9 
(6.1) Ons Ss / > Vago" Xx) = 
M ma-s 


1 1 
270 i Da fatuo SV (~ frat): Py(9, 9) + = frstel(”)* Pyarre(9, ) f= 


+ conjugué complexe} i 


L’opérateur o-V a une action simple sur les fonctions (1/r) f(r) D,,,:(0,9); 
on démontre que: 


Ò 1 2 2_ &2 | 
(6.2) oVao{o Ly —L ale 
dou: 
1 1 df J+4 l =J—te. 
(63) EVI Ould, (I fe) Camels 9s oy, 


en utilisant (5.4). 


pee 
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En appliquant (6.3), on a: 


d J 
de + ua: *) tamed > Dial Dim. + 


M 


1 d Jt#4 
na ri lora & pace a i fan pS Dr uire Prive + 
| a) fun E Dm Praise 


= 


oli : 
=a TE Fiam (È nni 5) feste” 2. Di ite‘ Om) + con). compl . 


Les sommes sur M sont des constantes: 


+4 jf se 1 
DI DÌ vm Pyar Ga > DF née ae 
(6.4) M=-J =-J 
y Din D JMWH8 — 0. 
M=—J 
Il ne reste ainsi que des fonctions radiales: 
È J+314d ; 
(6.5) Ong Di Tar — (fralm)te Fusite (M) + fr) fal?) . 


La source 0* produite par une couche complète est purement radiale; elle 
le reste naturellement si on additionne les contributions de N couches. Cette 
propriété est indispensable car nous avons postulé que V,= V,(r). Les couches 
incomplètes engendrent un 0* qui n’est plus de symétrie sphérique, mais elles 
sont en nombre très faible vis-à-vis des autres, et nous travaillons dans la | 
suite en négligeant leur faible apport. ou, ce qui est équivalent, en remplissant 
complètement toutes les couches. Ainsi, si 7 est indice des couches occupées: 


(6.6) or) = Ylt) 


Propriété: Les sources dues aux protons, d’une part, et aux neutrons, 
(autre part, s’ajoutent avec le méme signe. 
Dans une théorie scalaire, les sources produites par protons et neutrons 


: > x 
Se soustraient; si :=( "|: 
Xn 


o == xt = XX, = source de mésons scalaires neutres . 


Dans le cas pseudo-scalaire, il se produit le contraire; nous le démontrons 
pour les mésons 7° (neutres), en supposant qu’il n’y a pas de mésons nt 
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(chargés), pour simplifier. Le cas général est traité implicitement a la Sec- 
tion 7. 
Soient: V,=V,=0. L’éq. (3.1) devient: 


(6.7) (x, sh Lo-vre) = Ey. 


Xr 3 , 3 
Comme c2-( , on a respectivement pour yz, et x, les equations: 
Se AGH 


ti 


(6.8) (1, +4o-v,) Xo = EY , 


(6.9) (a, - Lo-vr.) (n= By, « 


On a vu que si une fonction y(x) satisfait (6.8), alors y(— x) satisfait (6.9), 


4 cause du caractère pseudo-scalaire de V(r). 


D’ou, réciproquement:: y,(x) = 7,(—x). Ainsi: 


oo + E Vizier. = a Vi~3(= x)o xo(— X)} = 03(— x) - 


Pour une couche complete: 0? = 08(r), Vou 0° = dì. 
0,70» 0p = n 


7. - Caleul d’un exemple. 


Les potentiels pseudo-scalaires V,(r) apparaissant dans (3.1) ne sont pas 


_ arbitraires, mais doivent satisfaire l’éq. (3.2) dont la source dépend des solu- 


tions de (3.1). Les sources construites avec les fonctions propres de l’hamil- 
tonien H,=— (1/2m) A+ U(r) sont identiquement  nulles, Pow SV, = 06 
on lui impose les conditions aux limites relatives au probleme. C’est la solu- 
tion triviale, et il faut s’assurer qu’elle n’est pas la seule. En d’autres termes: 
postulons dans (3.1) un certain V®; avec les solutions de cette équation on 
forme 0%, et calcule V satisfaisant (3.2) avec cette source. Puis, on introduit 
V® dans (3.1) et recommence l’opération. Il en résulte une suite V®; si cotter, 
suite converge vers une limite V,, et si V, 40, c’est une autre solution, celle 
qui a un intérét physique. Dans ce cas, on dit que la théorie est coherente. 

Le premier point è vérifier est existence générale de cette cohérence; le 
deuxième, que le V,(r) trouvé conduit è des « splittings » des niveaux J=1+4+% 
d’ordre de grandeur correct. On peut le faire directement sans recourir au 
procédé ci-dessus, en simplifiant le problème sans perdre en généralité et cal- 


 culant directement la solution V,. 


Rewer 
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Grace a l’approximation de premier ordre du calcul de perturbation, le 
calcul se décompose en trois parties semblables, une pour chaque composante 
V, («=1, 2, 3). Supposons un instant que nous tenons compte de tous les 
ordres du calcul de perturbation. Les sources dépendent alors des f, (6'8"), etc. — 
Pour se restreindre au premier ordre seulement, on ne conserve que les termes _ 


linéaires en 6. On obtient: 


fia kh a oy Boat | 
Dar ==> — non r? dr dr (0n1l7) DA. (Ban + Prin )In' ite r)), | 
î 
fJ+ild Be = . 
(7.1) Des “i Togni dr (Qui) 22 i SI nin 1S! "be T ‘aa . 
fJ+4id 
Cuts — dr ps Gp mlt) & (Bata + Bt) get) 


sources produites par une couche (n, J, 1=J—(e/2),t=1) de protons; pour 
les neutrons on a: 


PISS Sea A (ae } 
Oi wise Set n ora 72 dr (Guilt) 2 ( SE Pa Ban) Gn'itelr)) ’ | 
j 
[IS aed 
(7.2) Oped a 7 on = 72 dr (Gril) » A ae = = ae ) Inter) ’ 
i SEL LE | 
ia ve z E r Gp (Male digg rs i ie *)Gn't+e(7)) 2 : 


En vertu de (5.12): 


(pura Dia (04 p23) = = fonctionelle de V;(r), 


(BR + py) = pes (ses = fonctionelle de V,(r) , 


ip) pat ec) = fonctionelle de V,(r) (*). 


n' n nn n' n 
Posons: 


: a 
hg (pa na = 5 (Gat Lu ’ nia Sel nt Ba) . 


ee EE 


(*) Ce résultat montre que les sources 0* dues aux protons et aux neutrons sont 
identiques, non seulement pour les mésons r° comme nous I’ avons démontré a la Sect. 6 
mais aussi pour les mésons chargés. 
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Les &, dépendent de J et J, naturellement; ils valent: 


ij 1 | : 

Tg; (In) Inte) Valr) . 
E” esi Ei 

had + 1 J 


(7.3) Enin (1, J) = — 
Solutions de (3.1): 


Fak Li Ein 0 
Aninst=1(%) = (5) 7 Intl?) Pineal 9, gy) + D { 0 sli aa sp 


0 iL 
su | 2. ) 7 Intel”) Prrencs9, ~) . 
Equation pour V,(r): 


ER RD : 
(7.4) (A ir ia I i: a pe on yp dr (Guilt) 2 Errnlly J) Gn ite(”)) = 


Les fonctions y sont linéaires en é*, et la source de l’équation régissant V,, 
ne contient que des é*. Dans l’approximation du premier ordre le probléme 
se découple, et on peut traiter séparément les calculs relatifs aux trois poten- 
tiels V,_,,,- La théorie sera cohérente s'il existe des &,, 70, calculés 
par (7.3) et tels que (7.4) soit satisfaite. On voit que les calculs sont parfai- 
tement les mémes pour x=1,2 ou 3 (en conséquence, nous sous-entendrons 
dorénavant l’indice «). De plus les effets dans le splitting sont additifs; il 
suffit done de multiplier par trois le AE” trouvé dans un calcul pour avoir le 
splitting total. 

Le membre de droite de (7.4) contient une infinité de £%,; on peut en 
limiter le nombre, et rendre ainsi le calcul praticable. Les £,, sont des élé- 
ments de matrice (7.3) qui décroissent rapidement lorsque croit l’énergie sépa- 
rant les états entre lesquels ils sont caleulés. Nous n’en conservons qu’un 
(A une exception près), celui pour lequel | £°,,,, — £%,|= minimum, n’ variant 
(dans la numérotation habituelle, c’est pour n’= n). Il s’ensuit une com- 
pensation considérable dans les sources; considérons les g engendrés par les 
couches (n,l,J) et (n,l+e, J), dites associées (*): 


J+4 ed 
Orig dr È "e Enn(l, J) dr (GniGnite) è 
ER a 
Onites ee -! on : Enn(l + E, J) a dr (GniGnite) . 


(*) Voir remarque, Sect. 5. 
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On voit d’aprés (7.3) que é,,(l,J)=—énn(l+e, J), donc la somme des 
deux o est nulle dans approximation faite. Les seules contributions qui restent 
dans le 0,,, sont celles des couches près du niveau de Fermi qui sont associées 
à des niveaux non perturbés vides. 

La conséquence est qu’on peut ramener à un nombre très restreint le nombre 
des coefficients &,,, à déterminer. L’approximation que nous utilisons pour 
cela est assez radicale, mais pour calculer un ordre de grandeur de « splitting » 
et vérifier l’existence d’une solution, elle est suffisante. 


Données de l’exremple. — Nous prenons un noyau constitué de A = 120 
nucléons (N= 70 et Z= 50); le potentiel U(r) introduit est celui de l’oscil- —. 
lateur harmonique, qui permet de trouver une solution exacte 


A a ()- | 


U,= profondeur du potentiel = 59 MeV; 


a =rayon du noyau = 7.76:10- cm. 


Les fonctions radiales g,,(7) du système non perturbé doivent satisfaire: 


2 5 l(l 2 


| 2m dr? 2m x? 


La solution est bien connue: 


(7.5) Gril? j= E Pm"? exp (2/21 È pr (ati exp [—2]) , 
avec 
k=5(n Sees Ne e=(p si, yn = V2T(k+1+4 3)I(%+1)] 
Ex,=- U (1 =< ms "| s pt = 2ma?U, 


En tenant compte de la position qu’auraient les niveaux dans un potentiel 
en forme de cuvette et de celle qu’ils auront après la perturbation, le schéma 
des niveaux remplis s’établit comme indiqué a la Fig. 1. 

Les sources engendrées par les niveaux associés (reliés par une ligne sur — 
la figure) se compensent en première APRITE RIE Parmi ceux qui restent $ 
(indiqués 7 s) nous négligeons ceux de spin } et 3 car 0--J+4}. Nous retenons } 
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— finalement la source: 


(7.6) Q(7) = 2° 054,81) + Osos(1) + 043151). 
E EA 
Î neutrons protons 
N=6+ 8 ----- - --=---- N26) ------  -._---- 
5 aa aah mae 3) 5 b------ sssssas nana 
4 serene È cen TA era TAI ARA ae 
2 = 2 SS 
| 1 i uo 
ali ———_— I Sto dra A EA 1 L 1 salts 
<0 1 2 3 4 <0 1 È 3 4 


| Fig. 1. - Schéma de remplissage des niveaux non perturbés. Les couches notées a, b. c. 
sont retenues pour le calcul des sources. 


Fonctions d’onde des particules participant aux sources: 


Niveau a (n= 5, 1=4, J=%): 


JRE 1 
%aM)= + L Gs,a(") Paar 8(0, gy) + È n Yss(1) Ds, u,8(0, 9). 


Niveau b (n= 5, 1=2, J= 


Lolor 


ye 


— 
~] 
=] 

— 


al ol 
X(M) = vi È Y5.2° Dy u,5 i 7 (— E2903 + É3943)" Ds ws . 


Niveau c (n=4,1=3,J=1): 


lle yal Il 
Ye(M) == N, È Ys,3° Ds u 3 = È, È I5a° Ding) . 


avec: N,=V1+(6)?, N= V1+(&)?+(&)?, N.=V14+(f)*, les fonctions 
n’étant pas automatiquement normalisées car on a coupé dans le développe- 
ment en série. Pour simplifier la notation, nous avons posé: 


(7.8) È, = ACE 3) ’ E,=— E, 5(2, 3) ’ E,= Bake 3) ’ È, = &,4(3, 3) ® 


Le &,, suivant &,;(4, 3), par exemple, est é,,; il est proportionnel è 

(E°, — EB? ,)1=1(E},—;,), done au moins trois fois plus petit que SR 

Sa contribution au AF est i de celle de É,, on peut donc le négliger en 
première approximation. 

Par contre, & et £, sont du méme ordre du grandeur puisqu’il y a un niveau 

| (n=6, 1=3) en dessus et un niveau (n= 4, /=3) en dessous de db (n= 5, 

_ 1=2), tous deux à la méme distance. Il faut conserver les deux coefficients. 
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En vertu de ce qui précède ((7.4) et suite) nous ne différencions pas expli- 


a ee ARLEI 


i 
citement neutrons et protons, et sous-entendons dans les y les symboles | ) 
0 0 
ou (') car de toute fagon les sources sont identiques. En introduisant (7.7) | 


dans (7.6), on a: 


: p pla i 
(7.9) o(r) = Hey, > ae 496,5) = La (943954) + i 
r° dr 
Bi sia 14 si 
dieta Vara di (95.294,83) + Ys — PP de Lu). fi 
ou n par 
Y. a È, Di Y È, 
NA 0 IA (Gs)? + (E) 
CT, RANE 
9 14 (&)+(&)?’ “1+ (&)? 
Résolution. — Il faut trouver les &; tels que: 
1d dV 
7.10 — ji? omens (pant tt 
(7.10) CEE ie a >) [eV = o(r) 
avec le o(r) donné en (7.9), et, selon (7.3): 
(I) palla È 7, fard 
u2 
G(r) = 9659549 G,(7) = 913950» G3(7) = 9639509 G(r) = Is49a3 . { 


? 
A) Caleul de V(r) en fonction des &. Posons x=p(r/a), A=pa/p, È 
et changeons de variables dans (7.10); il vient: : 


=: 
ia 
En portant dans (7.9) les formes explicites des g,,, introduisant la nouvelle 
variable x et groupant les termes de méme puissance on obtient: i 


3 


(7.12) fay = —4 (BY (È vI+inti)z gy @ exp (2) + 


u\a 


4 sd 3 È 24 1 
+ gn >. = 16¥s) 5 


DA a 189 ib! P 
st (- Da YY. 7S V3 Ys) 5 di (a? exp [— 2}, 
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14) 9 E Tai 7 9 \-3 
=o (5151) , n=2(3) , Va = Va = (35) z 


L’équation différentielle est linéaire, ce qui entraîne que 


où 


(7.13) TATO 
È % a La vip ’ 
avec v, satisfaisant: 
| alia, ase i 
(7.14) er (e) — Ay, = — ae (a"*” exp [—#]). 


L’éq. (7.14), en posant v,= (1/x)w,, se ramène & 
” d 742v 
Ww,(@) — 22 w,(v2) = — — =" (a’*” exp[—@2]), 


qui se résout par la méthode de Green. En exigeant que v, soit continu è l’ori- 
gine ainsi que sa dérivée, et que 0, +0 (woo), on a: 


(7.15) ve) = = {exe [— Aw] dti È (t7+2” exp [— #?]) sinh At + 


0 


@ 


+ sinh Ax | dt 


® 


1d ; 
10 (rrexp[— #)fexp [- in| , 


Les intégrales donnant v, sont des cas particuliers de: 


5 
M(E, n, %, m) = | dt(nt”-? — 2t") exp [— (t+2)?], 


n 


qui s’intégre par partie, 4 une fonction d’erreur près. Mais le plus simple est 
de poser: 
: 
M(b, Ny % n) a |P.(c, t) exp FE (t + a)?]|5 = ba(a) fat exp i (t as a)? ] ’ 


n 


avec P,(x,t)= (polynòme en « et t)= > A®(«)t™. 


8 - Supplemento al Nuovo Cimento. 
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En dérivant M, on a: 


DÀ {(m a 1) Agta == 2aAm = 2A, att" — bn, (a) — nin? — in : 


dot on tire une relation de récurence pour les A%,(%). 

Cette méthode a l’avantage d’étre applicable 4 d’autres exemples sans chan- 
gement. 

On trouve finalement: 


Il 
(1.16) = eso =—37440(22, 2) exp [— 2°] + 


T 1 
+ LF bran(412) È (exp [Ae Ble + 4/2) — 11+ exp [- 20](P(e— 2/2) + 1) 
où 
v+2 
Gotan A[2 0) = Y Ages (A/2) at, 
m=0 
A‘ (1/2) = coefficients dépendant de A= ua/p, calculés è partir de la rela- 


amen 
tion de récurence, 


xt Al2 


9 : 
E(a + 4/2) = n | exp [—@]dt. (JAHNKE et EMpE, fonction d’erreur) , 


0 


(voir Fig. 2). 


Fig. 2. - Fonctions »,(x), » = 0, 1, 2 (formule (7.16)). 
Et V(r) explicitement en fonction des ¢;, par l’intérmédiaire des Y; (7.9): 
(7.17 Vines = i È 7 
17) (= FEY 0 (Mola =3 ti BY, ‘0,(0) + 


3 4 4 21 189 
ca (nr. =" oe v3 V3 + lai SÉ (- os YaYa + da Za) ee} . 


1100 


115 


MODELE SEMI-PHÉNOMÉNOLOGIQUE DU NOYAU ETC 


En portant dans la définition des € Vexpression 


i B) Calcul des €,. 
de V(r) calculée ci-dessus, on obtient quatre équations algébriques pour les 


quatre €; inconnus: 


da) f pe TI fan f p° I y! FINI 
E = u20; drG,(r) V'(1) i 20>. dG) Vie 
0 0 
at ee caren î dG, 
(J aU, ye D 048) {dn 1,(0) ss 
0 
Posons 
È dG, 
71 =| deva) - 
(7.18) | Ley) d 
0 
Il vient au total: 
43 2 
ESRI AI 


. 9 (ef È = 
(7.19) E=f 22a, "È, 


Toutes les grandeurs figurant dans ces équations sont connues; on calcule 
graphiquement les Ji, et laisse pour l’instant le choix de f ouvert. En posant 


“p 10 a i 
SE aU, nye li); Tr) 
i io é Ope ) Cay fp? 3 | 24 ; | 189 ) 
#2 = 72.0200, 2aU, Dee PEREGO : ada Jo a di dora ® 
(7.19) s’écrit: 
4 
(7.19’) Bo ad x, We) s 
vont 


Les « sont des coefficients numériques, et les Y,(€) les grandeurs définies 
(7.19') est une équation non linéaire 4 quatre inconnues € . 


eM (09). 
La solution est, pour f=1.1 (*) 
£0312) “e = 0050, E 06080; “2S 0.285 , 
(7.20) dou: 
O Y,= 0.284, Y:=0.034,. Y:= 0.464, Y,= 0.264. 


(*) Moyenne des valeurs généralement admises pour f 


gi 
1101 


116° J.-P. AMIET 


C) Forme explicite de V(r). En remplagant les Y, dans (7.17) par i 
leur valeur, on trouve: ; 


(7.21) V(r) = 3.82-0,(2) — 19.2+0,(@) + 42.7-0,(a) (MeV). 


| (Voir Fig. 3). 


1 i 
di PEA e “ago 
0 1 i 

i) 


Fig. 3. — Potentiel pseudo-scalaire V(r), r=x-a/p (formule (7.21)). 
Explicitement: 


Va) y= FE lot) (x) exp [— #7] + 
+ BÈ (exp [Ae] [Mle + 4/2) 1] + exp [- Z0]{E(e— 2/2) + 1)} 
Q(0) = 2.15 + 0.8807 + 0.26! — 0.013a° + 0.2505, 
C= p=, p= 3.65 A= 146., B= 465. 
Discussion de la cohérence. — La condition physique de cohérence est: 
(7.22) VE Wa 


où V est un potentiel postulé, et V le nouveau potentiel, calculé è partir des 
sources dépendant de V. Explicitons cette relation: 


2 2 
> tote) = Dion, e) 
v=0 v=0 


D’ot 


(7.23) 3, =.0,(E) = 6,(E). 
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Ces trois relations relient les anciens È, (i=1,...,4) et les nouveaux È , 
qui d’autre part sont définis par: 


(7.24) a > of Y,(é) 


k=1 


les xi étant les coefficients de l’éq. (7.19). 

(7.23) est évidemment satisfaite si on pose E,=&,, ce qui, avec (7.24), 
donne l’éq. (7.19’). Par contre, si on substitue (7.24) dans le membre de gauche 
de (7.23), on a trois relations pour les quatres inconnues £;: le problème a lair 
indéterminé. Ce n’est pourtant pas le cas car dét|ai|= 0, ce qui entraîne 
une dépendence linéaire des É,, ou des È,. Il ne reste alors que trois inconnues 
dans l’éq. (7.23); (7.22) et (7.19') sont bien des conditions équivalentes. La 
dépendence linéaire des £, est évidente si on considère (7.11): les quatre pro- 
duits G,(r) ne sont pas linéairement indépendants. Cela provient de ce que 
deux fonctions g,,(r) et g,,,(r) ne sont orthogonales pour n # n' que si J=U’. 
Il y a exactement autant de €,,,, linéairement indépendants dans le cas général 
qu'il y a de produits g,,,,(")"Yn:(") linéairement indépendants. Ainsi, dans tout 
probléme analogue au présent, on aura une dépendance linéaire entre les £,,,, 
à condition de prendre garde en approximant de ne pas négliger trop de coef- 
ficients. Ce fait important permet de tirer quelques conclusions sur l’existence 
d’une solution cohérente dans un cas quelconque. 

Partons d’un cas à deux coefficients: 


2 
oa, dét|a,|#0, i=1, 2. 
k=1 


On veut È;= €,; la solution se trouve a l’intersection des courbes &;= > agli 
k 


et &= Y a:Y,. Ces courbes ne se coupent qu’a Vorigine, sauf cas très parti- 
k 
culiers (Fig. 4). On a done presque toujours, dans ce cas, uniquement 


la solution triviale; il en va de méme si 
120, > 2. 6, 
Prenons maintenant <= 1, 2, 3: 


3 
(7.25) COST. 
Ka 
avec dét|x{|= 0, qui entraîne 
0 1 é, 
(7.26) &, = 0,8, + bb... Fig. 4. 
$ 9 - Supplemento al Nuovo Cimento. 
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On débute avec: 


2 
&= @%Y,+ 2 %Y,, 


k=1 


2 
di iy ant, = > Oils ’ 


k=1 


et résout en gardant Y, arbitraire (Fig. 5). 
Le lieu des solutions, noté en pointillé, va de 0 à une valeur < 1 (’) | 


lorsque Y, varie de 0 a 4, qui est sa valeur maximum (Fig. 6). La solution | 


6, a 
1 Tal 
2 
a5), | | 
I a i 
; È L 53 | 
Fig. 6. i 


de (7.27) se trouve à l’intersection des courbes fi(é3) = di &1(&3) + b,°&(&) et 
fo(€s) = & (Fig. 7). | 
L 


e domaine des solutions est maintenant beaucoup plus large. Il en existe 1 
une dès que df./dé3|;,-, > 1. . 
(Remarque: Le signe des a joue | 
peu de réle, car toutes les expres- 
sions ont une symétrie centrale.) : 
Il en va de méme lorsqu’on a | 
plus de trois coefficients, mais avec — 

des hypersurfaces. 

%, En résumé, nous voyons que la 
cohérence du probleme n’existe pas 
par coincidence, et pour une valeur — 
déterminée de la constante de cou- | 

plage f. Les quelques considérations faites ci-dessus sur les systèmes d’équat- 

ions intervenant montrent qu’en général l’existence d’une solution n’est pas 
du tout problématique, comme ce serait le cas si dét|ax{|0. On peut voir 
de méme que les solutions tendent vers zero si f >0. | 


(*) Les ai dépassent rarement 1 en valeur absolue, et seulement de peu. 
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Splittings. — Il n’est pas intéressant de calculer les splittings de tous les 
niveaux dans le présent exemple. Le potentiel U(r) choisi est trop loin des 
potentiels réels en forme de cuvette, et il donne lieu à des fonctions radiales 
9ni(r) trop ramassées vers l’origine. D’autre part, nous avons négligé des sources 
provenant des niveaux à petits 1, étant donné leur moindre importance quan- 
titative. On obtient ainsi un V(r) de pente plus forte vers l’origine, ce qui 
favorise les splittings pour / petit. Le fait important est que les splittings ont 
un ordre de grandeur satisfaisant; on trouve, par exemple, pour n= 5, 1= 4, 
J=4+41: 

AE=1.5 MeV. 


En multipliant par trois, pour tenir compte des contributions des trois 
types de mésons, on arrive a: 


AE, = 4.5 Mev . 


Le calcul de ce AF fournit une justification d’une des approximations 
| que nous avons faites. 


AE = AH,,.— AE,,, = différence des déplacements . 
Par exemple: 


AE, <a? > [Enm(4, 3) LA a7 Esa) ni (2? Si (077-2 A (Sgr ee piani = 
10,81 — 180/06 0,6 Mev. 


(Ep) P1021 ,)*5 (Exe)? 0a 


Les termes négligés pour le calcul de V(r), comme &,,, & ., sont nettement 
plus petits que les termes retenus, et particulicrement ceux contenant au dé- 
nominateur une différence |H°,,,,—H?,| plus grande que 9 MeV (qui est le 
minimum). 


8. — Niveaux d’énergie. 


Du point de vue du spectre d’énergie, le présent modéle pseudo-scalaire 
est semblable au modéle en couches avec couplage spin-orbite. On peut prendre 
le méme potentiel scalaire U(r) dans les deux cas, ce qui donne la méme ré- 
partition des niveaux. A part cela, le terme o:VV,7* lève la dégénérescence 
en J comme le terme L-S, de facon que le niveau J=/+4 se trouve en- 


1105 


120 J.-P. . AMIET 


dessous de celui 4 J=1—4; en effet, les déplacements des niveaux valant: 


(8.1) AG ease = >> iS. (0, J) |? ( oi E Ei) = Lenn |? (Brits = E ’ 


n'n 
7 
n 
x=1,2,3 


et étant donné que: 


(Ea E a) >0 ’ a Ea) 0 ’ 
on a AP joi 0, Ale Vea 


L’indice n” désigne les deux niveaux (n’, 1 +1) les plus proches du niveau 
(n, 1) avant perturbation (avec un U(r) en forme de cuvette, on a n°"= n). 

Les È, Sont suffisamment plus petits que €,,, pour que le résultat ne 
soit pas modifié si on les fait intervenir. L’exemple calculé à la Section 7 le 
montre, et dans un cas où le potentiel U(r) défavorise l’approximation en 
rendant équidistants les niveaux non perturbés. L’ordre de grandeur du split- 
ting que nous obtenons est de quelques MeV en moyenne, ce qui correspond 
aux valeurs observées. La question qui reste ouverte est la variation du split- 
ting en fonction des nombres quantiques n et Z. Il est très difficile d’y ré- 

pondre en général, car la variation des AH 

à travers du spectre d’énergie dépend con- 
--  sidérablement du potentiel U(r). Il est 
tout au plus possible d’avancer une remar- 
que spéculative. 

Considérons les niveaux les plus proches 
de la limite de Fermi, de différents J, nu- 
mérotés (a), (b), (c), etc., dans la Fig. 8. 
Les déplacements ont la forme (8.1). Pour 


Fig. 8. — La ligne pointillée dési- (a), toutes les différences H?.,,,— EH}, ont le 
gne la limite de Fermi. méme signe car aucun niveau (n’, 1-+1) ne se 


trouve en dessous de lui. Pour (b), (c), ete., 
on a des termes de signes différents, d’ot compensation d’autant plus grande 
qu'il y a de niveaux en dessous du niveau considéré. Et de méme pour 
(/—1). Tout dépend de l’importance des termes ¢,,,_,, par rapport à la va- 
riation de &,,,, en fonction de J, mais il est probable qu’on obtient en général 
une diminution du splitting avec / au voisinage de la ligne de Fermi. Et ce 
sont la les seuls splittings mesurables. 


9. — Parité. 


Etant donné une solution p(X, 7) = gt(X, n) + q-(X, n) de Véq. (3.1'), Pétat 
fondamental est de parité (+) et est représenté par Yt+=gt|+>+q-|—>, 
lorsque AH > 0; il est de parité (—) et représenté par Y= y+|—)+q@-|+4), 
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quand AZ<0. Nous allons montrer que la parité ainsi définie correspond è 
la parité donnée par le modéle en couches ordinaire. 

Le AF est donné par la formule (2.8). Il est démontré à la fin du para- 
graphe que le second terme est nul; il ne reste que: 


(9.1) AE = h{(g*, ¢*) —(@ @)}. 


Le facteur h ne concerne que les mésons, et (y+, y+) — (p-, g-) seulement 
les nucléons. Calculons d’abord l’ordre de grandeur des produits scalaires 
N== (y+, p+). Les fonctions g* sont des sommes de produits de fonctions 
relatives 4 une particule, du type: 


ik I 
Jom = 3 Inu(?) Piu (9, g) + Enia n In'ivet Dive - 


Les &,,, sont largement plus petits que 1, comme nous l’avons vu dans 


|. Pexemple calculé 4 la Section 7. Pour simplifier, nous les prenons tous égaux, 


au plus grand par exemple; cela ne change pas le signe de la différence 
(N*— N) qui seul importe. En posant v= (nlMJt), v'= (n, 1+e, M,J,t), et 
indiquant la parité de é,(x)= (1/r)g,,,y, par un indice + ou —, on a: 


Pour / pair CAI be Xp oy og 
l impair Vries) Ly Ge ele : 
La norme des È, est égale à 1. 


Remplissons un niveau de (nlJt) déterminés, et appelons g*(k) les fonctions 
correspondant a % particules sur le niveau: 


a) l pair: 
1 particule: gl) eye eles 
N+(1)~1, NS ces 
récurence: Orel) = Co ak) Rec. (kh): 
P(k +1) > ESTA + Oe (k). 
D’ou 


N+(k+1) = N*k) 21,  Nk+1) SN) Le. 


Done, pour une couche de | pair: N+ ~1, N7= &, quel que soitk<27+1. 
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b) l impair: 


1 particule: GUI, PL) Crs 
N+(1) = &, N-(1) 21; 
récurence: g*(k+1)= EE Vt (k) + Ca PIE) ; | 
dk +1) 2 EPM + EG, (hk). 
D’ou: | 
N*(k 4-1) = WN (k), Nk AVA) Ti 
et: | 


Nitya) e AC) 53° 815680 
Nk) ot ~~ N*(k) DE Resta pate 


Pour un niveau plein, N+~1 et N-~é& car 2J+1 est pair. 

La généralisation 4 un coeur formé de couches complètes plus une couche 
incomplète est immédiate; le coeur ne joue pas de ròle puisque dans son cas 
(N*+— N-) ~1—&?> 0; seule importe la parité du / de la couche incomplète, 
et de son occupation k. Si l’on suppose maintenant que h> 0 et utilise les 
résultats ci-dessus, on en déduit le tableau comparatif suivant, des parités 
calculées dans les modèles en couches ordinaire et pseudoscalaire: 


Modèle Su couches Modèle pseudo-scalaire 
ordinaire 

Toutes les couches | se AE=h1-&2)>0, état fon | 
complètes en damental de parité — + 1 | 
Ceur +k particules sur | oo AE> 0, 
niveau a ! pair [et Pao parité = + 1 î 
Coeur + % particules sur | he AE ~ h(1 — £?2)(—1)*, 
niveau à 1 impai PTO eee ité = (= 

| impair parité = (—1) 


La parité obtenue dans le modéle pseudo-scalaire est bien la méme que 
dans le cas du modéle en couches habituel, sous réserve que h> 0, ce que 
nous allons démontrer. 


Lorsqu’on ne suppose rien sur les normes individuelles de g* et |+), on a: 


n= CHI) _ IHIH) 
CIT GT | 
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Jusqwici nous avons laissé les kets | +) indéterminés, et travaillé directe- 
ment avec <M|W,|M) qui est donné par une équation aux sources (2.4), les 
fonctions g* étant inconnues. Le probléme est maintenant inversé, et les | +) 
sont les inconnues a calculer. 

En introduisant 7 = o(X, n): |M) dans (1.2), multipliant chaque membre 


par g*(X, n) et intégrant, on obtient: 


MODÈLE SEMI-PHÉNOMÉNOLOGIQUE DU NOYAU ETC. 123 


Nous prenons pour les vecteurs | +) et |—> la meme approximation que 
pour g* et g, et posons: |+)= partie paire de |M), |—>)= partie impaire 
de |M). Nous sommes en présence d’un problème résoluble exactement: la 
description des mésons virtuels engendrés par A nucléons dont le mouvement 
est connu et donné par les fonctions y+, calculées comme 4 la Section 7. 

Le terme d’interaction peut étre transformé en intégrant par parties: 


A 
Hy= (p(a, n) H'q(e, n)) = 2 fato bs (x), Ta Xa(x)V: V(x) = 
a asl 


wai 
Lu i > [d'eV(x)W(ztor%) - 


a=1 


Posons: 
1 AI 
| A ae Oe V= V2 (V,+ iV.) (*), 
(9.3) ‘ KG: 
| o*(x) = Dea, o(x) => a (este) 
a=1 a=1 a/2 
Il vient: 
(9.4). H, =— 1 [rato + oV+ o*V*). 
(22 


En introduisant dans (9.4) les formes explicites de V; et V (*), puis les 
transformées de 0° et o relativement 4 la base Wim: 


(9.5) im) = | a0? V inlky x), Vm(k) = | dao V imlk, X) ; 


Wim(Kk) =|d't0* Vim (k, X) ; 


(*) Voir Appendice. 
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on trouve: 


, OEMS; 2 Of dk 
(9.6) H, st Z Dil /9 {a;(k, Lg m) Uim(k) == a3 Wim = A:Vim ae Ai Vim == 
Be -  1=0} 2% 
—l<mSl 0 i ae AyWim "n ati} © 
D’autre part: 
H, = D dkmxaX(k, 1, ma,(k, l,m) . 
a 0 


L’hamiltonien H,, = H° +H), peut étre diagonalisé au moyen de la trans- | 
formation unitaire U= exp[iS], où 


foo} 


x [atogt a;ju— afu* + av — at0* + aj;w — afw*) . 
sm 


i 
V2 Di fa. | 


On vérifie aisément que: i 


[S, H,|}=—iH,,  [S, H,]=—72e, [S, ]=0, 8 
ou 
= 3 n) dal ATE v*ot w*w) . | 
2 LU Lm Ox 
0 
En conséquence: 
(9.7) UH, U = H® —e. | 


Et le vecteur d’état solution de (9.2) est: 


(9.8) |M) = exp [i8]|0) = DIE S"[0). 


UA 


En supposant, comme nous avons fait jusqu’ici, que 0*= 0*(r), les seules 
transformées (9.5) non nulles sont celles d’indices 1=0, et m=0. |M) ne 
contient plus que les mésons S, et ses parties paire et impaire sont: 


foo} 42” (ce) q2°+1 


(9.9) Co ee 


S2"(0 —\ = —_ i 
Zoni e Rin 


x 
» 
ae 


S?"|0) est une somme de vecteurs contenant un nombre pair de mésons, et 
S?"*1/0>) un nombre impair. | 
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En vertu de l’unitarité de U et de l’orthogonalité de |+> et |—>, on a: 
Gri SM Sa | 
D’autre part: Be 


Fee; 


car 
<M|H%,|M)= <0|U-+H%,U|0)=e, et <—|H]+)=0. 
D’ou: 
{|Hx|=>= se 
9.10 ra 
di Si cy 
avec 
= 1 
| ON _ 1\N S2w+1 3 
si r= 2, | a ete AS oi (Qn +1)! (Qn’+ 0)!’ 
; 12 1 
<— | Hy |—> o i ad DE 


2 la (Qn + 1)!(20/+ DI 
O <0 | 8241 Hi, Sensi a. 82*+1H%, S2r+1 |0) à 


Le calcul de <—|H%,|—) se ramène & celui des éléments de matrice 
<0|.8?"|0). En vertu de [S, H°]=— iH,, et [S8", H,,]=—i2meS”™ on a: 


ALE Senti LL UE Ae Sa ae ar ET San SE 
+ S2 H® Gente + 2(2n — 2n'+ 1)e 82" . 
On démontre alors aisément que: 
ERE e da Pra be S2nt1 — Santa's) ET + 
4 H Satta L 920, 1 )(2n Desens. 


Comme H°|0)=0, on a: <0|H3,8%++0|0) = <0|8%*++bH,|0)= 0; et 
finalement: 


_|H3|—) = e ¥ (—1)* |S” |0)- PESARE] 
CHA) = 2 BU OLS: È, Gata) 


Pour simplifier, calculons <0|S?”|0> dans le cas où les mésons chargés sont 
absents; il reste S= (i/V2) (f[u) [Ak wx°(a(k) — a* (Ke)) w(k) (en supprimant les 
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indices), 0° = 03(r) entrainant Win(k) = dwOmou(k). La généralisation au cas 
général est quasi evidente. 


<0| 827 |0)= (a) (2) [en OTe (:;)—a*(ke,)) | 0) w(k) (kan) 


IT x +++ Orgy) 
Posons: 
a,=a(k,), DECIO : 
a’ == < = 15 ’ 2N 
a,=a*(k), v=1, 


On peut alors écrire: 


SITI (ak) — a*@))| = (ayn Olay... a5] 0. 


Vi... , VaN=0 


Les relations de commutation se résument a: 
[ai, a; | = (v —v')dok —k.). 


Tous les termes se terminant à droite par a}, = @,y donnent zéro puisque 
a(k)|0)= 0. De plus, tous les termes tels que »,+...+7,,4 N donnent éga- 
lement zéro car ils contiennent des nombres inégaux d’opérateurs de création 
et d’annihilation. On calcule la contribution de chaque terme restant en ap- 
pliquant les règles de commutation; on prend le premier opérateur a; à partir 
de la droite, et on le déplace jusqu’a ce qu’il soit devant le ket |0>; à chaque 
commutation avec un a” on obtient un nouveau terme contenant une fonction 
o(k;—k,) et deux opérateurs de moins. On procéde de méne jusqu’a ce qu’on 
n’ait plus que des termes donnant zéro et des produits de N fonctions è. 
Jamais deux fonctions des deux mémes arguments n’apparaissent dans un 
méme produit car chaque a’ n’apparait qu’une fois dans un terme. De plus 
aucun des termes de départ ne peut s’obtenir d’un autre par permutation des 
opérateurs. Tous les produits de 6 seront différents. En outre l’intégrale mul- 
tiple est symétrique par rapport a toutes les variables k;. Le résultat est done 
forcément: 


(9.11) 1° S ROIO I 


perm ky ... kon 


Si p, est le nombre de termes de la somme (9.11) pour N, et Py, pour 
N—1, on a la relation: p, = (2N —1)p,_,; @’ou: 


Dy = Q2N —1)!! | (N > 2). 
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Finalement, chaque terme donnant la méme contribution, on a: 


eo 


ard N 
(9.12) <0 | 82" |0> = (2N — 1)!! (2) | a 101) 


0 


Si on prend les mésons chargés: 


r 2710) =. (€ pp 
(9.12) <0|S2|0) = (2N (A/S 


-{u2(k) + v(k) v*(k) + row} |= (2N —1)!!-1¥ 
avec v(k) = voo(k), wk) = Wolk). * 


Substituons ces valeurs dans les développements donnant <—|H},|—) et 


“GS 
2 il 
BRIO = (+ Zen pin. > ani): 


N=0 nina (20)! (2n 


e 1 
ea => (-I(2N4+ II > | 


N=0 n+n'=N 2n * 1) : (2n'+ 12 


En vertu de 


1 
vis On)! On)! — veya On + Dim + DI 2 ON! 


i Lala 

NONE 11.7%. oP gots è eS ae wee I 
Se sO (2n)!(2n')! 2 N! 

(—1)**1 


1 LÌ 
Sr — (QT)¥+1 
SR (2n +1)! (2n’+1)! 3 (N+1)! | ) 


(21)*, 


_ 1)(2N + 1)! IMA. 
( ( 


C'est, au 3 près, le N°" terme du développement en série de exp[— 27]. 
On a ainsi, simplement: 


(9.14) <—|H®|—) = }e(1 + exp([—27]), <+|Ht|+)=4e(1—exp[-27]), 
(9.15) <—|—> = }(1-exp[-21]), <+|+> =#(1+ exp[--21]). 
Et: 
4 exp [— 21] 
(9.16) =e 1— exp [— 4J] . 


Conclusion: h> 0, puisque e> 0 et [> 0. 
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Evaluation de h. — Caleulons d’abord È dans le cas ‘'o*= 0*(r). 

D’après la définition (9.5), v(k) = w* )=fa%x(01 (r) +i0*(r r))V oo(k, x), car 
Vd =F kre Gs p) = (42)-* kj, (kr). ie posant: u,(k )= fas 0% (1) oo(k, x) 
on a: 


(9.17) oe Al ASS ta an 


2 
‘ode + bs 


L’intégrale se calcule facilement, et on trouve finalement: 


3 
(9.18) pe 2a È farro) VES (o3)ie 
a=1 


A titre d’exemple, calculons e avec les valeurs trouvées à la Section 7; 
introduisons dans (9.18) le o(r) donné en (7.7), en prenant garde au fait qu’il 
contient déjà f/u. L’intégration est immédiate en vertu de (7.9). On trouve: 


(I 


12 
Sg oa (ONE: — 3Y.6,+3Y,6, + 4V 6) = 118 MeV. 


j 
: 
î 


Puisque nous calculons h uniquement pour avoir une idée de son ordre de 
grandeur, nous pouvons nous permettre d’évaluer I au moyen d’une approxi- 
mation grossière. La transformée w,(k) est vraisemblablement déjà très petite | 
quand k= wu (0x= 2u); aussi posons-nous: 


(k? + ua) ~ wk? is Ta: 


De sorte que: 


Et 7 


4 exp [— 2e/u] = 
= ==) 3) Mev. . 
£DE exp [— 4e/u] È 


Le AE est plus petit que h A cause du facteur |N*—N-|<1. Pour ne 
ne pas avancer de nombre sans justification, disons que sa valeur se mesurera 
en dizaines de MeV. 


Quelques valeurs limite. 


1) limI=lime=0, limh=y, lim+|+)=1, lim—|—)=0. 
f>0 f>0 f>0 


f>0 F>0 
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Si le couplage s’évanouit (constante de couplage f +0), le vecteur |M} 
tend vers |0) qui est de parité (+). Mais les niveaux (+) et (—) ne se 
réunissent pas car h >, de méme sue AE; la fonction g(X, 7) tend vers une 

fonction a parité déterminée. 


2) lim J =lime=oo, limh=0, lim <+|+>=lim «|> =}. 
fro fo 


fro foo foo 
Si le couplage croît fortement, le nombre des mésons a tendance a s’équi- 


librer entre les vecteurs |+ > et |—>, et les états (+) et (—) tendent a se 
rapprocher: 


AE({=0)=135 MeV, AZ(f=11)<80 MeV, AE(f=00)=0. 


Démonstration de (y-, H'g*) — (pt, Hip-) = 0. 


il 
roe [Val )»_ € IV.]+)), 


i oc (— DER 
ergy an genti o 
(9.19) <4 |Wx |> I Guten 1) yy S18 Vv, 82" 10) 
è f %,(k) 
(9.20) [S, a,(k, ti, m)] = [S, at(k, l, m)| = Ono mo: .. ws ’ 


où S est défini en (9.12). 
Soit 


D’ou: 
[sep = mE) 


Ce qui permet de dire: 
4<0| S2°W, S2"'+2|0> = (2n'+ 1)K,(r)<0| 82+"? |0> = fonction réelle; 


de méme pour W*. En conclusion, ¢<+|V,|—> est réel, puisque V, est une 
- combinaison linéaire des W, et W%, et Vi(r) = H,=0. 
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10. — Utilisation des fonctions d’onde; moments magnétiques. 


Il existe trois types d’opérateurs dont on peut avoir besoin de calculer la 
valeur moyenne pour un état Y du noyau. 


1) Opérateurs Oy n’agissant que sur la partie nucléonique de la fonction 
donde ¥?= ot |p)+ 9" |@ (p=+, a= F). 
CP? | Oy |P?) = <p|p}(p*, Ove") + (a|D (9, Ove) ; 


car Oy p(x, n)-|M)=|M}Oy; de plus, on peut choisir <p|p> = <q|@ =1, 
dou: 


(10.1) <¥?|O~|'P "e (X, 1) Ov gt(X, n)+ PAX, n) Ov~-(X, n). 


n 


En particulier, du point de vue du modéle en couches, on peut s’intéresser 
a un opérateur Of n’agissant que sur la fonction y, du nucléon d’état 
(MalaMatata); en laissant de còté l’antisymétrisation, on a: 


Pay Fe Coa is REG Be ’ 


ou YG, sont les fonctions combinées relatives au A—1 nucléons res- 
tants, et x, les parties paire et impaire de y,. Avec une normalisation adé- 


quate des fonctions, (10.1) donne dans ce cas: 


(10.2) | OY 


y= farts; Xa + Xa Ov Xa) - 


On remarque qu’il n’est pas nécessaire de connaitre les kets |p) et |q>. 


2) Opérateurs O,, n’agissant que sur la partie mésonique de YW. 


<P|O, |P) = (vt, vt) <p|O,|p> + (97, 7 )<g|O,|@ - 


3) Opérateurs O agissant sur les deux parties. 


(P| O|P> = (pt, Qt) + (97, 20) + (pt, 207) + (97, 2*o*), 


où 2°= <p|O|p, Q°= q|O|m, 2=<p|0|9. 


de ce type. 
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Les opérateurs en relation avec la parité et la charge, en particulier, sont 


| 


- ware. fy + 
eee 
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Moments magnétiques. — La formule (10.2) permet une constatation intéres- 
sante, relative au moment magnétique d’un noyau calculé dans le modeéle 
pseudo-scalaire. Si l’on suppose que le moment magnétique est entièrement 
dù aux nucléons, on tire facilement de (10.2) que la contribution des couches 
pleines est nulle. Ainsi, le modéle prévoit que les noyaux magiques ont un 
moment magnétique nul. Pour un noyau constitué d’un coeur fermé plus un 
nucléon: 


DE | Gacy | paral = ee, be) 


xi et x, sont fonctions propres de J’, et de L® pour les valeurs propres différant 
de 1: = J,+}. Ainsi, avec le modéle en couches ordinaire on aurait pour 
<p> soit ut=|x}|®(xi, wx}); soit wo =| yz |" (xa, #20), suivant que l= J, +3 
est pair ou impair. Avec le modéle pseudo-scalaire on a: 


delie ye (lyal?+ |x|? =1)- 


Le moment magnétique se trouve entre les lignes de Schmidt. 


En terminant, je tiens è remercier le Professeur K. BLEULER qui m’a sug- 
géré le sujet de ce travail. Ses directives éclairées et ses conseils me furent 
constamment d’une aide précieuse. 


APPENDICE 


A) Unités et symboles. 


Les unités utilisées sont telles que 7#= = 1. 


constante de couplage; 

masse d’un méson 7; 

— masse d’un nucléon; 

A=Z4+N= nombre de nucléons; 

X= (x, ...; x,) = ensemble des coordonnées spatiales des A nucléons; 
n = (8t,, ---, Sat.) =ensemble des coordonnées de spin et iso-spin; 

x = indice d’iso-spin. Toute sommation sur x est sous-entendue. 


Soo 
Il 


* désigne le conjugué hermitien; 


+ désignent la parité lorsqu’ils sont utilisés comme indice; 
~ est réservé à distinguer deux grandeurs différentes, désignées par la méme 
lettre. 
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Les opérateurs agissant sur l’espace de Hilbert nucléons-mésons x sont 
désignés par des lettres italiques grasses: H, V(x), II,... Les grandeurs se 
transformant comme des vecteurs de l’espace euclidien à trois dimensions sont 
notés en lettres grasses: x, V, o,... Les éléments de matrice et produits her- 
mitiens sont parfois notés de deux différentes maniéres suivant la forme expli- 
cite de l’espace de Hilbert: 


<M,|H|Mx), on (91, Hy.) 
<M,|M>, (913 Pa) + 


B) Champs d’opérateurs mésoniques V(x). 


Les opérateurs V, sont exprimés en représentation de Schrédinger, dans 
un système de coordonnées sphériques (x= (r, 0, g)): 


(ce) 


dk l 
= > Iz (a;(k, l, m): Vim(k, x) aî(k, l, m) SELE x)) ’ 
Ss 40x 


1=0 
—lSmSl 0 


foe) 


V(x) = (V,+ iV.) 


o a o (a, ( (k, l, m)* Him k, x) + at (2, k, m)* Wim ( k, x)) i 


Les fonctions de base ¥,,,(k, x) sont, explicitement: 
Vim (hy x) = kji(kr) Yim (9, @) , 


où j,(kr) = V(kir) J a(kr), J i43(2) = fonction de Bessel (JAHNKE et EMDE), 
Y,,(9, p)= fonction sphérique normalisée, Op=V ke + pw? =energie d’un méson 7, 
1, m= nombres quantiques orbital et azimutal d’un méson x. 

On a les relations d’orthogonalité: 


[ove Vk, x) Wm(k', x) = O11 Omm' Ò(k = k') ’ 


>) dk Vié(k, x) Vin(k, x') = 6(x — x). 


—lSmS1l 0 


Les opérateurs de création et d’annihilation a%(k, l,m) et a,(k, l,m) satisfont 
les relations de commutation: 


[a,(k, l, m), at: (E, U', m')] = bag'011'Omm' O(K — k’) , 
[a,(k, 1, m), a} (k’, l',m)}=0. 


Les vecteurs mésoniques |M} sont des combinaisons linéaires de vecteurs 
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_ |..-Mikim +> è nombres d’occupation fixes: 


az(k, l, m)a,(k, lg m) [es Me to = hei aS . ‘ 


C) Parité II. 
L’opérateur de parité II relatif au système nucléons-mésons a la propriété: 


(II, H] cas [II, Hy| a [II, Hy] = (II, H'| =0. 


II se décompose en un produit d’opérateurs lorsqu’il agit sur un produit de 
| fonctions des variables respectivement nucléoniques et mésoniques 


On a: 
I=ILgylIly,= 11,0, , NY = IIyg(X,n)I,{|M}, 


4 ° Oo 
DY U +L) Myym 


a 11,9(X, 7)=o(— X, n), Thy |. Mim --- = (— 1) ata 


eee artes 


Propriété: [Ily, H]40, [IL,, H']#0, bien que [ILyII,, H]=0. Pour sim- 
| plifier, IL, est noté x lorsqu’il agit sur une fonction concernant un seul nucléon: 


my(x) = y(— x). 


D) Quelques proprié!és utilisées des éléments de matrice de V(x). 


1) Soient |m, +) et |n, —) deux vecteurs quelconques mais de parité 
déterminée. Alors: 


<4,m|V,(x)|n, —>= fa(x)=  fa(-x) fonction paire de x. 
‘+, m|V,(x)|m, +> = ga(x) = — gal x) fonction impaire , 


<—, n |V.(x)|n, —> = ha(x) = — hal(- x) fonction impaire . 


| 2) Soient |.S,>, |S.>, deux vecteurs relatifs a des mésons S exclusivement, 
i c'est-à-dire que: a,(k, 1, m)|S,)=0 si 1#+0. Alors: 


| 48, |Velx) | Se) = fall xl) - 


3) Soient |+>, |—, deux vecteurs de parité déterminée et relatifs a 
des mésons S. Les propriétés 1) et 2) entrainent: 


<+1V,.(x)|—> = Val); <+|V.(x)|4>= 0. 


4) Si <m|V.(x)|n) = V(r), on a: (x, o-VV,(r)] #0. 
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